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THE USE OF INFRARED SPECTRA FOR THE DETERMINATION 
OF ABSOLUTE MAGNITUDES 


P. C. KeENAN* AND J. A. HYNEKT 
Washington, D.C. 
Received February 14, 19-45 


ABSTRACT 


On infrared spectra taken on a scale of 48 A/mm at the Perkins Observatory it has been found pos- 
sible to apply the method of spectroscopic absolute magnitudes to stars of type M. Lines of neutral and 
ionized elements serve as criteria which allow the assignment of luminosity classes Ia (very brilliant 
supergiants), Ib (ordinary supergiants), II (fainter supergiants), III (giants), and V (dwarfs). The most 
sensitive feature is the blended line at \ 8514, which is due, at least in part, to a member of a low- 
lying multiplet of neutral iron. 

The luminosity criteria should be useful chiefly in the study of the atmospheres of the cooler compo- 
nents of spectroscopic binaries in which there is a marked difference in temperature between the two 
members. 


The pioneer exploration of stellar spectra in the near infrared region extending from 
\ 7000 to X 8700 was carried out in 1934 by Merrill,! who charted the principal features 
found in stars of types B8-M6. Since that time the speed and quality of photographic 
plates sensitive in the infrared region have been greatly improved; in particular, the 
Eastman I-N emulsion now has consistently high sensitivity between A 7500 and A 8600. 
With these plates and efficient grating spectrographs, it is possible to obtain infrared 
spectrograms of late-type stars with exposure times comparable to those needed for the 
ordinary blue region. Thus the long-wave part of the stellar spectrum is now available 
for practical use in the comparison of atmospheric characteristics of red stars. 

In order to determine the usefulness of this region for the estimation of spectral types 
and absolute magnitudes, spectrograms of a number of standard stars were obtained by 
the second writer with the infrared spectrograph on the 69-inch reflector of the Perkins 
Observatory. The spectrograph employs a grating of 15,000 lines per inch and, as used, 
was fitted with a 250-mm camera, the combination giving a scale of about 48 A/mm near 
\ 8000. Most of the spectrograms were taken on Eastman I-N plates, although a few of 
the brighter stars were photographed also on Eastman IV-N plates in order to take ad- 
vantage of the higher contrast and finer grain of this emulsion. 


* On leave from the Yerkes Observatory for research in connection with the war. 


t+ On leave from the Perkins Observatory to the Applied Physics Laboratory, Johns Hopkins Uni- 
versity. 


1 Ap. J., 79, 183, 1934. 
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The stars used as standards are listed below in Table 1. In addition to spectral types* 
and luminosity classes,’ estimated absolute magnitudes are included; but it must be 
kept in mind that for the more luminous stars these are quite uncertain.‘ 

Comparison of the spectra of these reference stars showed that the infrared bands of 
TiO behave consistently with the blue bands and can accordingly be used to assign spec- 
tral types through the range MO-MO6. The most useful bands are the groups with heads 
at AA 7054, 7088, and 7126 for types MO-M3 and the weaker bands at AA 7743, 7820, 
8441, and 8450 for types M4-M7. These bands, which arise from or near the normal 
level of the 77 o molecule, were originally found in the spectra of M-type stars by Mer- 
rill’ and by Bobrovnikoff.6 The atomic lines from low levels of titanium, etc., also in- 


TABLE 1 


STANDARD STARS OBSERVED IN THE SPECTRAL 
REGION AA 6900-8700 A 


Star Type Luminosity M, 

ee KO III 0 
Ki Ilr — 0.2 

K2 III 0 

K3 Ib-II — 3 

K3 II 2 

K3 III 0 

RE K5 Ib — 4 
0-Ia — 6.5 

MO III 0 

Mi Ib — 4 
M2 Ta — 5.5 

M2 III 0 
BD+36°2147........ M2 Vv +10.5 
M3 Til — 0.5 

M5 III — 1 


crease in strength with diminishing temperature; but they are not so sensitive and 
consistent as the bands. 

The problem of estimating luminosities can best be approached by considering first 
the early M-type stars, among which there are known reference stars ranging in luminos- 
ity from dwarfs of the main sequence to the brightest supergiants. In Plate XXVI the 
spectra of a set of these stars have been arranged in order of decreasing absolute bright- 
ness from top to bottom. In the spectral region represented here the bands of atmospheric 


+ aeaaatia Keenan, and Kellman, An Ailas of Stellar Spectra, Chicago: University of Chicago Press, 
43. 


3 For the M-type stars the system is that of Mount Wilson (Adams et al., Ap. J., 81, 187, 1935). 


4 The luminosity scale for the M-type stars discussed by Keenan (Ap. J., 95, 461, 1942) has been 
modified slightly on the basis of work on the M-type variables at Mount Wilson (Joy, Ap. J., 96, 344, 
1942, and Wilson, Ap. J., 96, 371, 1942) and on the modulus of the double cluster in Perseus (Bidel- 
man, A p. J., 98, 61, 1943, and Morgan, unpublished). 


5 Bur. Standards, Sci. Papers, 14, 487, 1918. 6 Ap. J. 78, 211, 1933. 
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PLATE XXVI 


RW Cep 

M Cep 

5 Lac 

A Peg 

BD + 36° 2147 


— Na 8183 

Na 8195 

Ti 8426 

Ti 8435 
468 

— 8514 * 
— Fe 
— Fe 
Fe 8689 * 
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Call 8498 
Call 8542 
Ca Il 8662 


LUMINOSITY SEQUENCE FOR M-TypeE STARS IN THE INFRARED 


Lines marked with an asterisk (*) are enhanced in the supergiants 


(Note: What is indicated as atmospheric band O, in the region \ 8200 to \ 8300 is due to water 
vapor and not to oxygen.) 
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oxygen and water vapor, indicated in the reproduction, are prominent and hide some 
of the stellar features. Between them, however, the line spectrum of the stars is practi- 
cally undisturbed. 

Conspicuous in the dwarfs and supergiants are two pairs of lines of neutral potassium 
and sodium, arising, respectively, from levels of low excitation potential (Table 2). Both 
pairs lie in the neighborhood of atmospheric bands with which they are partially blended, 
as indicated in the table. Nevertheless, it can be clearly seen in the reproduction that 
they pass through a minimum in the giants. 

A number of lines of somewhat higher excitation potential, including the strong Ti 1 
lines AA 8364, 8378, 8382-3, 8397, 8412, 8435, and 8548, are markedly stronger in giants 
than in dwarfs and share the further strengthening which most lines exhibit in passing 


TABLE 2 
STRONG LINES OF LOW EXCITATION POTENTIAL 


7664. 907 
7698. 979} K 0.00 42S — 42Pe 7665.41 (Atm) 
8183.258) = 8183.1 (Atm) 
8194. 836/ Na 2.09 32D 8193.7 (Atm) 


TABLE 3 
INFRARED LINES WHICH STRENGTHEN WITH INCREASING LUMINOSITY 


Ip Cep ALab Source Transition 
$468.2. .... 3ns 8468 . 41 68.5, 20,300 
8662.4..... 15s Cau 2D 1.7 
8688.8..... 5s 8688 . 63 150 y 


from giants to supergiants. The Fe lines \\ 8388, 8582, and 8675 also belong in this 
group. 

Most important for the separation of supergiants from the other stars are a group of 
lines which are strongly enhanced in passing from low to high luminosities. These lines, 
several of which are listed in Table 3 and are marked with asterisks in Plate XXVI, arise 
chiefly from rather low levels of neutral iron or of ions of metals of low ionization poten- 
tial. The latter are rarer in the red region than in the blue part of the spectrum, since the 
first excited levels of these ions generally lie several volts above the normal level. There 
are, however, a few which are favorably situated for observation. With the exception of 
the members of the infrared Ca triplet, they are rather weak in stellar spectra and are 
usually blended with other lines on spectra of moderate dispersion. Consequently, the 
identifications in Table 3 are tentative only. The assignment of the atoms chiefly re- 
sponsible for each line has been made on the basis of the ‘“Multiplet Table” of Charlotte 
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268 P. C. KEENAN AND J. A. HYNEK 
M. Sitterly,’ with regard also to the elements found by D. M. Davis to be represented 
by lines in the blue region of the spectrum of a Sco. 

The essential point is that these observed lines, though frequently blends of several 
lines from different elements, strengthen progressively with increasing luminosity of the 
star. The strongest and most sensitive is \ 8514, which is indicated on Plate XXVI.° In 
dwarfs and giants \ 8514 is only slightly stronger than \ 8518 of 771, the line adjacent 
to the red. In passing to stars of higher luminosity, the 771 line remains at about the 
same intensity, while \ 8514 strengthens rapidly until in RW Cep it surpasses \ 8435, a 

‘blend of two of the strongest 77 1 lines. The tentative prediction of Merrill! that the Ca 1m 
triplet “‘may have an interesting relation to luminosity” is fully borne out, as a glance at 
the reproduction will show. Most of the remaining lines in Table 3 appear only in the 
brighter supergiants (on spectrograms having a scale of about 50 A/mm) and are of 


TABLE 4 
INFRARED CRITERIA OF ABSOLUTE MAGNITUDE (TYPES Mo-M,) 
BEHAVIOR 
SENSITIVE COMPARISON 
LINES LINES 
Dwarfs Giants Supergiants 
(K \ 7664.9) 
Strong Weak Strong <K 7699.0/ 
|Na 8183 Atm 8177 
|Fe 8327 Atm 8321 
Very weak Weak Rather strong | Fe-Ti 8468 Ti 8435 
8514 Ti 8518, 8435 
Fe 8689 Fe 8675 


value chiefly in distinguishing the very brilliant supergiants of luminosity class Ia (RW 
Cep, » Cep) from those of class Ib or II (a Ori, 5 Lac). 

The enhanced lines are consistently stronger in RW Cep than in » Cep. This result 
confirms the conclusion indicated by the lines in the photographic region‘ that, as far 
as the spectrographic evidence is concerned, RW Cep is more luminous than wu Cep and 
is one of the intrinsically brightest stars known. 

A summary of the luminosity criteria is given in Table 4, where the lines have been 
grouped according to their behavior. Numerical estimates of relative intensities have 
not been given, for they will depend upon the dispersion and contrast of the spectro- 
grams. The criteria have been described on the assumption that direct comparisons will 
always be made with a set of reference stars such as those of Plate XXVI. 


7 We are greatly indebted to Mrs. Sitterly for her kindness in comparing our approximate wave lengths 
with infrared lines listed in her revised but yet unpublished ‘“‘Multiplet Table.” 


® Ap. J., 87, 335, 1938; 89, 41, 1939. 


* It is worthy of note that three of the lines which show the strongest enhancement in the supergiants 
are \\ 8468, 8514, and 8600, which are all members of an Fe1 multiplet with an excitation potential of 
2.2 volts, not greatly different from that of other Fe lines which do not exhibit such marked enhance- 
ment. The great strength of \ 8514 in RW Cep suggests that the Fe line is blended with a line from an 
unidentified ion. ~ 
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It is important to keep in mind the limitations in the accuracy of spectroscopic abso- 
lute magnitudes. In the infrared the precision will be slightly less than in the photo- 
graphic region, chiefly because of the smaller number of sensitive lines. The “resolving 
power” in luminosity estimated from plates of moderate dispersion can be set at about 
2 mag. This means that we cannot expect to distinguish reliably between stars 
differing in absolute magnitude by less than this amount. What can be accomplished, 
however, is the separation of dwarfs, giants, and supergiants and the division of the 
supergiants into those that are moderately bright (M, = —3) and those of highest lumi- 
nosity (My = —6). 

One of the most useful fields of application of the infrared luminosity criteria is the 
study of the red components of spectroscopic binaries in which the blue region of the 
spectrum is dominated or disturbed by the light of the hotter component. A notable ex- 
ample of such double stars is VV Cep, the infrared spectrum of which is described in the 
paper which follows. 


é 
a 
a 
d 
hs a 
of 
Bey 


THE PRESENCE OF STRONG LINES OF O1 IN THE 
INFRARED SPECTRUM OF VV CEPHEI 


J. A. HynEK* P. C, KEENANT 
Washington, D.C. 
Received February 23, 1945 


ABSTRACT 
The spectrum of the red component VV Cep was photographed in the region \\ 7200-8700 in Octo- 
ber and November, 1944. No emission lines were observed; the absorption lines of the metals had in- 
tensities normal for an M2 supergiant and indicated a spectroscopic absolute magnitude (M,) of —4 or—5. 
The only anomalous feature was the great strength of the blended absorption lines of Or at \ 7773 
and \ 8446. The intensity of these lines was as great as that normally observed in supergiants of types 
B-F, and physical considerations suggest that the lines are produced by oxygen atoms in the normal 


atmosphere of the M-type star. 
The presence of oxygen atoms in the excited states from which these transitions arise may be due to 


resonance absorption of the intense ultraviolet diluted radiation from the B-type star near \ 1300. The 
dimensions of the system are consistent with an energy density sufficient to populate these levels. 


In the ordinary photographic region the spectrum of VV Cep is produced by light from 
both the red and the blue components of the system, the observed radiations from the 
two being equal near \ 4000. It is only near primary minimum, at intervals of twenty 
years, that the pure spectrum of the supergiant M-type component can be observed in 
this region. For this reason, it is important to turn to the infrared spectrum in order to 
study the atmosphere of the red star throughout the orbital cycle. 

In October and November of 1944 two spectrograms of VV Cep, covering the region 


from A 7200 to A 8700, were taken with the 69-inch reflector of the Perkins Observatory . 


by J. A. Hynek and N. T. Bobrovnikoff, respectively. One was taken on an Eastman 
I-N plate, the other on a IV-N plate; apart from contrast effects, the spectra appear iden- 
tical. The I-N spectrogram of October 7, 1944, is reproduced in Plate XX VII; for com- 
parison a spectrogram of a Ori is given, since the lines of neutral Fe, Ti, and the other 
metals were found to have nearly the same relative intensities in VV Cep as in early 
M-type supergiants. 

The spectral type is close to M2, in agreement with V. Goedicke’s estimates from the 
spectrum in the blue region.’ The appearance of the lines which were found, in our pre- 
ceding paper,” to be sensitive indicators of luminosity is consistent with a visual absolute 
magnitude of —4 to —5, since the enhanced lines were found to be slightly stronger than 
in a Ori or a Sco, (Ib), and slightly weaker than in wu Cep, (Ja). The closeness of the 
match with a Ori is apparent in the reproduction.* 

It will be immediately noticed, however, that there is one striking difference between 
the two spectra. The blended lines of neutral oxygen at \ 7773 and \ 8446, which are ab- 
sent (or, if present, barely perceptible) in a Ori, form two of the strongest features in the 
spectrum of VV Cep. In the latter they have approximately the same strength as in 
supergiants of types A and F, where they normally reach their maximum development.‘ 

* On leave from Perkins Observatory to the Applied Physics Laboratory, Johns Hopkins University. 

t On leave from the Yerkes Observatory for research in connection with the war. 

1 Pub. Mich. Obs., Vol. 8, No. 1, 1939. 2 Ap. J., 101, 265, 1945. 

3 The similarity of the blue regions of VV Cep and yu Cep was noticed by Goedicke at the time of pri- 
mary minimum in 1936 (op. cit.). 

4 The variation of the O 1 lines through the spectral sequence can be followed in the reproduction given 
by Merrill, A p. J., 79, 183, 1934, Pl. II. The lines are conspicuously strong in the spectra of a Cyg, 6 Ori, 
and ¢« Aur. They show a marked luminosity effect. 
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Since the O11 lines are weak or absent in all other late-type stars which have been 
photographed in the infrared, their strength in VV Cep must be associated with the 
proximity of the high-temperature (B-type) component of this binary. They cannot be 
part of the normal spectrum of the B-type star, since its radiation is negligibly small com- 
pared to that of the red star in this part of the spectrum. The alternative that the O1 
lines arise from the projection of the B-star’s ordinary atmosphere against the disk of the 
M-type star is improbable because of the low central intensities of these lines; absorption 
of this strength could be produced only if almost the entire disk of the red star were oc- 
culted by the atmosphere of the B-type star. This would be inconsistent with the gen- 
erally accepted picture of the system according to which the red star has a radius at least 
one hundred times greater than that of the blue component. Furthermore, the spectro- 
grams were taken not at secondary minimum but more than a year earlier, and no change 
was noticeable in the interval of seven weeks between exposures. 

There remain two possibilities: the excited O 1 atoms responsible for the lines are situ- 
ated either in the atmosphere of the M-type star or in a diffuse cloud surrounding one or 
both of the components of the binary system. The absence of even a trace of emission 
associated with these strong lines (these lines show emission in P Cyg, y Cass, and ¢ Tau) 
speaks against the second alternative. The evidence from radial velocities is inconclusive 


TABLE 1 
RADIAL VELOCITIES OF VV CEP FROM INFRARED LINES 
VeLocity (In Km/SEc) 
PLATE DATE l 
TiLines | Fe Lines | Ca Triplet | O 1 8446 
10- 7-1944 —29.4(7) —18.7 (4) | — 37.7 (3) —37.7 
Snr | 11-25-1944 —25.4 (8) — 20.8 (8) | —31.9 (3) —53.1 


at present, because at the time of observation the system was approaching secondary 
minimum. 

The velocity of the M-type star (predicted from the velocity-curve given by Goedicke) 
for the epoch of our observation is about — 15 km/sec; for the B component, about — 25 
km/sec. The measures of the present spectrograms of the red component of VV Cep 
are summarized in Table 1. 

The discrepancy between the observed velocities and those predicted for the red star 
are most readily explained by the known irregularities in its velocity-curve. The sys- 
tematic differences between the velocities of different elements are possibly an effect of 
stratification. 

On the assumption that the oxygen lines are produced in the atmosphere of the red 
star, we have the problem of explaining why they have intensities to be expected only 
in the spectrum of a much hotter star, while the lines of 77, Fe, etc., appear with strength 
normal for an M-type spectrum. 

The answer may lie in the combination of two circumstances: the intense radiation 
of the distant B-type star in the far ultraviolet and the differences between the oxygen 
and metallic atoms in their potentials of excitation and of ionization. The lower levels 
of the infrared lines of Fe and 77 have excitation potentials lying between 1 and 3 volts, 
and the relative distribution of the atoms among the excited states will be governed en- 
tirely by the radiation of the M-type star. The only effect of the B star may be expected 
to be some decrease in the intensities of all of these lines, since the numbers of neutral 
atoms will be decreased if there is increased ionization from radiation in the region 
AA 1500-1700. Actually, these lines from neutral metals appear to be very slightly weak- 
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ened, relative to the Ca 1 lines in VV Cep, as compared to a Ori or other similar super- 
giants (cf. Ti 1 \ 8435 with Ca 1 8498 in Pl. XXVII). 

The situation of the oxygen atoms is essentially different, as can be seen from the 
Grotrian diagram, of which the lower part is shown in Figure 1. The observed infrared 
lines originate from the 38S? and 3°S° levels, which are 9.11 and 9.48 volts, respectively, 
above the normal 2°P level of the neutral oxygen atom. Transitions up to the 3°S? level 
give rise to the ultraviolet triplet \A 1302-6. These are the ultimate lines of oxygen, since 
the other low-lying levels, 'D. and 'So, are metastable. 


Odd Even 


3 5, 3 5 


10.69 


9.46 3°S 
9.0 


A.17 


1.96 


Fic. 1.—Lower energy-levels of O1. Permitted transitions are shown by full lines; a few forbidden 
transitions are included and are shown by broken lines. 


If it is assumed that in the atmosphere of the M-type star the energy density of dilute 
radiation from the B-type star is sufficiently great near ) 1300, the population of the 
33S? level will be governed by this radiation quite independently of the radiation of the 
M-type star itself. When once excited to this state, the atoms will be ready to absorb 
light at wave-lengths given by the series 3°S°—m*P, of which the first member is the 
observed pair at \ 446, where the intensity of radiation from the red star is high. It 
will be observed that the other infrared triplet, AA 7772-5, is the first member of the series 
3°S°— m5P, Since there are no quintet levels lying below 35S, this level can be reached 
from below only by the transition 2*P — 3°S°. Such intercombination lines can be expect- 
ed to occur in the spectrum of O 1,5 and this transition has been observed as the pair of 


6 Condon and Shortley, The Theory of Atomic Spectra, p. 218, Cambridge University Press, 1935. 
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lines at \ 1355.7 and \ 1358.7, having emission intensities approximately one-half those 
of the permitted triplet lines. In VV Cep the relative intensity of the \ 7774 and \ 8446 
blends is very nearly the same as in a normal early-type supergiant, such as 6 Ori. 
To test the reasonableness of this interpretation of the observed spectrum, a rough 
calculation of the radiation fields can be made on the basis of the orbit derived by 
Goedicke! from spectroscopic and photometric data. Admittedly, the dimensions of the 
system are quite uncertain, but the general agreement of our independent estimate of 
the luminosity of the M-type component with the values given by Goedicke suggests 
that the order of magnitude of the results will not be far wrong. Out of the several pos- 
sible alternative sets of values given by Goedicke for the relative dimensions and tem- 
' peratures (allowing for a diffuse atmosphere of the red component) the following were 
chosen as reasonably consistent with the spectroscopic luminosity results: 


18,000°K | 13.0R© 
3,000 | 1220 RO 


A = semimajor axis = 33 A.U. = 7095 RO 


Consider an atom in the atmosphere of the M-type star. It is assumed to be sufficiently 
close to the surface so that the relative dilution factor, W,, for radiation of the red star 
is 3. The computation is then independent of the radius of this star. Since the M-type 
star fills a considerable part of the orbit of the system, the atoms in different parts of its 
atmosphere will be at a wide range of distances from the B-type star. Assume, as a rea- 
sonable mean, that this distance is A — Ry = 6000RO. The ratio of the dilution fac- 
tors, W,/W,, for radiation from the two stars is then 2.7 X 10~*. In Table 2 the rela- 
tive contribution of the two stars to the field of radiation surrounding the atom is given 
for different wave lengths on the assumption that both stars radiate as black bodies. 

It is evident that the blue component dominates the radiation field in the far ultra- 
violet to such a degree that the picture will not be substantially altered by any reason- 
able changes in the dimensions of the system. The spectral region within which the light 
from the two stars approximately balances is, however, quite sensitive to the tempera- 
ture of the B-type star, which may be anywhere between 12 ,000° and 30,000°. For the 
crude model assumed above, the intensities from the two stars approximately match 
near A 3100, and the radiation of the M- type star is relatively negligible at its own sur- 
face for all wave lengths shorter than about \ 2000. 

The computed flux of dilute radiation at \ 1300 surrounding the oxygen atoms is 
2.9 K 10" ergs/cm?/sec/cm of wave length. The same energy density of radiation 
would be found at the surface of a star radiating as a black body at a temperature of 
slightly less than 5600° K. If the ultraviolet radiation of the stars were actually given by 
the Planck curve for the star’s effective temperature and if the lines at \ 1300 were the 
only transitions determining the population of the 3°S level, the intensities of these 
infrared lines should be the same as in an F7 supergiant, for which the effective tem- 
perature is in the neighborhood of 5800°. The close agreement with the intensities of 
these lines and those in the spectrum of the FO supergiant e Aur is therefore of consider- 
able interest. The observed intensities, however, also match those in considerably hot- 
ter stars; this is not surprising, since neither of the above assumptions holds true exactly. 
The essential conclusion to be drawn from this crude calculation is that the absolute in- 
tensity of ultraviolet light from the B-type star is of the right order of magnitude to 
give rise to observable lines. 


6 Hopfield, Ap. J., 59, 114, 1924. 
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If the assumption examined in the previous paragraphs proves correct, a more detailed 
treatment of the problem will be important, for the relative intensities of the lines of 
different excitation potentials provide a set of conditions relating 73, 7, Rg, and A, in 
addition to the data available to Gaposchkin and Goedicke in their original solutions for 
the dimensions of the system. 

The first need is for a thorough study of the entire red and infrared spectrum of VV 
Cep in order to compare intensities of lines from as many elements as possible. Among 
the important lines are the members of the Paschen series of hydrogen. The lines from 
P13 to the limit of the series fall within the range of the available plates; but other strong 
absorption features, including the three lines of the Ca 11 triplet, fall at nearly the same 
wave lengths and prevent the definite identification of the Paschen series in VV Cep. 
The appearance of the spectrum indicates that they cannot be markedly stronger than 
in other late-type stars and are not present with intensities comparable to those found 
in B Ori, a Cyg, or other hot stars. The lack of strengthening in VV Cep is not necessarily 
surprising but requires explanation. 


TABLE 2 


RATIOS OF THE RADIATION FLUXES FROM THE BLUE AND RED 
COMPONENTS AT A POINT IN THE ATMOSPHERE 
OF THE RED STAR 


| | 
| 5.6X10? || 3000A............ 
2.71073 


Nitrogen is somewhat analogous to oxygen in its arrangement of energy-levels. The 
infrared lines, the strongest of which lie at AA 8680-86 and 8703-18, originate from a 
‘Pe level which is reached from the normal 4S level by transitions in the neighborhood 
of } 1200. These lines, however, are less easy to observe than the infrared lines of O 1, 
for the number of permitted transitions from the N 1 ‘P° level is much greater and none 
of the lines is comparable in intensity to the strongest oxygen lines.’ None of the V1 
lines has been detected with certainty on the Perkins spectrograms of VV Cep. The 
question of whether they show any strengthening in VV Cep, as compared to other late- 
type stars, can be settled only with spectrograms of higher dispersion. 

The interpretations of the strength of the O 1 lines as a result of direct excitation by 
the radiation of the B-type star leads to the prediction that the oxygen lines should dis- 
appear completely during the total phase of primary eclipse (occurring next in 1957), 
since at that time none of the atoms on the observed side of the M-type star is exposed 
to direct illumination from the hotter component. Whatever residual excitation occurs 
would be caused by diffused light. Actually, it is possible that the absorption will be re- 
placed by emission from the illuminated rim—resulting in a selective “flash’”’ spectrum. 

If this interpretation is correct, one might expect similar phenomena in related systems 
like ¢ Aur. A plate of this star taken on February 11, 1945, does not show the O 1 lines. At 
this epoch, however, the B-type star is approaching eclipse (two months prior to ingress) 
and is unfavorably placed for the excitation of the red star atmosphere that is in the 
light of sight to the earth. Infrared observations of both the VV Cep and ¢ Aur systems 
at other epochs should prove of considerable interest and will clarify the tentative inter- 
pretation of the O 1 lines presented here. 


7 The relative intensity of the Or and V1 lines can be compared in Plate II and Table 2 of Merrill’s 
paper, Ap. J., 79, 183, 1934. 
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AND SWASEY OBSERVATORY 
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Received March 6, 1945 


ABSTRACT 


In this paper the description of the mounting, optics, and objective prism of the 24-inch Schmidt-type 
telescope is presented, together with numerical data pertaining to it and its performance. 


The program of research and instruction in astronomy at the Warner and Swasey Ob- 
servatory of the Case School of Applied Science, in Cleveland, Ohio, was considerably 
advanced recently by the acquisition of a 24-36-inch Schmidt-type telescope. It is in- 
tended to present here the essential details and characteristics of this instrument. 

In 1919 Mr. Worcester R. Warner and Mr. Ambrose Swasey built and equipped an ob- 
servatory to serve the basic educational needs of the Case School in astronomy. The 
principal equipment consisted of a 93-inch refractor, two transit instruments, a zenith 
telescope, and the associated recording apparatus, such as precision clocks, chronographs, 
etc. Later a 3-inch Ross-lens astrographic camera was added. 

When an appeal was made in 1939 for funds to enlarge and expand the observatory, 
many friends were ready to help. Gifts totaling over a quarter of a million dollars were re- 
ceived, principally from persons closely associated with the original founders of the ob- 
servatory and of the Case School. The additional facilities include a Schmidt-type tele- 
scope, a new building to house the telescope, an auditorium and an exhibit hall in mem- 
ory of Mr. Warner, a library, offices, a photometer room, and a shop. The building was 
completed in October, 1940, and the telescope was put into operation during the month 
of November, 1941. The Warner and Swasey Company of Cleveland, Ohio, designed and 
constructed the mounting and the optical elements of the instrument. 

In future publications concerning observations made with this instrument we shall re- 
fer to it as the “Burrell telescope,” a designation in memory of the late chief engineer of 
the Warner and Swasey Company. 


MOUNTING 


The mounting is of the cross-axis type, along the same lines as the large reflectors built 
by this company (PI. XXVIII). 

Polar axis.—This axis is composed of three castings: a cubical center carrying the dec- 
lination axis, an upper section provided with radial bearing resting on the north pier- 
head, and a lower section having a radial and thrust bearing resting on the south pier. . 
The extended portion of the lower section is made of steel forgings pressed into place and 
carrying the right-ascension driving gear and circles. 

Declination axis.—This axis consists of a hollow shaft machined from a forging. It is 
carried on antifriction bearings of ample capacity for the radial and thrust loads. These 
bearings are mounted on the center section of the polar axis. 

Telescope tube.—The tube is made up of a number of welded-steel sections, ,*,;— 3%; inch 
in thickness. The diameter of the lower section is 4 feet, with the mirror cell bolted to the 
lower end. The upper section has a diameter of 2.8 feet and supports the correcting lens 
and a removable dewcap of 17 inches in length. Diaphragms, arranged to prevent inter- 
reflection, together with ribs extending the length of each section, add to the transverse 
rigidity of the tube. This construction assures rigid alignment not only between the opti- 


275 


pe 
. 
| 
aa 
— 
4 
: 


276 J. J. NASSAU 


cal parts and the plateholder inside the tube but also between them and the guide tele- 
scope on the outside. 

Mirror cell.—This cell has an over-all diameter of 44 inches. The bottom supports for 
the mirror consist of three spiders, each having three arms terminating in flat horizontal 
rings of 6 inches in diameter. The rings are so placed that they support equal portions of 
the weight of the mirror. The back of the mirror is smooth-ground to avoid excessive fric- 
tion between the ring pads and the glass. Four side supports provide means for collimat- 
ing. In order to facilitate handling the mirror, the cell is made in two parts, the inner part 
being removable to permit the attachment of a lifting band around the circumference of 
the mirror (Pl. X XIX). 

Lens cell.—This cell is mounted in a flanged cell provided with three sets of push-and- 
pull ‘screws for “squaring.”’ This cell is mounted in a ring provided with collimating 
screws; three bolts attach the ring to the tube. 

The shutters.—There are two shutters inside the tube of the telescope, one just above 
the mirror and the other just below the lens. Each is made of two cylindrical pieces of 
metal, which, when open, lie along the inner sides of the tube. One shutter is used prima- 
rily to protect the mirror when the telescope is in actual operation, and the other is to con- 
trol the exposure. Both are operated by hand wheels at the lower end of the tube, close to 
the guide telescope. 

A third shutter has been added for use in cases of short exposure where more accurate 
timing is necessary. It is a circular, occulting metal sheet 9 inches in diameter, placed 
just in front of the photographic plate. A lever device at the lower end of the tube actu- 
ates it. 

Plateholder mounting.—This mounting is supported by four thin vanes attached to 
the tube 4 feet below the correcting lens. The 2-foot arm which holds the adaptor for the 
photographic plateholder is pivoted 2} inches from the intersection of the vanes, which 
permits swinging the plateholder attachment toward the outside of the tube, as shown 
by the dotted arc line in Plate XXIX. This enables the observer to change plateholders 
and to alter the focus without reaching into the tube. Changes of focus are made by a 
screw-and-slide mechanism. A graduated scale and a large micrometer dial permit direct 
setting of the focus to 0.001 inch. A clamp-and-latch mechanism aligns the plate ac- 
curately in the focal plane after loading. 

Circular plateholder (P|. XXX).—This plateholder is made of aluminum and is de- 
signed to accommodate circular glass photographic plates. These flat plates are bent to 
conform to the spherical focal surface by the curved back plate of the plateholder, which 
is pressed firmly against the back of the photographic plate and clamped in place. The 
plateholder has proved highly satisfactory; there has been no breakage at this stage of 
handling the plates, although more than one thousand plates have been taken. Clips are 
provided for holding color filters directly in front of the photographic plate. 

Guide telescope.—This is a 93-inch refractor of 123 inches in focal length. The usual 
right-angle adaptor, eyepieces, and illuminated cross-wires are provided. 

Visual observing.—This type of observing may be done on occasion by interposing 
a 7-inch flat mirror in the convergent light-beam to bring the focal surface to the side of 
the tube in the manner of the usual Newtonian arrangement. This mirror is mounted in 
an arm which is ordinarily folded against the inner tube wall (this feature is not shown in 
Pl. XXIX). By turning a knob on the exterior of the tube the observer may throw the 
mirror into position for use. 

Motions in declination.—Four speeds in declination are provided. For guiding, the 
speed is 2 minutes of arc per minute of time, and for setting the guide star at the inter- 
section of the cross-wires of the guide telescope the speed is 15 minutes of arc per minute 
of time. These speeds are obtained by means of a two-speed reversible motor, driving 
through a gearbox having two trains of gears and a solenoid-operated clutch, arranged 
so that the same motor may operate either the slow or the quick motions, the quick mo- 
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PLATEHOLDER ATTACHMENT SHOWING PLATEHOLDER OUTSIDE OF TELESCOPE TUBE 
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PLATE XXXI 


PHOTOGRAPH OF A PORTION OF THE NORTH AMERICAN NEBULA (left) AND DENEB MADE WITH THE Bur- 
RELL TELESCOPE WITH AN EXposuRE OF 10 MINUTES 

The irregular double ellipse just above the nebula is the ghost of Deneb. The long streak is due to the 

reflection from one of the vanes supporting the plateholder. 
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tions being 46.0° and 7.6° per minute of time. A clamp with magnetic release connects the 
driving mechanism to the.declination axis. 

Motions in right ascension.—The sidereal rate of driving is obtained by the use of a ;',- 
h.p. synchronous motor. As a partial compensation for differential refraction, the rate is 
slower than the theoretical rate by 1 second of time per hour. A clamp with magnetic re- 
lease connects the worm wheel to the polar axis when the instrument is driven. The 60- 
cycle A.C. current provided for the synchronous motor is furnished by the Cleveland 
Electric Illuminating Company. Its frequency remains sufficiently’ constant that no 
guiding is necessary during a 10- or 15-minute exposure. 

To increase the efficiency of the telescope, particularly when it is used for the photo- 
graphing of stellar spectra, a frequency-control device and a frequency regulator are 
planned and will be secured when war requirements for materials are relaxed. 

The two slow speeds in right ascension corresponding to the slow speeds in declination 
are applied to the driving worm by means of an additional motor. The telescope is also 
provided with two quick motions in right ascension, one 39. 0° and the other 8.5° per min- 
ute of time. These motions are achieved by a two-speed motor, driving a spiral pinion 
which engages a spiral gear mounted on the polar axis. A solenoid-operated clutch re- 
lease, interposed between the spiral pinion and the motor, leaves the pinion free to turn 
when the slow motions are in operation. 

For both right ascension and declination motions, the clutches and the slow-motion 
clamps are electrically interlocked so that both quick and slow motions cannot be en- 
gaged at the same time. However, both may be released at the same time to permit man- 
ual moving of the telescope. This release, which is controlled by a switch near the eye-end 
of the guide telescope, affords convenient means of rapidly setting the telescope and of 
adjusting the balance. 

The electrical control—For convenience in setting, the motions referred to above are 
controlled from portable push-button stations. Push buttons are arranged in two groups, 
one group for quick motion and circle setting and the ether for centering the guide star 
and guiding. The switches for turning the dome are placed in the box with the first group. 


THE DOME 


The 28-foot dome is constructed of structural steel covered by insulating material and 
made waterproof by an outside covering of sheet copper. It is provided with a double 
shutter opening, 6.5 feet wide, making possible 3-hour exposures at the equator without 
movement of the dome. Two wind curtains are provided to cover that part of the slit- 
opening not in use by the telescope. The dome, shutter, and curtains are operated by 
motors, 


OPTICAL PARTS 


The 36-inch Pyrex-glass mirror blank.—This blank was manufactured by the Corning 
Glass Company and is remarkably free from bubbles or other defects. The mirror was 
ground and polished at the Warner and Swasey Company under the direction of Mr. C. 
A. Robert Lundin. Its optical surface was then coated with chromium-aluminum by Dr. 
Robley Williams. After three years the surface shows no sign of loss in its reflecting pow- 
er. 

Correcting lens.—This lens, also made under Mr. Lundin’s direction, has a diameter of 
24.5 inches. The distance of the depressed zone from the center is 8.48 inches, which ful- 
fils the condition for the minimum amount of glass to be removed to make the aspherical 
surface of the correcting lens. 

O ptical testing.—The spherical mirror was tested in the usual manner at the center of 
curvature. For the testing of the lens the spherical mirror was mounted in a vertical posi- 
tion, with the lens placed at its center of curvature, and a 60-inch-diameter optical flat 
mirror, placed a few feet beyond the correcting lens for a collimator. The curved surface 
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of the lens was turned toward the mirror, and a scale 24 inches in length, with the inch di- 
visions well marked, was placed horizontally just in front of the curved side of the lens. 
A carriage supporting the pin-point light-source and the knife edge was placed at the 
focus of the lens-mirror system. A diagonal flat, appropriately placed, brought the focus 
away from the optical axis. For each setting of the carriage the boundaries of the shadow- 
gram were read by means of the scale in front of the lens, thus securing quantitative 
measures of the diameter of the circle of confusion and means for ascertaining the zone 
of the correcting lens, which required further figuring.. This method of testing was sug- 
gested by Dr. James G. Baker. The accepted lens yielded a circle of least confusion of 25u 
in diameter with a negligible amount of astigmatism. Numerical details of the optical 
system are given in Table 1. 
TABLE 1 


SUMMARY OF NUMERICAL Data 
Mirror (Pyrex): 


Diameter (clear aperture)................., 36 inches 
Correcting lens (Vita glass): 
Distance of depressed zone from center of lens........... 8.48 inches 
Distance from mirror to plane of lens (plane surface outside) 168.04 inches 
Photographic plate: 
Ratlius'‘of ‘curvature of plate. .... 84.00 inches 
Sagitta of bent plate at center....................0.... 2.34 mm 
Over-all diameter of 8.56 inches 
Objective prism: 
Length of spectrum from H6 to 332 mm 


THE 24-INCH OBJECTIVE PRISM 


The Bausch and Lomb Company furnished the disk for the 24-inch objective prism. 
It was manufactured from one piece of dense flint (DF-2) glass weighing 260 pounds, 
molded to shape in a deep furnace utilizing a pot design to produce the wedge shape. The 
disk was free from deep striations and air bubbles. The diameter of the finished disk, be- 
fore being reduced in the optical shop, was 26.75 inches; the thickness varied from 3.0 
to 4.3 inches. 

The diameter of the finished prism is 24.5 inches, the graduated thickness from 0.75 
to 2.5 inches, producing an angle of 4°. The two surfaces were figured flat to within an ac- 
curacy of one-tenth of a wave length of sodium light. The annealing of the glass seems 
to have been excellent. This became apparent from the quality of the spectra secured on 
first trial, no additional local figuring of the prism being necessary. The spectra (PI. 
XXXII) compare favorably with slit spectrograms of low dispersion. Spectral classifica- 
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tion can be made with considerable accuracy, and absolute-magnitude characteristics are 
readily visible. 

The prism-cell mounting is so constructed that when in place it may be easily rotated 
through any desired angle in a plane perpendicular to the optical axis of the telescope. 
The cell with the prism forms a symmetrically balanced mass of 150 pounds. A 26-inch, 
150-pound ringweight, situated in front of the correcting lens, is first removed from the 
telescope when the prism is to be mounted; thus any rebalancing of the instrument is 
avoided. Optical data for the prism are given in Table 1. 

Exposure and limits.—With Eastman 103a-O emulsion and 10 minutes’ exposure, the 
limiting magnitude is 18.5. Using Eastman 103G emulsion with a No. 15 Wratten filter 
and with an exposure of 20 minutes, the limiting photovisual magnitude obtained is 15.5. 
A No. 22 Wratten filter and Eastman 103-E emulsion with 20 minutes’ exposure gave for 
a limiting red magnitude 16.8 (estimate). Usually plates are developed in a tank with 
D19 at 65° F. for 5 minutes. 

Spectra of stars 0.2 mm wide, of limiting magnitude 12.5, have been obtained with 
Eastman 33 emulsion and exposure of 30 minutes. 


In planning this instrument I am indebted to many friends for inspiration and advice, 
and it is fitting to express here my sincere thanks to all. Particularly, I wish to thank Dr. 
Theodore Dunham, Jr., Dr. Walter Baade, and Dr. J. A. Anderson for their generous 
help during the preliminary stages of planning; Dr. James G. Baker for his generous as- 
sistance in the optical tests; Dr. Robley Williams for aluminizing the mirror; and my 
colleague, Dr. S. W. McCuskey, for his all-round help throughout the duration of the 

roject. 
, I am deeply indebted to Mr. Charles Stilwell, Mr. Warner Seely, and Mr. William 
Burger, of the Warner and Swasey Company, for making available for the construction 
of this instrument all the talent and long experience of the company in the field of tele- 
scope-making. Their only concern was to produce an instrument as nearly perfect 
in its optics, as complete in its design, and as convenient in its operation as was humanly 
possible. To attain this goal they intrusted Mr. E. N. Jennison with the design of the 
mounting, Mr. C. A. Robert Lundin with the optics, and Mr. H. L. Cook with the as- 
sembling and erection of the telescope, all highly skilful workers and outstanding leaders 
in their respective fields. , 
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Received March 7, 1945 


ABSTRACT 


The section of the Milky Way in Monoceros, between galactic longitudes 176° and 186°, represents one 
of the most regular parts of the Milky Way. The present investigation is based on spectral types and col- 
ors to m = 11.5 and on general star counts to m = 14.5. The general absorption in the least-obscured 
regions is found to be of the order of Apg = 0.6 mag/kpc, for distances up to 2000 parsecs. The final 
results and conclusions are presented in Section VII at the end of the paper. The star densities exhibit 
nowhere marked increases with increasing distance from the sun. For the stars at large, they are close to 
unity to a limiting distance of 2500 parsecs. The most notable negative density gradients are found for the 
F stars. The spectral density-curves indicate that there exist marked variations in the shape of the general 
luminosity function over the first 500 parsecs. 


The present paper represents the first of a series on spectral types, colors, and star 
counts for regions at low galactic latitude not affected appreciably by irregular local 
obscuration. The section between galactic longitudes 176° and 186° is notable for its 
uniformity of stellar distribution. If we avoid here the regions south of the galactic 
circle (where the extensions of the dark nebula in Orion can still be detected) we find 
nowhere marked regions of local obscuration. 

In the next few years students of galactic structure may well wish to pay special atten- 
tion to the more uniform sections of the Milky Way. The researches of the last fifteen 
years have made it possible to delineate the more important dark nebulae; and it has 
further been shown that, for the smooth and relatively unobscured sections of our Milky 
Way, one may derive from observations of color excesses reliable information concerning 
the average coefficient of absorption to a distance of 2000 parsecs. In a first attempt at a 
study of the regular features in the stellar distribution for the nearer portions of our 
galaxy, we should now concentrate on those regions for which the absorption is probably 
small and not subject to sizable fluctuations. 

The results to be presented here and in some of the papers to follow are based on pub- 
lished and unpublished spectra classified by Miss Cannon, on color indices derived from 
blue and red photometric comparisons, and on star counts from blue-sensitive plates. 
Since our work depends to such a large extent on spectral types and color indices, we 
shall first examine the principal sources of error which affect the analysis of data of this 
sort. 

1. Systematic errors in spectral classification.—During recent years it has become more 
and more evident that most systems of spectral classification are subject to systematic 
errors, frequently dependent on the apparent magnitude. Miss Cannon’s system of 
classification is probably affected as little as any system by errors of this sort, but even 
here there has been some indication of the presence of systematic errors depending on 
apparent magnitude. We have tried to eliminate the effects of such errors by the com- 
bination of data for several neighboring spectral subdivisions. 

2. The lack of distinction between giants and dwarfs.—This represents one of the most 
serious problems in the use of the HD or HDE spectra of the later spectral types. We 
have tried to predict the approximate numbers of dwarfs from data on the bright stars 
and thus to obtain at least approximate information on the distribution of the giants. 

3. Uncertainties in the magnitude and color systems.—In spite of continuous checking 
and cross-checking we have found no indications, in the blue or red standard magni- 


* A paper presented at the symposium during the dedication of the Observatorio Astrofisico Nacional 
de Tonanzintla in Mexico, in February, 1942. 
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tudes, of systematic errors in excess of 0.2 mag. It is believed that the final color indices 
are in no case affected by systematic errors as high as 0.15 mag. 

4. Errors in the computed absorptions.—The computed absorptions are subject to error 
for two additional reasons. The first is that there exists for several spectral types some 
uncertainty about the valye of the average color index for unreddened stars. The errors 
from this source could hardly be much over 0.1 mag. in the blue-red color system. The 
second and more serious source of error arises from lack of definite information about 
the value of the conversion factor from color excesses to total absorptions. We have used 
the factor 3.0 for our blue-red system, which agrees substantially with the 1/)-law of 
scattering, but we realize that future work may make it necessary to use a different 
value. A comparison between color excesses and galaxy absorptions in intermediate 
latitudes, however, makes it seem unlikely that the factor is in error by as much as 
30 per cent; this is further confirmed by color studies for stars viewed through isolated 
dark nebulae. Since we have limited ourselves largely to regions of slight obscuration, 
our computed absorptions should be in error by not more than 0.5 mag. 

5. Errors in absolute magnitudes and dispersions.—It is extremely difficult to assign 
average absolute magnitudes and dispersions for each HD spectral class. The establish- 
ment of an accurate two-way system of spectral classification, in which the absolute 
magnitude and spectral class are read off simultaneously, represents one of the most 
urgent needs. Pending the completion and large-scale operation of such a system, it 
would, however, be a mistake not to make the best possible use of the available spectral 
data. For the present investigation the uncertainty in the average assumed absolute 
magnitudes is large, but we should remember that an analysis based on HD and HDE 
spectra is already much more accurate than one based on general star counts alone. 

If we consider the total effect of these various sources of error, we find that the per- 
centage error of the computed star densities for distances within 1000 parsecs of the sun 
is generally of the order of 30 per cent and that in extreme cases it may be as high as 
50 per cent. It does not seem that at this moment any investigation of galactic structure 
can claim a higher accuracy, but even with this uncertainty we can decide about the 
general characteristics of the galactic density distributions in what is generally referred 
to as the “local system.” From the accumulating observations we should soon be able to 
decide whether our sun is in a region of excessive, average, or subnormal star density 
with respect to its surroundings. 

In a paper presented at the symposium for the Jiediention of the McDonald Observa- 
tory,! the stellar distribution for the Milky Way in Monoceros was investigated on the 
basis of the spectral types, colors, and star counts that were available at that time. We 
present here the complete data and analysis for the entire section. . 


I. SPECTRAL TYPES AND COLOR INDICES 


Miss Cannon made available to us in advance of publication the HD charts for the 
section between a = 6"22™ and 7°08™ and 6 = —5° to +3°. The 1939 paper! was largel 
based on the color indices of 1054 stars in two regions, one centered at a = 645™, 6 = 0° 
the other at a = 6"57™, 6 = +1°, which together cover an area of 27.8 square degrees. 
The color indices for these two regions were determined from blue and red polar com- 
parison plates, taken with the 8-inch Ross (IR) camera at the Oak Ridge station of the 
Harvard Observatory. The blue series were made on Cramer Hi-Speed plates without a 
filter; the red, on Eastman 1-C special plates with a ciné-red filter. Dr. Cuffey and Miss 
Cherry were responsible for the taking of the plates and for their measurement and re- 
duction. As standards of magnitude, the polar sequence on the International Photo- 
graphic Scale? and the red polar sequence determined by Mrs. Payne-Gaposchkin* were 
used. 

1A p. J., 90, 249, 1939; Harvard Reprints, No. 180. 

2 Trans. I.A.U., 1, 71, 1922. 8 Harvard Ann., Vol. 89, No. 5, 1935. 
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Dr. MacRae, Mr. Fleischer, and Dr. Mowbray took the photometric plates for 
four more centers. The plates for the six centers, together, provide the material for the 
determination of color indices and magnitudes for the area of 96 square degrees covered 
by the Henry Draper Extension in this section of the Milky Way. The extremes of galac- 
tic longitude and latitude‘ are 176° and 186°, and —4° and +7°. For the four centers of 
the second group, the blue series were again taken on Cramer Hi-Speed plates; but for 
the red series, Eastman 103 E plates and a ciné-red filter were found to give better re- 
sults. In addition to the pole, Harvard standard regions C3, C4, and C5 were frequently 
used as standards. As photographic standards, we used those of Miss Leavitt,® corrected 
according to HB 781; and as photo-red standards, those of Mrs. Payne-Gaposchkin,® 
reduced to the color system of the IR, Eastman 103 E plates, and a ciné-red filter. The 
photo-red magnitude found from Eastman 1-C special plates are on the same color 
system as those of H.A., Vol. 89, Nos. 5 and 8, but the magnitudes from the 103 E 
plates refer to a somewhat longer effective wave length. The color equation of the 103 E 
system of magnitudes was determined by first plotting the curve relating the scale meas- 
urements and the standard magnitudes for the polar stars with color indices <0.40 
mag. The magnitudes of the polar stars with C.I. > 1.20 mag. were then read from this 
curve, and the differences were computed between these magnitudes and those listed in 
H.A., Vol. 89, No. 5. The results showed that the magnitudes of H.A., Vol. 89, No. 5, 
needed a correction of 

Am = —0.08XCLI., 


in order to reduce them to the system of the 103 E plates. This correction was applied 
to the magnitudes of the standard stars in the drawing of the reduction curves for the 
standard 103 E plates. All measurements and reductions were carried out by Mrs. Arons. 
For each center, five blue and five red comparison plates were taken. The total num- 
ber of color indices available is 2800. 

All plates in a series were preferably taken at similar zenith distances; but where it 
was necessary to apply a correction, we used the formula 


Am = 0.40 (sec ¢—1) 
for the extinction correction to the blue magnitudes and 
Am = 0.20 (sec ¢ — 1) 


for the red magnitudes. All measurements were confined to the central areas of the plates, 
for which no correction for distance from the center is necessary. The remaining error 
should not exceed 0.10 mag. 

As a check on the consistency of the blue magnitudes for the six IR centers, these 
magnitudes were intercompared on one of Dr. Baker’s star-count plates (see next section) 
with an exposure time of 30 minutes. This investigation showed that the blue magnitudes 
in five of the six centers were consistent to within 0.10 mag. but that the blue magni- 
tudes of one region seemed to differ systematically from those for the other centers. 

It happened that there was considerable overlap between the IR plates for the region 
with the questionable magnitudes and a neighboring region for which there was no 
indication of errors in the photographic magnitudes. The existence of a systematic differ- 
ence was confirmed by measurements made on the overlapping portions of the IR plates. 
A comparison between the IR sequence and the sequences used for the star counts also 
tends to confirm the presence of a small systematic error in the photographic magnitudes 
of this particular one of the six centers. 


4 All galactic latitudes and longitudes refer to the Harvard pole; a = 12540™, § = +28°. 
5 Harvard Ann., Vol, 71, No. 4, 1917. 6 Harvard Ann., Vol. 89, No. 8, 1938. 
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The only doubt about the justification for the application of a systematic correction 
(which at maximum amounts to +0.20 mag. for m = 10.5) arises from the good system- 
atic agreement between the average residuals for the five IR plates originally used to 
derive the standard sequence. We have used the following corrections to the blue magni- 
tudes for the doubtful center: 


Mog Am Mog Am 


These corrections have been applied to the color indices listed in Table 2 for regions 
9, 10, 11, and 12. By applying the full correction we may make the color indices for 
regions 9-12 slightly too red, since the corrections are admittedly the maximum values. 

There are no indications for the presence of systematic errors in the photo-red magni- 
tudes; but, because of the lack of red check plates, this result is not too certain; there 
remains the possibility of systematic errors of the order of 0.10 mag., in the photo-red 
magnitudes. 


II. STAR COUNTS 


The plates for the star counts were very generously provided by Dr. R. H. Baker, of 
the University of Illinois. Ten plates were needed to cover the region between a = 5545™ 
and 715", 6 = —9° and +9°. The centers of the plates are, on the average, 4° apart; 
and, since the countable area free from distance correction is close to 60 square degrees, 
the overlap between the counts from one plate to the next is considerable. 

As standards, we used the magnitudes in S.A. 98, which falls well within the region. 
The magnitudes were taken from H.A., Vol. 101; and the corrections listed by Seares’ 
were used for the reduction to the International Photographic Scale. Because of the con- 
siderable overlap from plate to plate, we encountered no difficulties in the transfer of 
sequences. Even for the plates centered farthest away from S.A. 98, it is unlikely that 
the standard magnitudes could be in error by as much as 0.2 mag., and the errors are 
certainly smaller for the nearer plates. Twenty stars were considered a norm for each 
sequence. The smallness of the systematic errors in the magnitudes was confirmed by a 
comparison with the standard blue magnitudes for the color program that were obtained 
by direct polar and C-region comparisons (see preceding section). All photometric cali- 
brations and counts were carried out by Mrs. Arons. 

Two checks were made on the accuracy and internal consistency of the counts. First, 
the counted numbers of stars to different limits of apparent magnitude were compared 
for all overlapping parts of the plates. It was found that the blue magnitudes used as 
standards in the star counts and in the work on the colors are generally consistent to 
within +0.10 mag. but that a few deviations as large as 0.15 mag. exist. Second, at the 
conclusion of the program Mrs. Arons recounted for each of the subdivisions A-R (see 
following section) an area of 4 square degrees. Small corrections were applied to reduce 
the star counts to one single uniform system. 

The absence of an appreciable constant error of counting was proved by a compari- 
son of the star counts for S.A. 98 published by van Rhijn® with those for the same center 
made by Mrs. Arons. For the brighter magnitudes a further check was obtained by 
comparing the star counts with the results of counts in the same area made in the cata- 
logue of photographic magnitudes for the HDE stars. It may be of interest to mention 
at this point that the Henry Draper Extension is essentially complete to mpg = 11.5 in 
this part of the sky. 

The section covered by the star counts was divided into eighteen areas, over each of 
which the stellar distribution appeared to be reasonably uniform. This division was 


TAp. J., 61, 303, 1925; Mt. W. Contr., No. 289. 8 Groningen Pub., No. 43, 1929. 
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made on the basis of the counts to various magnitude limits, from an inspection of long 
exposure plates reaching to faint limits, and with constant reference to Ross’s Milky 
Way Atlas. Figure 1 shows the division in areas A- R. 

In Table 1 we have listed for each area the final values of log V(m), N(m) being the 
average number of stars per square degree brighter than apparent photographic magni- 
tude m for each particular area. 


Ill. THE ANALYSIS FOR SPACE REDDENING 


The spectrum-color survey covers a smaller area than the star-count survey. We felt 
that it was desirable at first to subdivide the region with known colors and spectral 
types—which is really the heart of this section of the Milky Way—somewhat more finely 
than shown in Figure 1. The region with known colors was therefore divided into twenty 
subregions. 


TABLE 1 
VALUES OF LOG N(m) FOR PHOTOGRAPHIC MAGNITUDES IN THE MONOCEROS REGION 
| 

Area l b 8.0 | 8.5} 9.0 | 9.5 | 10.0 | 10.5) 11.0 | 11.5) 12.0 | 12.5) 13.0 | 13.5) 14.0 | 14.5 
180°} + 7°| 0.25/0.39| 0.58/0.84| 1.12]1.31) 1.49/1.64 1.79;1.98) 2.11/2.42) 2.49)..... 
180 | + 5 48) .69} 1.53/1.71) 1.85)/2.00 2.16|2.45 2.68) 2.93 
180 | + 2 38) .73} .69/0.87) 1.19)1.43) 1.61/1.81) 1.95)2.12) 2.29)2.48) 2.83)..... 
185 | + 3 38) .74| .83/0.98 1.19]1.39 1.61/1.79| 1.94/2.11) 2.28/2.51| 2.78)2.99 
475°) 73/0.95} 1.13}1.32} 1.51]1.72) 1.88'2.03} 2.19|2.40 2.71 
180 | — 3 | 0.25) .59/0.87| 1.11]1.32| 1.50)1.69) 1.84/1.98) 2.15/2.35| 2.61)..... 
H .48} .73/0.96] 1.20)1.35} 1.59)1.70) 1.88/2.02) 2.19)... .).....J..... 
K 0.42} 0.77}1.01) 1.17|1.37) 1.54/1.69) 1.87/2.07| 2.34)..... 
| 183 0.49} .46/0.86} 1.12/1.35) 1.53)1.69| 1.84/1.94) 2.14/2. 34) 2.60/2.89 
M 188 .91/1.08} 1.28]1.48} 1.64/1.82) 1.92/2.07| 2.22/2.41) 2.67)..... 
\ | 1.00} 1.21]1.37| 1.49)1.62) 1.76|1.90) 2.05|2.21 2.58) 
.30/0.62} 0.83}1.10) 1.23)1.40) 1.55)1.72) 1.84/2.04) 2.39)..... 
0.02/0.56} 0.81/1.09) 1.30)1.46) 1.51/1.73) 1.83/2.03) 2.29)..... 


The data on the color indices and excesses in color areas 1-8 have already been repro- 
duced in the paper presented at the time of the symposium for the dedication of the 
McDonald Observatory.' The resulting color excesses were summarized in Table 2 of 
that paper. The new material gives similar data for twelve additional regions. The result- 
ing color indices are shown in Table 2 for the stars with spectral types B8-A7 and ap- 
parent photographic magnitudes between 10.0 and 11.5. 

When the averages for the separate regions were combined, it was found that the 
results for regions 9-12, 13 and 14, 15-17, and 18-20 were strictly comparable; and in the 
final analysis we have therefore combined the colors for these areas into single averages. 
In agreement with the earlier investigation we have assumed the following values for 
the normal colors of various spectral types: 


The final color excesses are listed in Table 3. 
In spite of the somewhat larger base line of the new color system we are still assuming 
that a good approximation to the total absorption, K, is obtained from the relation 


K=3XC.E. 
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The distances in the last column have been computed from the relation 
Slog r= (m—-M)+5-K. 


We have assumed throughout that M equals 0.0 for the B8—A0 stars and that M equals 
+1.2 for the A1—A7 stars. 

The total number of early B-type stars for the region was small, but the observed 
colors of these stars have been taken into consideration in deciding upon the most likely 


TABLE 2 
MEAN COLOR INDICES IN THE MONOCEROS REGION 
es 9 10 11 12 
186° 186° 185° 185° 
+5° +3° +2° +1° 
Sp. B8~A0 A1-A7 B8-A0 A1-A7 B8-A0 A1-A7 Bs-A0 A1-A7 
m 
+0.16(10)| +0.20 (7)} —0.01(18)} +0.15(19)| +0.01(21)} +0.04(16)} +0.04(21)} +0.14(14) 
+ .10(18)} + .22(16)} + .08(28)} + .22(22)} + .07(30)} + .10(26)| + .03(27)| + .11(22) 
+0.16(20)} +0.30(24)} +0.18(27)| +0.22(34)} +0.06(31)} +0.17(30)| +0.10(22)| +0.19(32) 
a 183° 184° 184° 182° 
+1° —1° —? 
Sp. B8-A0 A1-A7 B8-A0 A1-A7 B8-A0 A1-A7 B8-A0 A1-A7 
m 
+0.01(7) | +0.09 (9)} +0.02(11)} +0.22 (5)} +0.02 +0.19 (7)} +0.03(11)} +0. 18(10) 
00(5) | + .11 + .06(16)] + .20(18)} + .09(17)} + .22(14)] + .10(15)| + .24(14) 
+0.09(5) | +0.23(12)} +0.05(16)| +0.28(17)} +0.31 (9)| +0.20(11)} +0.21(22)} +0.22(26) 
17 18 19 | 20 
—4° —4° 0° 
B8-A0 | A1-A7 | B8-A0 | A1-A7 B8-A0 | A1-A7 B8-A0 Al-A7 
™m 
+0:05 (9)|/ 40.24(11)| —0.27(1) —0.07(18)} +0.09(14)} —0.19 (2)] +0.10 (6) 
+ .12(11)} + .22(14)} — .O5(1) |........... + .04(14)} + .04(13)} — .04(11)| + .09(13) 
+0.28 (8)| +0.33(23)} +0.16(2) | +0.16(4) | +0.17(12)} +0.17(12)} +0.01 (6)| +0.17(16) 


values for the interstellar absorption for each area. There are ten stars of spectral class 
BS and earlier for which Stebbins, Huffer, and Whitford’ have measured color indices 
with the aid of the photoelectric cell and fifty-nine such stars for which our plates yield 


color indices. 


The stars with photoelectric color indices fall almost entirely in regions 17 and 19. 
The five stars in region 17 indicate no color excess at r = 400 parsecs and one amounting 
to +0.50 on our scale at 2000 parsecs, which is in general agreement with our results in 
Table 3. The three stars in region 19 are considerably reddened, the average color ex- 


Ap. J., 91, 20, 1940; Mt. W. Contr., No. 621. 
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cesses on our scale being of the order of +0.50 at r = 1000 parsecs. This is larger than 
the excesses according to Table 3 but, as we shall see presently, is in agreement with 
our colors for the early B stars. 


TABLE 3 
COLOR EXCESSES IN MONOCEROS 
B8-A0 A1-A7 
Mog 
Cz. | K | r C.Ei K r 

Regions 9-12 (/=185°, b=+2°) 
+0.20 0.6 850 +0.08 0.2 580 
.22 0.7 1020 690 
+0.37 1 1070 +0.17 0.5 790 

Regions 13 and 14 (J=184°, b=0°) 
+0.16 0.5 890 +0.07 0.2 580 
SR te +0.22 0.7 1290 +0.11 0.3 870 

Regions 15-17 (J=183°, b=—2°) 
+0.18 0.5 890 +0.15 0.5 500 
+0.42 980 +0.20 0.6 760 

Regions 18-20 (/=178°, b=—2°) 
+0.02 0.1 1070 +0.04 0.1 600 
4 1180 .03 Be 760 
+0.25 0.8 1230 +0.12 0.4 830 


The early B’s are especially numerous in regions 9-12 and 18-20. In regions 9-12 the 
average color indices for the stars with spectral types BO-BS are: 


Mog n 

+ .02 8 


The average unreddened color index for the BO-B5 stars is probably of the order of 
—0.35 mag., and their mean absolute magnitude close to —2.5. The early-type stars 
suggest, then, for regions 9-12, a total photographic absorption of the order of 0.6 
mag/kpc, in agreement with the value derived from Table 3. 
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In regions 18-20 the early B stars are quite reddened. For ten stars with apparent 
magnitudes 9.5—10.5 we find: 


53 2 


These data suggest a total photographic absorption of the order of from 1.0 to 1.5 mag/ 
kpc, which is definitely larger than the value in Table 3 but in agreement with the photo- 
electric colors. 

The results on the reddening for the Monoceros—Canis Major region are confirmed, 
in a general way, by the work of Cuffey'® on colors in three galactic clusters in Monoc- 
eros and Canis Major. One of Cuffey’s clusters, NGC 2324, which shows a color excess 
of 0.30 mag. at a probable distance of 3320 parsecs, falls in region 7 of our color survey. 


TABLE 4 
TOTAL PHOTOGRAPHIC ABSORPTIONS IN MONOCEROS 
Area i i | 3 ‘hes 
180° 0° | 1.1. mag/kpcl| 13-14...... 134° | 0° | 0.5 mag/kpc 
Sergey 181 +3 0.9 } 15-17...... 183 —2 1.0 
185 | | | 
| | 


The second cluster which interests us here, Messier 50, lies directly on the galactic 
equator 3° south of the southern border of our field. The observed excess of 0.30 mag. 
at 1210 parsecs agrees well with the averages for regions 13 and 14 in Table 3. 

By a combination of the results from various sources we determine the total photo- 
graphic absorptions as listed in Table 4. The values of Apg refer generally to limiting 
distances of from 1500 to 2000 parsecs from the sun, depending on the value of the total 
obscuration. The only regions for which there exists some conflict between the color 
excesses derived from the A and B stars are regions 18-20. We have here assumed the 
mean of the value of A,, suggested by the early B’s, Apg = 1.2 mag/kpc, and that for 
the B8—AO stars, Ap, = 0.5 mag/kpc. 

In our discussion of star counts and spectral distributions we have used a somewhat 
coarser system of subdivision than for the color work. Subdivision C of Figure 1 covers 
regions 1, 2, and 4 and half of 3 and 6; B covers 5, 7, and 8 and the remaining parts of 
3 and 6; D takes in 9, 10, 11, 12, and 13; F includes 14, 15, 16, 20, and half of 19; and 
G, 18 and the other half of 19. The small subdivisions O, P, and H all fall in region 17. 

The color data lead to the variations of K (the average total photographic absorption) 
with distance from the sun shown in Table 5. The colors of the distant stars in region C 
are apparently generally redder than those of D. Since we have already applied the 
maximum positive correction to the blue IR magnitudes in D, it is unlikely that the true 
reddening is in excess of that listed here. By using the value Ape = 0.6 mag/kpc (see 
next section) in region D we are therefore certain that we are applying the maximum 
allowable absorption correction for this region. The value Apg = 1.0 mag/kpc appears 
to represent an upper limit to the photographic absorption in region C. 


10 4 p. J., 94, 55, 1941. 
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The star-count data suggest that region F is affected by an obscuring cloud absorbing 
about 0.3 mag. The resulting extra reddening would be small and is apparently lost in 
the errors in the color indices. This same cloud covers region G, for which the total ab- 
sorption will be found to amount to 0.65 mag. from star-count data. This absorption is 
quite apparent from a comparison between the colors in G and C and D. The colors sug- 
gest that the obscuring cloud lies probably at between 500 and 700 parsecs from the sun. 

Regions B and H are further away from the galactic equator. For B at 6 = +5° the 
absorption amounts to Apg = 0.6 mag/kpc to a distance of 1250 parsecs. The small 
color excess of NGC 2324 suggests that there is little absorption beyond 1500 parsecs in 
the region, but this needs further confirmation. In region H at b = —5° the color excesses 
show an average absorption of A,, = 1.0 mag/kpc to 1500 parsecs. 


TABLE 5 


VARIATION OF TOTAL PHOTOGRAPHIC ABSORPTIONS, K ) WITH DISTANCE 
FROM THE SUN (FROM COLOR EXCESSES) 


REGION 
DISTANCE 
Cc D F G B H 
(180°, +2°) (185°, +3°) (180°, —3°) (175°, —3°) (180°, +5°) (183°, —5°) 
0.4 0.1 0.2 0.0 0.3 0.3 
0.5 0.3 0.4 0.3 0.5 
0.9 0.5 0.5 6 0.7 
| tA 0.6 0.6 0.9 


IV. DENSITY ANALYSIS FOR REGION D 


We found in the preceding section that, for the area covered by the color survey, the 
smallest reddening at the galactic equator apparently occurs in region D. Spectral types 
and color indices are available for an area of 20.7 square degrees in D; and, since these 
stars all fall between galactic latitude 0° and +5°, we find D an ideal region for a detailed 
investigation of the variations of star density with increasing distance from the sun. 

The photogaphic magnitudes that have actually been used are those corrected for 
the systematic error discussed in Section I. The absorption characteristics ofthe region 
are assumed to be described satisfactorily by the value Ap, = 0.6 mag/kpc to a limiting 
distance of 2000 parsecs; because of our correction to the original photographic magni- 
tudes, this coefficient represents an upper limit. The available information on the sur- 
face distribution of the stars in D is summarized in Table 6. 

The numbers of stars with photographic magnitudes fainter than 8.5 are generally 
sufficiently large for a satisfactory density analysis, but the total area covered by D is 
so small that the numbers of brighter stars are too much affected by random fluctua- 
tions. We have therefore strengthened the data for the brighter stars in two ways: 

1. For the stars with visual magnitudes brighter than 7.0 we have derived values of 
log A’(m) from a paper by Seydl."' The counts represent the averages over all galactic 
longitudes for the zone with galactic latitude 0° to +10°; A’(m) is again the number of 
stars with photographic magnitudes m — } to m + } reduced to an area of 100 square 
degrees. The smoothed values of log A’({m) are listed in Table 7. 


1 Pub. Nat. Obs. Prague, No. 6, 1929. 
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2. For the stars brighter than mpg = 8.5 special counts were made in the Henry Draper 
Catalogue for the region between a = 6"25™ and 7"05™; 6 = —7° and +3°. These counts 
are reproduced in Table 8. 

The values of log A’(m) from these three different sources were plotted against the 
apparent magnitude for each successive interval of spectral type, and smooth curves were 
drawn to represent the counts. The final values for log A’(m) are listed in Table 9. 


TABLE 6 


THE SURFACE DISTRIBUTION OF THE STARS IN REGION D 
(Total Area= 20.7 Square Degrees) 


SPECTRA 
Myg 
BO-B2 | B3-BS | B8-AO | Al-A7 | FO-F2 | FS-F8 | GO-GS | KO-K2 | KS-M 
<r 0 5 8 2 2 1 0 5 0 
8.0- 8.5... \o l 2 15 3 1 0 0 8 0 
J 24 13 9 1 1 3 
9 28 10 8 14 4 
9.5-10.0.... f f 41 29 6 3 11 13 5 
10.0-10.5.... \t \o0 64 47 15 10 19 20 8 
10.5-11.0.... f f 102 81 20 15 29 32 20 
137 96 15 19 28 25 17 
TABLE 7 
VALUES OF LOG A’(m) FOR THE ZONE 6=0° TO +10°, ACCORDING TO SEYDL 
Mog 
Sp. 
4.25 4.75 5.25 5.75 6.25 6.75 7.25 7.75 
B8-A0O..... 9.30 9.68 0.03 0.39 0.70 
Al-A7..... 9.40 9.54 9.78 9.97 0.22 
9.45 9.67 9.90 0.22 0.50 0.70 


The addition of the data in Tables 7 and 8 makes it possible to compute from Table 9 
the density distributions from the immediate vicinity of the sun out to the limit of our 
survey. 

Before we can compute the star densities for the stars of spectral types G and later, 
it is necessary that we attempt somehow to separate the giants from the dwarfs. Because 
of possible fluctuations in the galactic plane of the numbers of giants at large distances 
from the sun, it is impossible to depend on average percentages of dwarfs, such as are 
frequently used in studies for high and intermediate latitudes. Instead we shall try to 
predict, in various ways, the probable values of A’(m) for the dwarfs alone. Since scarcely 
any of the dwarfs with photographic apparent magnitudes brighter than 11.5 are beyond 
100 parsecs from the sun, we can assume, without too much risk, that the star density in 

. the galactic plane does not vary appreciably for these stars. The predicted numbers have 
been computed in three different ways: 
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1. From Seydl’s" tabulations we find the average numbers of GO-G5, KO—K2, and K5- 
M stars for m S 7.0 for the zone b = +30° to +70°. Since the percentages of dwarfs 
among these stars have been given by van Rhijn and Schwassmann,” we can compute 
the probable numbers of dwarfs for qur region. The predicted numbers are represented 
by a formula of the type : 


log A’(m)a=C+0.60 (mp, — 8.0) . 


TABLE 8 


VALUES OF A’(m) FOR A REGION IN MONOCEROS 
COVERING 100 SQUARE DEGREES 


Mpg 
Sp. 
6.75 7.25 7.75 8.25 
oS! 4 11 18 49 
1 0 + 12 
0 2 1 3 
1 1 0 6 
0 0 4 1 
+ eae 1 5 3 4 
1 0 0 1 
TABLE 9 
ACCEPTED VALUES OF LOG A’(m) IN REGION D 
SPECTRA 
Mog 
B8-A0 Al-A7 F0-F2 F5-F8 G0-GS K0-K2 K5S-M 
0.39 9.97 9.50 9.50 9.41 9.67 9.05 
0.70 0.22 9.80 9.72 9.66 9.90 9.28 
0.90 0.46 0.03 9.95 9.90 0.20 9.52 
0.71 0.30 0.20 0.18 0.50 9.77 
ye ee 1.45 0.93 0.55 0.43 0.45 0.76 0.07 
1.68 0.80 0.67 0.75 1.02 0.37 
Stites 1.93 1.42 1.05 0.88 1.07 ee 0.68 
BY eee 2.13 1.71 1.32 1.13 1.34 1.60 1.03 
2.30 2.03 4-55 1.65 1.90 1.33 
jo, 2.48 2.33 $277 1.62 1.93 2.05 1.70 
| LY ae 2.68 2.56 2.00 1.83 2.16 2.13 1.81 
2.82 2.70 2.00 2.03 2.26 2.20 2.05 


where C = 0.18 for the GO-GS5 stars and C = 9.37—10 for the KO-K2 stars. A’(m)a 
represents here the average number of dwarfs per 100 square degrees with apparent pho- 
tographic magnitude m — } to m + 4. The number of K5 and M dwarfs proves to be so 
small as to be negligible for all practical purposes. 

2. A study of the percentages of dwarfs in intermediate latitudes for the stars in the 
Yale zone catalogues'* shows that at b = +30° the average percentage for mpg = 


2 Zs. f. Ap., 10, 161, 1935. 18 Harvard Circ., No. 400, 1935. 
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9.0 — 10.0 is close to 55 per cent for the GO—GS stars and 10 per cent for KO-K2 stars. 
We find for the most probable values of C: C = 0.28 for the GO-G5 stars and C = 
9.40—10 for the KO—K2 stars. 

3. In the third method the numbers of dwarfs are predicted directly from the Jumi- 
nosity function for G and K stars published by van Rhijn and Schwassmann.” Again the 
basic assumption is made that the star density is equal to unity to the limit of our sur- 
vey for the dwarf stars. The resulting values of C are: C = 0.21 for the GO-GS stars and 
C = 9.54—10 for the KO-K2 stars. 

The agreement between the three different values of C is quite satisfactory. In the 
following calculations we shall use the formula 


log A’ (m) 0.25 +0.60 (my, — 8.0) 


for the GO—-G5 dwarfs and 
log A’(m) a= 9.45 + 0.60 (my, — 8.0) 


for the KO-K2 dwarfs. We neglect the dwarfs among the K5-M stars brighter than 
Mp, = 11.5 in the galactic plane, their number being too small to be significant. The 
extreme possible values for C are 0.10 and 0.35 for the G’s and 9.40 and 9.60 for the 
K’s. It is unlikely that the predicted numbers of dwarfs could be in error by as much 
as 20 per cent. 

By the use of the above formulae, together with the totals for giants and dwarfs com- 
bined from Table 9, we find the probable values of log A’(m), for the giants alone; these 
are listed in Table 10, together with the corresponding percentages of dwarfs. 

The variation of star density with distance from the sun has been computed with the 
aid of the method of Schalén and Lindblad.'* The mean absolute magnitudes and dis- 
persions around these means were assumed to have the values given in Table 11. The 
values listed there are the means per unit volume; the mean values M,, that are used 
in the Schalén-Lindblad method were computed with the aid of Malmquist’s formula.’ 
The Malmquist corrections were generally of the order of —1.0 mag. 

The computed star densities are given in Tables 12 and 13. In Table 12 we list the 
densities uncorrected for absorption effects; in Table 13, the corresponding values cor- 
rected for a photographic absorption of 0.6 mag/kpc. The densities are expressed as 
numbers of stars per 1000 parsec* as a unit. 

The densities in Table 13 represent the most likely values, though we admit that, 
because of the application of the maximum possible correction to the blue magnitudes, 
the value Ay, = 0.6 mag/kpc is decidedly an upper limit. We note that the densities 
for the B8-AO stars appear to increase slightly to a maximum of 0.18 at a distance of 
200 to 300 parsecs from the sun, that the star density between 500 and 700 parsecs is 
again very nearly equal to that near the sun, and that the star density drops rapidly 
between 700 and 1000 parsecs. There is little variation in the star density for the red 
giants, but a density decrease seems established for the A1—A7 stars. The most striking 
result is probably the steep negative gradients for the F stars, for which the density at 
250 parsecs is one-fourth that near the sun. 

We do well, however, at this point to remember the limitations of our andlysis as 
outlined at the beginning of this paper. The density analyses are very sensitive to un- 
certainties in the assumed values of the mean absolute magnitudes. In the absence of 
absorption all densities will be doubled if one chooses for the assumed mean absolute 
magnitude a value less bright by 0.5 mag. without making any change in the assumed 
dispersion. Since at present the assumed mean absolute magnitudes can be trusted only 


14 See Uppsala Medd., No. 37, 1928. 
16 Lund Medd., Ser. II, No. 22, 1920. 
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to the nearest half-magnitude, we see that the computed densities may possibly be in 
error by 50 per cent. 

We note further that we must be very careful in the application of the correct value 
of the Malmquist correction” at each apparent magnitude. Frequently the slope of the 
curve, relating the apparent magnitude, m, and the observed log A’(m)’s, decreases as 
we go to fainter magnitudes. If, in the analysis of the Schalén-Lindblad method, we 


TABLE 10 
VALUES OF A’(m)g AND PERCENTAGES OF DWARFS IN REGION D 
GO-GS5 Stars KO0-K2 Srars 
Mog 
Percentage Percentage 
A'(m), of Dwarfs of Dwaris 
0.88 42 3.06 3 
1.56 45 5.6 4 
11.9 46 38.2 4 
24.7 40 76 4 
40: | 80: 133 16 
TABLE 11 
MEAN ABSOLUTE MAGNITUDES 
AND DISPERSIONS 
Sp | Mo (pg) | a 

+1.0 | +10 

4.5 | 1.2 

+2.0 | +1.0 


apply then a single average Malmquist correction, we are in effect using too bright a 
value of M,, for the faint apparent magnitudes; and the densities for large distances 
will then come out smaller than the actual values. By following this procedure we might 
derive erroneous negative density gradients. 

In cases where it is difficult to apply Malmquist corrections because of irregularities 
in the (m, log A’(m) ) curves, it is generally preferable to do the analysis with the aid of an 
(m, log 3) table for each spectral interval. It is not difficult to construct the basic (m, 
log 7) tables so that the densities, found by trial and error, represent the average numbers 
of stars per unit volume of 1000 parsec* at each distance. In this method it is not neces- 
sary to compute the Malmquist corrections, since these are automatically accounted for 
in the construction of the basic tables. 
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With respect to the results of Table 13, we may doubt the reality of the excess density 
for the B8—A0 stars at 200 to 300 parsecs, but the drop in star density beyond 700 par- 
secs seems too large to be accounted for by uncertainties in the analysis or other sources 


TABLE 12 
SUMMARY OF STAR DENSITIES (NUMBER OF STARS PER 1000 PARSEC®) 
IN MONOCEROS REGION D 
(No Absorption Correction) 


| 
| SPECTRA 
d | | | 
B8-A0 Al-A7 FO-F2 | F5-F8 | gG0-gM 
deere 0.12 0.43 1.3 3.5 | 0.25 
12 35 0.33 |. 
0.47 0.8 22 
TABLE 13 


SUMMARY OF STAR DENSITIES (NUMBER OF STARS PER 1000 PARSEC®) 
IN MONOCEROS REGION D 
(Corrected for Apg=0.6 Mag/Kpc) 


| SPECTRA 
d 
| B8-A0O | A1-A7 | FO-F2 F5-F8 gG0-gM 
| 0.12 0.43 1.3 3.3 0.27 
| 15 34 0.93 1.9 27 
| 16 32 0.76 28 
18 0.55 0.8 27 
| | | 


of error. The negative density gradient for the A1-A7 stars is slight, but in all probability 
real. The rapid drop in star density for the F stars is too large to be accounted for by 
errors of any sort. Finally, for the red giants there is little indication of the presence of 
large variations in star density over a distance of 500 parsecs. It seems unlikely that 
the remaining uncertainties would be sufficiently large to change the character of the 
above general conclusions. 
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V. THE CHANGES IN THE LUMINOSITY FUNCTION 


For all analyses of general star counts we are forced to make use of an assumed general 
luminosity function. Not only is our knowledge of this function for the vicinity of the 
sun far from complete; but we have come to realize, more and more, that it may vary 
considerably as we proceed outward from the sun. The existence of such variations for 
the directions of the galactic polar caps is well established, and their approximate values 
have been derived. Until recently very little information was available on the variations 
in the galactic plane. From spectral data, such as those analyzed in the preceding sec- 
tion, however, it is possible to study these variations directly, and we shall attempt to do 
so for the direction of region D in Monoceros. 

The star densities of Table 13 are used as a basis for the calculations. With the aid 
of the mean absolute magnitudes and dispersions of Table 11 we can compute the 
contribution to the general luminosity function made by the stars of successive spectral 
classes at each distance. By adding the numbers for the spectral subdivisions we can 


TABLE 14 


CHANGES IN THE GENERAL LUMINOSITY FUNCTION 
FOR REGION D IN MONOCEROS 
(Values of log ¢’(M)) 


d (PARSECs) 
Mye | 

Van Rhijn | 50 100 | 200 300 500 

8.18 8.20 | 8.26 8.23 8.15 

na... 8.52 3.53 | 8.56 8.57 8.46 

8.88 8.86 8.85 8.83 8.72 


obtain the number of stars, ¢’(M), between absolute photographic magnitudes M —} 
and M + } per 1000 parsec* at various distances from the sun. We have reproduced 
these figures in Table 14. For comparison we have added the values of y’(M) according 
to van Rhijn’s standard function," corrected to refer to the same interval in absolute 
magnitude and unit volume. 

We note that there exists a satisfactory agreement between the standard function 
and that derived from our material at a distance of 100 parsecs from the sun. The only 
marked excess of our numbers over those represented by the standard function occurs 
at M = 0; apparently, we find here for near-by regions an excess in the number of A 
stars over the average included in the standard function. -- 

As we proceed outward, the shape of the general luminosity function changes ap- 
preciably. The values of log y’(M) at M = 0.0 remain steadily in excess of the standard 
value, but the log ¢’(M) values at M = +2.0 and +3.0 are well below the standard 
values at 200 parsecs from the sun. Table 14 also shows clearly how very limited is the 
information, on changes in the luminosity function that can be derived from spectral 
data complete to the eleventh or twelfth apparent magnitude. We shall need accu- 
rate spectral classifications to at least the fourteenth magnitude before we can be. 


16 Groningen Pub., No. 47, Table 6, 1936. 
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gin to investigate changes in the luminosity function to distances of 1000 parsecs with 
any degree of completeness. 


VI. THE ANALYSIS OF THE STAR COUNTS 


The data on the star counts have already been given in Table 1. For purposes of 
analysis it is desirable to have also available the values of A(m), the average numbers 
between m — 4 and m + 3 per square degree of the sky. The adopted values of A(m), 
computed from slightly smoothed values of N(m), are printed in Table 15. We shall now 
proceed to analyze these A(m)’s in detail. For a brief résumé of conclusions the reader 
is referred to the summary at the end of the paper. 

1. Regions C and D.—The two fields at 6 = +2° are close to the galactic equator; 
and the absorption coefficient there seems to be, on the average, equal to Apg = 0.7 
mag/kpc. With this value for the absorption assumed to hold to a distance of 2500 


TABLE 15 
ADOPTED VALUES OF A(m) IN MONOCEROS 


| | 

Area | l b | 9.0 | 10.0 11.0 12.0 13.0 14.0 
1 180° + 7° 4.28 i335 25:2 49.8 119 296 
ee 180 +5 4.87 13.4 29.6 53 157 500 
sc 180 + 2 5.56 18.4 36.1 69 199 601 
D... 185 + 3 5.26 15.0 35.8 68 202 624 
oe 175 + 3 5, 54 12.4 30.0 56 133 466 
Eee ee 180 — 3 4.89 135 28.0 49 123 457 
175 8.9 20.9 40.4 95 297 
H 183 —- § Bett 10.2 22.3 37.8 111 161 
177 7.9 21.4 32.4 90 360 
K 177 aa), 4.5 13.2 26.7 68 229 
183 0 5.24 14.7 30.6 42.9 119 361 
188 0 8.0 16.5 54.6 140 318 
176 12.9 18.2 37.8 91 392 
181 — 5 2.85 13.9 26.9 58 280 
183 16.5 24.3 55 169 
i 173 —10 1.34 6.3 12:3 16.6 45 151 


parsecs, we can now proceed to analyze the general star counts for these regions, and 
also the published star counts® for S.A. 98, which falls in region C. Because of the 
lack of information about the changes in luminosity for distances beyond 500 parsecs, 
we shall use the standard (m, log 3) table for b = 2° from the paper by Bok and MacRae 
presented at the symposium of the New York Academy of Sciences on “The Fundamen- 
tal Properties of the Galactic System.’’!” The density analysis with the aid of the (m, 
log w) table, corrected for absorption, yields the average densities listed in Table 16. 

The star densities in Table 16 show an excess within 500 parsecs of the sun; but 
beyond that, to the limit of 2500 parsecs for which our absorptions are more or less 
valid, we find a constant average star density very nearly equal to that near the sun. 
Because of the considerable variations in the shape of the general luminosity function 
in the galactic plane which were brought to light in the preceding section, it is not simple 
to interpret the densities of Table 16. The incipient decrease in the star density for the 
F stars (Table 13) apparently does not persist to great distances—or, at least, does not 
affect similarly the densities for the G and K dwarfs. A careful study of the distributions 
of the F, G, and K dwarfs between m = 11.0 and m = 15.0 for this part of the sky 
promises to yield exceedingly interesting results. 


17 Ann. N.Y. Acad. Sci., 42, 219, 1941. 
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2. Regions G, F, and E.—These three regions are at the same low latitudes as regions 
C and D. The deficiencies in the A’(m)’s, therefore, strongly suggest the presence of a 
dark nebula of nonuniform density, with a small absorbing power in F and E and a some- 
what larger one in G. The general star counts in F and E are very similar; and since no 
color excesses are available in E, we shall limit our analysis to-F and G, making the tacit 
assumption that the results for F apply in general to E. 

From the deficiencies in A(m) for G we find that the total photographic absorption of 
the extra dark. nebula amounts to 0.65 mag. The general star counts suggest that the 
dark nebula is relatively close to the sun, but the total area covered by the nebula is too 
small to yield a dependable value for its distance. We have found previously (see Table 5) 
that the color excesses show that the stars beyond 600 parsecs from the sun are abnormal- 
ly reddened by an amount corresponding very nearly to the total absorption found 
from the star counts. The numbers of stars are rather small in G, and random fluctua- ” 
tions in the counted numbers render it inadvisable to make an elaborate distance analysis 
from the numbers of stars of successive intervals of spectral class. The deficiencies are 
apparently mostly in the B8-A0 stars; the later A’s and the G-M stars are deficient 
beyond m = 10.0, and the F’s show little effect even at 11.0. The spectral-distribution 
data place the dark nebula at a distance of roughly 500 parsecs from the sun. 


TABLE 16 


AVERAGE DENSITIES CORRECTED FOR Apg=0.7 MAG/KPC 
FROM GENERAL STAR COUNTS 


Distance Region D S.A. 98 Distance Region D S.A. 98 


The total photographic absorption of the dark nebula affecting F is only 0.30 mag. 
from the star-count data. The effects of this small absorption do not appear clearly in the 
color excesses, but the similarity of the spectral distribution curves in G and F make it 
seem likely that the distance of the nebula is 500 parsecs or slightly less. It seems reason- 
able to suppose that the dark nebulae overlying G, F, and E are all part of the same 
system and that we may assign an average distance of 500 parsecs to this nebula. The 
nebula is probably a part of the system of dark nebulae associated with the bright nebula 
in Orion. 

3. Regions L and M.—Both regions are comparatively small in area. Region L, which 
lies along the galactic equator, is one of star deficiency. It should naturally be compared 
with D and F. The run of the star counts in L agrees very closely with that in F; and we 
conclude from the A(m)’s in Table 16 that, overlying region L, in addition to a general 
absorption A), = 0.7 mag/kpc, there is present a dark nebula absorbing between 0.3 
and 0.5 mag. at a distance of the order of 500 parsecs. 

The deficiencies in region M are (with the exception of that for m = 14.0) less than 
in region L. Region M is toward the edge of our counted region; but the counts, so far as 
they go, indicate that the distribution of stars and dark material to a distance of 2000 
parsecs does not differ markedly from that for D. 

4. Regions B and A.—These two fields at b = +6° and +8° fall in the relatively un- 
obscured regions north of the galactic circle. No galaxy counts are available at these 
latitudes, but at 6 = +10° an average photographic absorption of the order of 1.0 mag. 
is suggested. The colors in B show that the absorption in this field is probably repre- 
sented best by Ap, = 0.6 mag/kpc to a limiting distance of 1500 parsecs; Cuffey’s color 


. 
h 
its 
: 
W 
r 
) 
| || | 
| 
nd 
he 
cs, 
ae 
m 
ut 
ss 
in. = 
ion 
ple 
the 
t 
10 
ons 
k 
KY 
if 
~ 


298 BART J. BOK AND JEAN M. RENDALL-ARONS 


data for NGC 2324 render it unlikely that there is much absorption beyond 1500 parsecs 
for this direction. 

The (m, log x) tables for 6 = 6° and 6 = 8° by Bok and MacRae" were corrected 
for Apg = 0.6 mag/kpc to limiting distances of, respectively, 1500 and 1200 parsecs. 
The computed average densities are found to be very nearly equal in both fields. In 
judging these densities we should keep in mind that they are expressed in terms of the 
average densities at corresponding heights above the galactic plane in the direction of 
the galactic pole and that the usual drop in star density with increasing distance from 
the galactic plane has already been taken into account. The densities are: 


Distance Density Distance Density 


Again we find an excess for the regions within 400 parsecs of the sun. The deficiency 
beyond 500 parsecs seems well established by the counts, but we should remember that 
the densities are rather sensitive to the assumed value for the absorption coefficient for 
these directions. The deficiencies would disappear if the total absorption were increased 
by 0.5 mag. A careful study of the distributions of faint galaxies for the region between 
l = 180° and 190° and 6 = +5° and +15° would be well worth while. 

5. Regions J, H, and Q.—These three regions are south of the galactic plane at the 
same distance from the galactic circle as B and A. Color indices and excesses are avail- 
able for H, in which field the observed excesses suggest an average value Ap, = 1.0 
mag/kpc to 1500 parsecs. 

A comparison between the star counts in H and BSuggests that in H we have, in addi- 
tion to the same general absorption that overlies B, a dark nebula absorbing 0.7 mag. at 
a distance of not more than 500 parsecs from the sun. Because of the lack of a sufficient 
number of stars brighter than m = 10.0, it is impossible, from the general counts alone, 
to fix the distance of this nebula more carefully. The absence of appreciable color excesses 
for the A1-A7 stars with m = 9.5 to 10.5, together with the appreciable excesses for 
the B8—-AO stars with m = 9.5 and fainter, makes it very likely that the dark nebula is 
at 400-500 parsecs from the sun. 

The counts in J show a somewhat smaller deficiency with respect to B and A than did 
those in H. The average absorption of the dark nebula is of the order of 0.5 mag.; but, 
because of lack of spectral and color data, we can say little more than that it is within 
500 parsecs of the sun. Region Q falls at almost the same galactic latitude as does A; 
the deficiencies in A(m) are large, and the total absorption of the dark nebula here 
amounts to nearly 1.0 mag. 

6. Regions R and K.—Lack of data on spectra and colors makes it very difficult to 
say anything definite about these two fields. We are further hampered by the lack of a 
suitable comparison field north of the galactic circle at the same latitudes (9° and 11°). 
Also, the apparent absence of galaxies shows that we have here a rather heavily ob- 
scured region. Region R (in which division falls S.A. 97) is apparently most obscured and 
the counts here seem to demand a dark nebula with a total absorption of approximately 
1.0 mag. 

7. Remaining regions.—Region N is a small patch, somewhat obscured, in E; the 
excess absorption is less than 0.5 mag. Regions O and P are both very small. The counts 
and colors show that in all probability the absorption characteristics in these fields re- 
semble those in H very closely. Region J is at the edge of the counted area. Apparently 
the stellar distribution and the structure of the obscuring material in J do not differ 
materially from what we have already found for G. 
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VII. RESULTS AND CONCLUSIONS 


The investigation deals with the study of stellar distribution for a region of the Milky 
Way in Monoceros, where the surface distribution of the stars is quite uniform. The 
analysis is based on spectral types (HDE) and color indices for an area of 100 square 
degrees and star counts for 400 square degrees. It is found that the computed star densi- 
ties are subject to uncertainties of the order of 30 per cent of their values because of 
errors arising from difficulties in the classification of spectral types, photometric errors, 
errors in the derived absorptions, and lack of information about mean absolute magni- 
tudes and dispersions around these means. 

The basic data are given in Table 1 (see Fig. 1) and Table 2. The main results of the 
analysis are the following: 

1. Reddening is present in the apparently unobscured regions of the Milky Way 
between / = 178° and / = 188°, the average coefficient of absorption amounting to 
Ayg = 0.6 to 1.0 mag/kpc to a distance of 2000 parsecs. 

2. Star counts and observed color excesses suggest the presence of extended regions 
of obscuration south of the galactic equator. The total absorption is slight, varying be- 
tween 0.3 and 1.0 mag. The obscuration is probably caused by the low-latitude extensions 
of the Orion dark nebulae and is apparently located at 400-500 parsecs from the sun. 

3. The regions at b = +6° to +8° are comparatively clear, with an average Ay, = 0.6 
mag/kpc to limiting distances of 1500-1200 parsecs. ; 

4. The density analysis for a region of minimum obscuration near / = 185°, 6 = +2° 
(see Table 13), shows that— 

a) The densities of the B8-A0 stars are either the same or slightly in excess of those 
for the vicinity of the sun to r = 500 parsecs and then drop to one-half the value near 
the sun at r = 1000 parsecs; 

b) The densities of the A1—A7 stars drop to 60 per cent of that near the sun at r = 500 
parsecs; 

c) Steep negative gradients are exhibited by the FO-F2 and F5-F8 stars (density 
25 per cent) of that near the sun at r = 250 parsecs; 

d) The densities of the gGO-gM stars are almost constant over the first 500 parsecs, 
with a tendency to decrease. 

5. The star counts to m = 14.5 in the clear regions, combined with the available data 
in S.A. 98 to m = 18.0, show that for / = 180°, 6 = 0°, the average star density for the 
vicinity of the sun is in excess by 20-50 per cent with respect to the prediction by the 
standard luminosity function of van Rhijn and that the average density for r = 1000 
to 2500 parsecs is close to unity. 

6. The spectral density curves show further that the general luminosity function for 
the direction / = 185°, 6 = +2°, exhibits pronounced variations in shape over the first 
500 parsecs (see Table 14). 
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THE STELLAR DISTRIBUTION FOR TWO SOUTHERN FIELDS* 
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ABSTRACT 


The present paper is a continuation of the investigation for the Milky Way in Monoceros, published 
in the preceding paper. The fields are at / = 261°, b = +1° (S.A. 193), and at / = 249°, b = +1° (Field 
II). Selected Area 193 is apparently free from obscuration to a distance of 1600 parsecs, where as in Field 
II the average absorption to 1200 parsecs amounts to Apg = 0.5 mag/kpc. The final results and con- 
clusions are presented in Section VI at the end of the paper. The star densities in Field II are similar to 
those in the Monoceros field. The same holds for S.A. 193 with the exception that for this field the early 
B stars do not exhibit any marked decrease of space density with distance from the sun. Further evi- 
dence is presented for the existence of marked changes in the shape of the general luminosity function 
within 200-500 parsecs from the sun. 


During the last years of her life Miss Cannon spent much time and effort on the 
classification of the spectra of faint Milky Way stars in Carina and Centaurus. Since it 
is extremely important that we should obtain magnitudes and colors for most of these 
stars, the observers at the Boyden Station, under the direction of Dr. Paraskevopoulos, 
have taken, in recent years, many photometric series of plates for this part of the sky. 
A special star-count program is intended to supplement the work on the colors, magni- 
tudes, and spectra of the brighter stars.! Because of world conditions, it may be many 
years before we shall be able to complete the collection of the necessary plate material 
for the entire area, and we have therefore decided to measure the available plates for 
those areas for which enough material has already accumulated. The measurement and 
reductions were all done by Miss Wright, who collaborated with Bok in the analysis. The 
carrying-out of the research was made possible in part by a grant from the Draper Fund 
of the National Academy of Sciences. 

We shall report here the results for two centers, the first being S.A. 193 at a = 
11"24™, 6 = —60° (J = 261°, 6 = +1°), and the second, a field at a = 10°13™, 6 = 
— 56° (1 = 249°, b = +1°). The bright nebula near » Carinae lies about halfway be- 
tween these two centers, at / = 254°, b = 0°. Selected Area 193 falls in the rich but fair- 
ly uniform section of the Milky Way between the Carina nebula and the Southern Coal 
Sack (J = 270°, 6 = 0°). Field II falls in a uniform region directly west of the Carina 
nebula. In spite of the close proximity of the Carina nebula, the stellar distribution for 
this center differs rather markedly from that around the Carina nebula, or S.A. 193. 
The surface distribution of the stars in Field II is remarkably uniform, but the surface 
density of the stars is there very much less than for the section between / = 252° and 
1 = 270°. A comparative study of S.A. 193, the region near » Carinae, and Field II 
promises to be very interesting. 

The colors for S.A. 193 have been investigated by Dr. and Mrs. Gaposchkin,’? who 
found there a surprising preponderance of stars with negative color indices. The work 
by the Gaposchkins was done before Miss Cannon’s spectra had become available. Be- 
cause of the presence of systematic errors in the Potsdam spectral classes,* Mrs. Gaposch- 


* A paper presented at the symposium held during the dedication of the Observatorio Astrofisico 
Nacional de Tonanzintla in Mexico, in February, 1942. 
1 Trans. 1.A.U., 6, 282, 1938. 
oo p. J., 90, 321, 1939; Harvard Reprints, No. 181, 1939; also “Harvard Mimeo. Ser.,” Vol. 2, No. 3, 
3 Cf. Vyssotsky, A p. J., 93, 425, 1941. 
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kin did not attempt a complete density analysis; but her results for S.A. 193 show clearly 
that the light of distant stars in this region is not appreciably reddened. 


I. Cotor INDICES AND SPACE REDDENING 


The total number of available color indices in S.A. 193 is 732; that for Field II, 333. 
In each case we have covered an area of 6 square degrees. 

The photographic magnitudes for S.A. 193 were determined from two 30-minute 
plates with the 3-inch Ross camera (RB) and from two 10-minute exposures with the 
8-inch Bache camera (B). The photo-red magnitudes were determined from two 120- 
minute and two 15-minute exposures with the 24-inch Bruce camera (A). As standards, 
the Gaposchkin photographic and photo-red magnitudes‘ for S.A. 193 were used. The 
magnitudes were measured only on the central portions of the plates free from distance 
correction. The only correction that was applied to the magnitudes was one for color 
equation to reduce the photographic magnitudes of the late-type stars measured on the 
RB plates to the International Photographic system. Grating images were found useful 
for the measurement of the magnitudes of the brightest fifteen stars in the field. 

The system of the photographic and photo-red magnitudes used here is, of course, the 
same as that in the study of the Gaposchkins, since both are based on the same standard 
sequences. Further checks of this standard sequence will be possible when the Carina- 
Centaurus programs are completed, but the plates that are available lead us to believe 
that the system of colors should agree to within 0.20 mag. with the blue-red system for 
the C regions at +15°. 

The probable error of-an average final photographic magnitude is of the order of +0.06 
mag.; that of a final photo-red magnitude, of the order of +0.08 mag. The average 
probable error of the color indices in this investigation is therefore close to +0.10 mag. 

The color indices in Field II were obtained in very much the same way. The blue 
magnitudes were measured on one 60-minute and two 30-minute RB plates, calibrated 
with the aid of a transferred sequence from S.A. 193, and two 10-minute plates taken 
with the 10-inch Metcalf (MF) telescope in series with the standard regions E4 and F2. 
The magnitudes were all consistent, with the sole exception that a small correction was 
applied to the standards in E4. The red magnitudes were found from one 15-minute and 
two 20-minute A plates and two 15-minute MF plates; S.A. 192 served as a standard. 
Small corrections for extinction were applied, but there was no need to consider distance 
corrections in any part of the program. 

The results of the color analyses are summarized in Tables 1 and 2. The numbers in 
parentheses directly following each average are the numbers of contributing stars. The 
most striking effect in both tables is the preponderance of negative color indices among 
the early-type stars of all magnitudes. 

There always remains, especially for far southern declinations, the possibility of sys- 
tematic errors in the photometric system. We shall therefore inspect first the average 
color indices of those stars for which the excesses, if present at all, would be expected to 
prove rather small. We list for that purpose in Table 3, side by side, the normal colors 
for Harvard spectral types adopted by Mrs. Gaposchkin? in her study of red indices in 
southern Selected Areas and some averages for S.A. 193 and Field II. 

Table 3 shows, first of all, that the averages for S.A. 193 and Field II are remarkably 
similar. A comparison between the observed averages and the norma] colors indicates 
that, especially for the faintest stars, the observed colors allow a maximum positive cor- 
rection of, at most, +0.20 mag. This result, however, is uncertain because of possible 
small systematic errors in the spectral classes and, further, because of lack of informa- 
tion concerning percentages of dwarfs and average true color indices for unreddened 
regions. 


4 Harvard Ann., Vol. 89, No. 9, 1937. 
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The data on the color excesses of the early-type stars are summarized in Table 4. The 
spectral classes and the assumed true color indices are at the head of the table. The en- 
tries are the color excesses derived from Tables 1 and 2 and, in parentheses, the numbers 
of stars contributing to each average. 


TABLE 1 


OBSERVED COLOR INDICES IN S.A. 193 
(l=261°, b=+-1°) 


Mog 
Sp. 
8.00-8.49 8.50-8.99 9.00-9.49 | 9.50-9.99 10.00-10.49 10.50-10.99 11.00-11.49 
| —0.18(1) | —0.30( 1) | —0.14( 1) | —0.15( 1) 
—0.60( 4) | —0.55( 3) | —0.30( 6) | —0.20( 2) +0.19( 1) 
—0.53(1) | —0.36( 1) | —0.44( 4) | —0.45( —0.35( 9) | —0.32( 6) | —0.23( 1) 
Bs-A0..... —0.44(6) | —0.37(11) | —0.35(28) | —0.32(42) | —0.30(63) | —0.39(90) | —0.26(90) 
Al—A2..... —0.05(1) | —0.03( 2) | +0.21( 2) | —0.15( 3) | —0.14( 9) | —0.28(18) | —0.11(23) 
+0.01( 1) | +0.09( 2) |} +0.16( 2) 0.00( 3) | —0.11(12) | —0.08(24) 
+0.14( 1) | +0.37( 2) | +0.29( 4) | +0.09( 5) | —0.16( 2) | —0.28( 2) 
RSeBee ie elicit. +0.16( 3) | +0.21( 1) | +0.48( 2) | +0.35( 6) | +0.36( 2) | +0.33( 7) 
+0.95( 2) | +0.82( 1) | +0.73( 4) | +0.60(11) | +0.87( 9) | +0.77(15) 
K0O-K2. +1.32(2) | +1.27( 4) | +1.21( 6) | +1.33( 3) | +1.20( 8) | +1.39(11) | +1.11(15) 
K5and later} +1.37(1) | +1.56( 3) | +1.95( 1) | +2.04( 3) | +1.75( 5) | +1.89( 6) | +1.74(11) 
TABLE 2 
OBSERVED COLOR INDICES IN FIELD II 
249°, b= +1°) 
Mog 
Sp. 
8.00-8.49 8.50-8.99 9.0€-9.49 9.50-9.99 10.00-10.49 10.50-10.99 11.00-11.49 
B8-A0....... —0.14(8)| —0.36(7)| —0.10(8)| —0.14(15)} —0.12(25) | —0.16(42) | —0.19(32) 
Ny +0.11(3) | +0.28( 1) | +0.06( 2) | —0.09(21) | —0.07( 4) 
A8-F2....... +0.53(2) | +0.10(1) | +0.23(2) | +0.52( 1) | +0.40( 6) | +0.22( 7) | —0.03( 1) 
+0.39( 1) | +0.67( 3) | +0.29( 6) | —0.01( 1) 
+1.10(1) | +1.15(1) | +0.82(2) | +0.54( 2) | +0.33( 2) | +0.60( 6) | +0.13( 3) 
+1.44(4) | +1.27(4) | +1.23( 2) | 4+1.30(10) | +1.45( 9) | +1.56( 2) 
K5 and later. +2.11(1) | +1.94( 2) | +1.98( 5) | +1.77(10) | +1.81( 9) 
| 


We shall first analyze the excesses for S.A. 193. We note from Table 2 that in S.A. 193 
all excésses for the B8-AO and A1—A2 stars come out negative. Their average is not far 
from —0.20 mag.; and, in order to reduce them to zero on the average, we shall have to 
apply the maximum correction of +0.20 mag. allowed by the colors of the late-type 
stars. We therefore reach the conclusion that in S.A. 193 the B8-A2 stars down to the 
limit of 11.5 show no sign of space reddening and that interstellar reddening does not 
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he affect the colors of the stars up to a distance corresponding to a value of the modulus 
n- m — M = 11.0. The corresponding distance amounts to 1600 parsecs. This conclusion 
ors is confirmed very nicely by the observed color excesses for the B3 and B5 stars for 
which m — M < 11.0. 


TABLE 3 


AVERAGE COLOR INDICES IN S.A. 193 AND FIELD IT 
COMPARED TO NORMAL VALUES 


NorMAL CoLors S.A. 193 Fretp II 
Sp. 
49 Main Gi m=9.C to | m=11.0 to | m=9.0 to | m=10.5 to 
aa Sequence wee 10.0 11.5 10.0 11.0 
1 
0.00 | —0.11 | +0.05 | —0.08 
+0.12 | —0.08 | +0.15 | —0.09 
1) ox. +0. 32 —0.28: | +0.32 +0.22 
+0.45 +0.80 +0.39 +0.33 +0.39 +0.29 
3) re +0.70 | +1.00 | +0.75 | +0.77 | +0.68 | +0.60 
4) A) +1.05 +1.50 |. 41.25 +1.11 +1.26 +1.45 
2) 
7) 
5) 
TABLE 4 
| 
) OBSERVED COLOR EXCESSES IN S.A. 193 AND FIELD II 
| | 
BO B2 B3 | BS B8-A0 Al-A2 
Normal C.I.........} —0.50 —0.45 —0.35 | ~0.20 —0.10 +0.05 
| 
Mog S.A. 193 
$.0= 9.0. 2a +0.26(2) | +0.09(1) | —0.25(4) | —0.24( 2) | —0.29( 17) | —0.09( 3) 
a 9.0-10.0........} + .14(4) | + .15(1) — .03(9) — .25( 9) — .23( 70) — .06( 5) 
49 10.0-11.0........) + .23(2) | +0.15(1)-| +0.28(3) — .14(15) — .25(153) — .28(27) 
| Field II 
32) | | | | 
1) 
3) 
3 The more distant stars in S.A. 193 show, however, very definite signs of space red- 
a dening. We assume the following values for the mean absolute magnitudes: 
BO = —4.0 B3 = —2.2 
193 B2 = —3.0 BS = —1.6 


fat | These agree closely with those suggested by Stebbins, Huffer, and Whitford.® We find, 
2t0 | after the application of the correction of +-0.20 mag. to the excesses of Table 4, that the 
ype average excess for m — M = 11.8 equals +0.17, for m — M = 12.6, +0.37; and for 


not 5 Ap. J., 90, 209, 1939; Mt. W. Contr., No. 617. 
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m — M = 14.0, +0.40. The uncertainties for the absolute magnitudes for the small 
number of stars involved render it inadvisable to analyze these data any further, but 
we conclude that for S.A. 193 space reddening does definitely set in beyond a distance 
of 1600 parsecs and that it amounts to +0.40 mag. at a distance of 3000-4000 parsecs. 

The total photographic absorption equals probably three times the blue-red color 
excess. The color data at hand suggest that for S.A. 193 there is little or no absorption 
(less than 0.5 mag.) for distances up to 1600 parsecs. Apparently, some absorption exists 
between r =_ 1600 parsecs and r = 4000 parsecs, but the average coefficient is probably 
not more than 0.5 mag/kpc. Because of the application of the maximum allowable cor- 
rection of +0.20 mag. to all observed colors, this coefficient is an upper limit. 

The data for Field II are rather meager in comparison, since so very few early B stars 
enter in Table 4. We shall, in this case, have to depend almost exclusively on the B8-A0 
stars. If we assume that the measured color indices require no correction for systematic 
error, we can conclude from Table 4 that no perceptible reddening effects are present in 
Field II up toa distance of 1600 parsecs. If, instead, we assume that the colors require 


TABLE 5 


VALUES OF A’(m) FOR A REGION OF 100 SQUARE DEGREES 
CENTERED AT S.A. 193 


Sp. 

2.75 8.25 8.75 
35 30 87 197 
6 8 13 35 
eer 1 1 9 17 
0 1 9 12 
1 1 7 9 


the maximum admissible correction suggested by Table 3, we find that the excess for the 
B8-A0O stars at mpg = 11.0 will be, at most, +0.20 mag. The corresponding total ab- 
sorption would then be of the order of 0.6 mag. at a distance of 1200 parsecs. The colors 
for the B8—AO stars indicate, therefore, a maximum photographic absorption of 0.5 
mag/kpe for Field IT. 

Very little can be said about the photographic absorption beyond r = 1200 parsecs 
in Field II. The very few stars that are observed (see Table 2) suggest the presence of 
considerable reddening at greater distances, but it would be unwise to attempt an esti- 
mate from the few scattered colors of early B stars. 


II. Density ANALYSES FOR S.A. 193 


We found in the preceding section that the stellar distribution in S.A. 193 is apparent- 
ly unaffected by interstellar absorption to a limiting distance of 1200 parsecs. It is, 
therefore, a simple matter to compute the variations in space density with increasing 
distance from the sun. 

The data on the surface distribution of the fainter stars in an area of 6 square degrees 
have already been reproduced in Table 2, where the numbers in parentheses represent 
the total numbers of stars for separate intervals of apparent magnitude and spectral 
class. For the brighter stars we have proceeded in very much the same way as for the 
Monoceros region. The data for the surface distribution of the stars brighter than 7.0 
mag. are those of Table 7 of the preceding paper. We reproduce in Table 5 the data on 
the surface distribution of the HD stars with apparent magnitudes between 7.0 and 9.0. 
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The final smoothed log A’(m) curves were obtained from diagrams in which the log 
A'(m) values from the three different sources were plotted against the apparent photo- 
graphic magnitudes. The adopted values of log A’(m) are listed in Table 6. 


TABLE 6 
ACCEPTED VALUES OF LOG A’(m) IN S.A. 193 
SPECTRA 
Mog 

B8-A0 A1-A7 FO-F2 F5-F8 G0-G5 K0-K2 K5-M 
0.02 9.77 9.20 9.28 9.20 
C5 ee 0.35 0.00 9.50 9.50 9.45 9.68 9.05 
a: a ee 0.65 0.23 9.79 9.73 9.68 9.94 9.28 
0.95 0.47 0.07 9.95 9.93 0.20 9.52 
yey | are 1.28 0.75 0.35 0.18 0.20 0.48 9.77 
ya: ee 1.62 1.02 0.63 0.41 0.47 0.78 0.02 
Ae. Seale 1.97 1.27 0.92 0.70 0.77 1.09 0.50 
2.28 1.5% 1.23 0.98 1.10 1.41 0.95 
ae 2.80 2.03 1.78 1.47 1.67 1.90 1.63 

3.02 2.28 1.88 1.67 1.93 2.09 1.90 
3-25 2.53 1.82 1.87 2.18 2.26 2.45 
1S ay) 3.50 2.78 1.50 2.10 2.48 2.41 2.37 
TABLE 7 
ACCEPTED VALUES OF LOG A’(m)g IN S.A. 193 
GO-G5 Stars KO-K2 Srars 
Mog 
log Ame | warts | 

9.08 25 9.40 2 

ay 9.30 28 9.67 2 

9.51 33 0.02 3 

9.72 38 0.23 

9.98 40 0.47 3 

0.23 43 0.77 3 

0.53 43 1.08 3 

0.88 40 1.40 3 

1.15 42 1.67 3 

1.43 44 1.88 4 

1.65 47 2.07 5 

1.86 52 2.23 7 

2.16 52 2.37 10 


The separation of giants and dwarfs for the G—M stars was effected in the same way 
as for the Monoceros field.* The results are given in Table 7. 
For the computation of the star densities we used the values of Mo and a as listed in 
Table 11 of the preceding paper. The computed star densities, expressed in the numbers 


6 See the preceding paper. 
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of stars per 1000 parsecs? as a unit, are given in Table 8. No corrections for absorption 
were applied. 

The most striking result of Table 8 is that it is very similar to Table 13 for the Mono- 
ceros field. There are no significant differences between the density gradients in the two 
fields. We may have here a real condensation of B8-AO stars, which reaches its maximum 


TABLE 8 


SUMMARY OF STAR DENSITIES IN S.A. 193 
(Numbers of Stars per 1000 Parsecs*; Apg=0.0 Mag/Kpc) 


| SPECTRA 
d 
| B8-A0 A1-A7 F0-F2 FS-F8 | gG0-gM 
16 .30 0.9 1.4 .26 
0.6 1.0 .25 
TABLE 9 
COUNTED NUMBERS OF STARS FOR AN AREA OF 
50 SQUARE DEGREES CENTERED ON FIELD II 
Mog 
Sp. T 
7.25 7.75 8.25 8.75 
B8—A0. 4 24 31 46 
0 3 6 20 
0 2 5 9 
3 4 13 15 
0 1 3 1 


density at between 200 and 500 parsecs from the sun, where the star density exceeds 
that for the vicinity of the sun by more than 30 per cent. At a distance of 100 parsecs 
from the sun the density of the B8—A0 stars is, in contrast to our results for Monoceros D, 
about equal to that for the immediate vicinity of the sun. The slight negative gradients 
for the densities of the A1—A7 stars and the red giants and the steep negative gradients 
for the F’s are again in marked contrast to each other. 


III. Density ANALYSES FOR FIELD IT 


The density analyses for Field II proceed along the same lines as those for S.A. 193 
and the Monoceros field. Again we use the data of Table 7 of the preceding paper as 


i 
} 
2 
3 
| 
| | | 
i 
° 


STELLAR DISTRIBUTION 307 
n representative of the distribution of the stars brighter than m = 7.00. The distribution 
of the stars between the seventh and ninth magnitudes is shown in Table 9. The counts 
)- refer to an area of 50 square degrees, centered upon Field IT. In the choice of this area, 
) care was taken not to have it cover the rich regions near Carinae; the square lies be- 
n tween galactic longitudes 246° and 253° and latitudes —3° and +4°. 
TABLE 10 
ACCEPTED VALUES OF LOG A’(m) FOR FIELD II 
SPECTRA 
Mog 
| B8-AO A1-A7 F0-F2 F5-F8 G0-G5 K0-K2 K5-M 
0.05 9.78 9.20 9.32 9.18 9.38 8.75 
sf Se | 0.38 0.00 9.50 9.54 9.40 9.64 9.05 
6.25. 6445, | € 70 0.22 9.78 9.76 9.62 9.92 9.30 
1.00 0.44 0.04 9.98 9.88 0.21 9.55 
4.25 0.68 0.30 0.20 0.18 0.51 9.80 
1.57 0.96 0.60 0.45 0.52 0.83 0.08 
8.25... 1.80 1.22 0.90 0.70 0.90 1.14 0.38 
2.00 1.50 1.20 0.95 1.22 1.45 0.75 
2.20 1.80 1.42 1.20 1.40 1-22 
Ree | 2.40 2.08 1.67 1.45 1.58 1.92 1.50 
| 2.61 2.38 1.90 1.70 1.74 2.10 1.90 
1a. as | 2.82 2.70 2.10 1.97 1.90 2.23 2.20 
TABLE 11 
ACCEPTED VALUES OF LOG A’(m)yg FOR FIELD II 
GO0-G5 Stars KO-K2 Srars 
Mpg 
lox Amo | Dwarts | | 
9.23 32 9.63 2 
9.42 38 9.90 4 
yc Se 9.64 42 0.20 3 
9.94 42 0.50 3 
0.31 38 0.82 3 
S25... 0.73 32 1.13 3 
1.06 30 1.44 3 
s 9.25... 1.18 40 1.70 3 
Ss 9.75... 1.26 53 1.97 3 
), 1.18 74: 2.08 5 
s 


Table 10 gives the final smoothed values of log A’(m) from the three different sources 
combined. The separation of giants and dwarfs for the G-M stars was effected in the 
same manner as for the Monoceros field. The results are given in Table 11. 

3 The star densities were again computed with the aid of the values for Mo and a listed 
s in Table 11 of the preceding paper. The results are in Tables 12 and 13. The entries in 


“gd 
A Ae 
— 
| 
7 
4 


308 ; BART J. BOK AND FRANCES W. WRIGHT 


Table 12 are the densities computed for the case of no absorption. In Table 13 the densi- 
ties were corrected for the absorption of 0.5 mag/kpc suggested by the color data. 

A comparison between Table 13 of the preceding paper and Tables 8 and 13 of the 
present paper shows that there exists a great deal of similarity between the star densities 


TABLE 12 


SUMMARY OF STAR DENSITY IN FIELD II 
(Numbers of Stars per 1000 Parsecs*; No Apg Applied) 


SPECTRA 
y | 

B8-A0 A1-A7 F0-F2 F5-F8 gG0-gM 

0.12 0.45 1.4 3:5 0.30 

ris 27 0.8 1.6 | 

TABLE 13 


(Numbers of Stars per 1000 Parsecs*; Corrected for Apg=0.5 Mag/Kpc) 


SUMMARY OF STAR DENSITY IN FIELD II 


SPECTRA 
| B8-A0 A1-A7 F0-F2 F5-F8 gG0-gM 
0.13 0.47 1.4 3.5 0.30 
.14 -35 152 2.4 
| .18 30 1.0 1.8 .30 ( 
in Monoceros D, S.A. 193, and Field IT. In all three fields the B8—A0 stars show an initial ( 
increase in average density, large enough, probably, to be real, followed by a decrease. | 
But, whereas in Monoceros D the density, at 1000 parsecs is decidedly less than that for 
the vicinity of the sun, we find for S.A. 193 and Field II an average star density for the 
B8-—AO stars again almost equal to that for the solar neighborhood. t 
The red giants in all three cases show almost constant density over the first 500 par- c 


secs, with an indication of a gradual drop in average density; the gradient is least for 
Field II, steepest for S.A. 193, and intermediate for the Monoceros field. The F stars 
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behave similarly in the three fields; there is apparently no escape from a rapid decrease 
in the star density for these stars. 


IV. THE DISTRIBUTION OF THE BO-B5 Stars 


It is of interest to compare the distribution of the early B stars for the field in Mono- 
ceros, S.A. 193, and Field II. The relevant data are reproduced in Table 14. The Mono- 
ceros numbers represent the total for regions C, D, and F of the preceding paper. 

The large number of early-type stars in S.A. 193 is quite noticeable. According to 
Seydl’s tabulations,’ we observe, on the average, two BO-B2 stars and seven B3-B5 
stars between mp, = 4.5 and mpg 5.5 per 500 square degrees near the galactic circle. In the 
absence of absorption, and for constant average space density, we should then expect 
about three BO-B2 stars and ten B3-B5 stars between m,, = 8.0 and my, 9.0 for the area 
of 6 square degrees covered by S.A. 193. Since absorption is present in S.A. 193 for dis- 
tances beyond 1600 parsecs, we conclude from Table 14 that the observed numbers of 
BO-B5 stars are compatible with a constant star density equal to that for the vicinity 
of the sun up to a limit of 2500 parsecs. 


TABLE 14 
OBSERVED DISTRIBUTION OF EARLY B STARS 


MOoNOCEROS S.A. 193 II 

(67 SQUARE (6 SQUARE (6 SQUARE 

DEGREES) DEGREES) DEGREES) 

Myg 
BO-B2 B3-B5 | BO-B2 B3-B5 BO-B2 B3-B5 

1 5 3 6 1 
9.0-10.0...... 5 23 5 18 0 1 
10.0-11.0...... Z 33 3 18 1 0 


Quite the reverse is true for the fields in Monoceros and Field II. Even if we consider 
that for Monoceros Ay, may be as high as 0.8 mag/kpe and that in Field IT Ayg = 
0.5 mag/kpc, we find that the numbers predicted for constant density are several times 
larger than the observed ones; and we conclude that the reality of negative density gra- 
dients for the BO-BS stars at / = 185° and / = 249° is established. These results, together 
with the observed high frequency of early B stars in the region of » Carinae’ at / = 254°, 
make it seem likely that the region of constant density of the early B stars has a sharp 
boundary near / = 252°. 

When we compare these results with those obtained for some other sections of the 
Milky Way, especially with those investigated by Schalén, Wernberg, and Wallenquist, 
we find that the star densities drop sharply as we move away from the sun for / = 73° 
and/ = 135°. In Cygnus we find a very slight decrease at / = 40°. In the direction of the 
Sagittarius and Scutum clouds, low densities are observed at 1000 parsecs from the sun; 
but for the Scutum region, the space densities of the early B’s at 2000 parsecs are again 
quite high. 

V. THE CHANGES IN THE LUMINOSITY FUNCTION 


The values of log ¢’(M) were computed by the same procedure as was followed in 
the construction of Table 14 of the preceding paper. Because of the similarity in density- 
curves for the various spectral classes in the three regions, it is not surprising that there 


. ™ Pub. Nat. Obs. Prague, No. 6, 1929. 8 Harvard Reprints, No. 77, 1932. 
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appears to be considerable resemblance between Tables 15 and 16 of the present paper 
and Table 14 of the Monoceros paper. 

It may be of interest to compare these results with data in the changes in the general 
luminosity function that may be derived for other sections of the Milky Way from spec- 
tra and magnitudes determined by Swedish astronomers. Two investigations by Scha- 


TABLE 15 
CHANGES IN THE GENERAL LUMINOSITY FUNCTION FOR S.A. 193 
(Values of log ¢’(M)) 
d (PARSECS) 
My 
Van Rhijn 50 100 200 300 500 
8.38 8.53 8.50 8.52 8.49 8.43 
8.90 8.85 8.82 8.78 8.69 
9.04 9.15 9.09 9.03 
TABLE 16 
CHANGES IN THE GENERAL LUMINOSITY FUNCTION FOR FIELD II 
(Values of log ¢’(M)) 
d (PARSECS) 
My. 
Van Rhijn 50 100 200 300 500 
8.38 8.57 8.53 8.59 8.62 8.54 


lén®’ provide the necessary data for fields in Cygnus and in Auriga. For a region in 
Cepheus we have the material of Wernberg,!! and for the Sagittarius cloud we have made 
use of Wallenquist’s” spectral data. Since for most of these fields we have little informa- 
tion on the densities of the A stars within 200 parsecs of the sun and since existing sur- 
veys tell us little about the densities for the F stars beyond this distance, we shall limit 
ourselves here to the computation for every field of the value of log ¢’(M) for the dis- 
tance of 200 parsecs from our sun. 


Uppsala Ann., Vol. 1, No. 4, 1941. 
2 Bosscha Ann., Vol. 5, No. 5, 1939. 


9 Uppsala Medd., No. 37, 1928. 
10 JbiJ., No. 55, 1932. 
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For the region in Cygnus, we deduce from Schalén’s colors, and from the photoelectric 
colors of Stebbins, Huffer, and Whitford,'* that at / = 41°, b = +1° to +3°, the ab- 
sorption is very nearly equal to A,, = 1.0 mag/kpc. For the region in Cepheus, / = 70° 
to 75°, b = 0° to —3°, studies by Wernberg and Miss Risley" indicate an average value 
Ay, = 0.33 mag/kpc. For the Sagittarius region we assume from the colors of Stebbins, 
Huffer, and Whitford that in Wallenquist’s Region II the absorption is represented by 
Aig = 1.0 mag/kpc. 

The results are collected in Table 17. The figures in the last column represent the 
general luminosity function of van Rhijn,'° which represents the best available reference 
curve for the vicinity of the sun. 

We have not included in our survey Schalén’s material for the Scutum cloud,'® but 
from the data at hand it seems that the values of log y’(M) for Scutum resemble those 
for the Sagittarius region. 

TABLE 17 
THE GENERAL LUMINOSITY FUNCTION AT 200 PARSECS FROM THE SUN FOR 


VARIOUS DIRECTIONS IN THE GALACTIC PLANE, COMPARED WITH 
VAN RHIJN’S STANDARD FUNCTION VALUES OF LOG ¢’(M) 


FIELD AND GALACTIC LONGITUDE 
M 
Cygnus Cepheus Auriga | Monoceros| Field II S.A. 193 | Sagittarius | Van Rhijn 
40° 7 135° 185° 249° 261° 332° Standard 
8.67 8.45 8.30 8.26 8.25 8.18 |... 
Un) aie ee 8.88 8.64 8.67 8.56 8.59 8.52 8.23 8.38 
PROBE. 9.25 9.17 8.96 9.04 9.10 9.03 8.65 9.04 
9.40 9.35 9.13 9.17 9.26 9.18 
9.49 9.48 9.26 9.26 9.36 9.28 8.95 9.47 
9.57 9.43 9.30 9.47 9.36 9.04 9.56 
9.60 9.58 9.49 9.34 9.50 9.40 


The most disturbing result of Table 16 is that the general luminosity function ap- 
parently undergoes considerable variations in its shape over comparatively small dis- 
tances in the galactic plane. The existence of such variations has been suggested pre- 
viously. Maxwell!’ gave evidence for variations in the shape of the general luminosity 
function for the direction of the Cygnus cloud. Later, Seares'* gave evidence for an ex- 
cess of highly luminous stars in the vicinity of the sun, and the possible effects of such 
changes were discussed by Bok.'® More recently, Miller and Hanna”? have computed 
the variations for the Cygnus cloud, for which they find an excess of A and K stars and 
a deficiency of B and F stars, as compared to the standard function. 

Luyten”! and Schilt” have given some results that show that our knowledge of the 


13 Op. cit. 

14 Ap. J., 97, 277, 1943. 

15 Groningen Pub., No. 47, Table 6, 1936. ; 

16 Uppsala Medd., No. 61, 1935. 7 Lick Obs. Bull., No. 390, 1927. 

18 4p, J., 74, 312, 1931; Mt. W. Contr., No. 437. 

19 Harvard Circ., No. 371, 1931. 

20 4.J., No. 1139, 1942; Swasey Obs. Reprints, No. 8. 

2 Ann. N.Y. Acad. Sci., 42, 201, 1941. 2 Ann. N.Y. Acad. Sci., 42, 259, 1941. 
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general luminosity function for absolutely faint stars is yet far from complete. The 
frequencies of stars with M 2 +6, however, do not concern us especially in the analysis 
of the stellar distribution for fields in low galactic latitudes. 

Our results in Table 17 indicate that the variations are quite erratic and large, even 
for absolute magnitudes M = 0 to +4, for which interval it had generally been supposed 
that the variations would be mild. We shall apparently have to be very careful in inter- 
preting the densities derived from general star counts. These densities are at best repre- 
sentative for the most “effective” absolute magnitude. It is well to remember that this 
“effective”? magnitude becomes fainter as we pass to fainter apparent magnitudes. The 
observed variations in the general luminosity function render it all the more necessary 
that we should supplement our star-count analyses with studies based on spectral dis- 
tributions. 

In a way, the most striking result of Table 17 is represented by the column for the 
Sagittarius cloud. We are likely to think of the direction to the center of galaxy as one 
of very high star density; and instead we find that at 200 parsecs from the sun the values 
of log ¢’(M) for 1 = 332° are smaller than those for any other galactic longitude, includ- 
ing the anticenter directions. The average density at 200 parsecs for / = 332° is just 
about equal to one-quarter of that at the same distance for / = 40°. The distance of 200 
parsecs to which the results of Table 17 apply is so small that it is unlikely that uncer- 
tainties in the absorption characteristics to that distance might account for the deficien- 
cies at / = 332°. We note, moreover, that the result is confirmed in a general way by the 
data for the Scutum cloud. 

Apparently, there is no escaping the conclusion that for the direction of the Scutum 
and Sagittarius clouds the star densities drop at first and recover again at larger dis- 
tances. This result has already been foreshadowed by an earlier investigation, based on 
the analysis of general star counts'**, Studies based on spectra and colors of faint stars 
for the bright patches in the central directions of the Milky Way not only should yield 
important information on the stellar distribution in the more populated distant central 
regions of our galaxy but should further contribute important data bearing on local 
structure. 


VI. RESULTS AND CONCLUSIONS 


The investigation deals with the stellar distribution in two southern fields: one at 
1 = 261°, b = +1° (S.A. 193); the other at / = 249°, 6 = +1° (Field II). In each case 
the area covered is 6 square degrees; the basic data are given in Tables 1 and 2. The main 
results of the analysis are the following: 

1. Selected Area 193 is apparently free from obscuration to a distance of 1600 par- 
secs; beyond that the absorption is represented by Apg = 0.5 mag/kpc to r = 3000 
parsecs, 

2. In Field II a value A,g = 0.5 mag/kpc to r = 1200 parsecs accounts for the ob- 
served color excesses; there is no information for greater distances. 

3. The density distribution for S.A. 193 (see Table 8) are similar to those of the 
Monoceros field. We add the comments: 

a) The excess density for the B8—AO stars persists to r = 500 parsecs; there is no 
marked decrease between r = 500 and r = 1000 parsecs; 

b) The densities for the Al—A7 stars drop to 35 per cent of the value for the vicinity 
of the sun at r = 500 parsecs; 

c) The FO-F2 and F5-F8 stars show the same steep mapetirs gradients as were ob- 
served for the Monoceros field; 

d) The gG0-gM stars yield slowly decreasing star dinalties. 


23 The Distribution of the Stars in Space, pp. 120 and 123, Chicago: University of Chicago Press, 1937. 
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4. The density distributions in Field II (see Table 13) are characterized by— 

a) High densities for the B8-A0O stars, with the density at r = 800 parsecs equal to 
that near the sun; 

b) A slight negative gradient for the Al—A7 stars; 

c) Steep negative gradients for the F stars; 

d) Slowly decreasing densities for the red giants. 

5. For the fields in Monoceros and for Field II the predicted numbers of BO-B5 stars 
(computed from the observed averages for m = 5.0 along the galactic belt on the as- 
sumption of constant average space density) fall far above the numbers actually ob- 
served. This suggests a rapid decrease in star densities for the early B stars in the direc- 
tions / = 185° and/ = 249°. For S.A. 193 (as for the neighboring region near y Carinae) 
we find, in contrast, fairly high and constant average densities for the early B stars. It 
is of interest to compare these results with those obtained for other sections. The be- 
havior for the Monoceros field and Field II appear to be characteristic for the whole re- 
gion between / = 70° and / = 252°. In Cygnus the densities behave not unlike those for 
the Carina-Crux-Centaurus region. In Sagittarius and Scutum we have apparently an 
initial steep decrease, followed by a rapid rise in star density at greater distances. 

6. Tables 15 and 16 show that for Field II and S.A. 193 the dependence of the shape 
of the general luminosity function on distance is again quite pronounced. Table 17 shows 
that large variations exist for almost all directions of the Milky Way. The behavior of 
the Sagittarius region is especially noteworthy. 
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AN INVESTIGATION ON DIFFERENTIAL GALACTIC ROTATION 


Paris Pismis AND AGUSTIN PRIETO 
National Astrophysical Observatory, Tonanzintla, Mexico 
Received January 18, 1945 


ABSTRACT 


Solutions for 7A and K have been carried out, using Oort’s equation and the radial velocities of 707 
stars of spectral types O-B9, common to the catalogue of radial velocities of O- and B-type stars by Plas- 
kett and Pearce and the catalogue of 1332 stars of Stebbins, Huffer, and Whitford. It appears that A 
does not show any variation with distance; its mean value is 0.011 + 0.002 km/sec-pc, much smaller 
than the generally adopted value of 0.017 km/sec: pc. 

Solutions carried out by dividing each distance group into two subgroups, one containing stars of 
spectra] types O-B2 and c stars and the other, spectral types B3-B9, give A = +0.016 and A = 0.009 
km/sec-pc, respectively. It is therefore concluded that the small value of A obtained from the general 
solution is due to an overestimation of the distances for the B3-B9 stars. This overestimation is caused, 
very probably, by the adoption (by Stebbins, Huffer, and Whitford) of too bright absolute magnitudes 
for those types. 

The K-term practically vanishes in the nearest group. By dividing the spectral subgroups into 
quadrants in longitude it is shown that the disappearance of the K-term for the nearest stars is a real 
effect, whereas the large K-term obtained hitherto for the nearest early-type stars is probably due to the 
nonuniform distribution in longitude or to loca] irregularities. 


Numerous investigations have been made to determine the coefficient A since the 
time when Oort first introduced the formulae for differential galactic rotation. A reliable 
value of this coefficient is desirable, since it is used quite often as a means for estimating 
distances. 

Any determination of A requires a knowledge of the mean distances of the stars on 
which the determination is based; for, in the well-known formula for differential rotation, 


rAsin2(l—h)co? b+ K=V/, (1) 


where V denotes the radial velocity of a star after the solar-motion component is elimi- 
nated, A enters with a factor r, the distance of the star. In general, for a group of near-by 
stars for which reliable distances are available, the term 7A will be indeterminate, and 
hence the value of A will not be reliable; while for stars far enough to give reliable values 
for A one is confronted with another difficulty, namely, the lack of reliable distances. The 
distances for the farther groups are, in general, photometric. Distances of this kind and, 
hence, A derived from them will be affected by the adopted value of the absorption 
coefficient. 

Since the patchy distribution of the absorbing matter in the Galaxy is being empha- 
sized more and more, the question arises as to how far a mean absorption coefficient rep- 
resents the effect of the nonuniform distribution of the absorbing matter. R. E. Wilson 
has shown in a recent paper! that the nonuniformity of the absorbing matter will give 
rise to a spurious decrease of the coefficient A with increasing distance. If one were to 
force the derived values of A for different distance groups to agree by adjusting a mean 
coefficient of absorption, the result would be a value of A much higher than the true one. 

The best value of A, therefore, can be obtained by using distances corrected for non- 
uniform absorption. In an earlier paper, one of the authors? made a solution of equation 
(1) with the aid of the radial velocities of O- and B-type stars determined by Plaskett and 
Pearce’ and the distances given by Stebbins and Huffer.‘ Least-squares solutions were 


1Ap. J., 92, 170, 1940. 

2 Paris Pismis, Pub. Istanbul Univ. Obs., No. 10, 1938. 

3 Pub, Dom. Ap. Obs., Vol. 5, No. 3, 1933. 4 Pub. Washburn Obs., Vol. 15, Part 5, 1934. 
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carried out for four distance groups. The material was, however, too scanty and not well 
distributed in galactic longitude. 

At present the more extensive catalogue of Stebbins, Huffer, and Whitford offers us 
the possibility of extending the earlier investigation. The data for the present solution 
consist of stars common to the catalogue of 1010 stars of types O and B by Plaskett and 
Pearce’ and the catalogue of 1332 stars by Stebbins, Huffer, and Whitford.® 

The distances given in the latter catalogue are based on two assumptions: first, that 
of a constant ratio of general to selective absorption (A/E = 7) and, second, that of an 
average absolute magnitude for each subtype. The recent investigation on the six-color 
photometry of stars by Stebbins and Whitford’ seems to support the high ratio of general 
to selective absorption; they further show that this ratio does not vary with the region of 
the sky considered. This implies that very probably no systematic errors depending on 
distance or direction will be introduced in the distances based on. the assumption of the 
mentioned ratio. 

The material consisting of 707 stars was divided into five distance groups. An equa- 
tion of condition based on equation (1) was formed for each star. The solar motion was 
assumed to be 20 km/sec in the direction a = 271° and 6 = 28°, and its component was 
eliminated by the use of the graph by Pearce and Hill.* The longitude /) was assumed to 
be 325°. The equations were solved for 7A and K in five distance groups by the method of 
least squares. The results with their probable errors appear in Table 1. 


TABLE 1 
SOLUTION FOR 7A AND K USING 707 O- AND B-TYPE STARS 
Distance Limits Mean No. of rA K A 
Group (Pc) Distance | Stars (Km/Sec) (Km/Sec) (Km/Sec -Pc) 
| r<250 182 128 + 1.36+0.8 +0.36+6.3 +0.007 +0.004 
2......| SOO 358 292 4.83+0.8 +0.81+0.5 .013+ .002 
500<r< 750 589 5.62£1.2 +2.81+0.8 .002 
750<r< 1000 855 57 10.40+2.7 —2.72+1.9 .003 
ae 1000<r 1443 53 +25.744+2.9 —7.95+2.3 +0.018+0.002 


A does not show any variation with distance if we leave out the last group (this omis- 
sion will be justified in the discussion to follow). The weighted mean over all groups gives 
A = +0.011 + 0.001 km/sec - pc. Excluding the first group, for which 7A may be in- 
determinate, the weighted mean A is 0.012 + 0.003. Both of these values are much small- 
er than the usually accepted one, 0.017 km/sec - pc. ‘ 

The K-term is very small for the nearest two groups. This result was a, real surprise, 
for solutions made hitherto had always yielded a K-term of the order of 4-5 km/sec for 
the nearest stars of types O and B. Moreover, as has been emphasized by one of the au- 
thors,? some earlier investigations indicate that the A-term decreases with distance. In 
general, there appears a decrease of the K-term in the farthest group, but the variation 
is by no means marked and smooth. The number of stars in the nearest three groups is 
sufficiently large (and the probable errors small enough) to make the determination of 
K and also of 7A fairly reliable. 

The material of the five groups is quite uniformly distributed between galactic longi- 


‘tudes 10° and 190°, whereas there exist large gaps at longitudes corresponding to southern 


declinations. Taking only the stars between / = 10° and / = 190°, we have made a new 
solution of equation (1) for the five distance groups. (This interval covers a complete 


5 Pub. Dom. Ap. Obs., Vol. 5, No. 2, 1931. 7 Ap. J., 98, 20, 1943. 
® Ap. J., 91, 20, 1940. +’ Pub. Dom. Ap. Obs., Vol. 6, No. 4, 1931. 
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wave of the differential rotation.) By this procedure one can avoid to some extent the 
effects of the nonuniform distribution of the velocities. The results obtained in this way 
should not differ essentially from those for the complete groups, provided there are no 
asymmetries (such as stream motion) in the groups considered. The results of the new 
solutions are given in Table 2. Again A shows no variation with distance. The weighted 
mean value of A is 0.013 + 0.001, a value smaller than that found in previous investi- 
ations. 

Z To give more weight to the determination of A, we have made a new solution, assum- 
ing K = 0. The results are given in Table 3. The weighted mean value of A is +0.012 
+ 0.001 km/sec - pc. 


TABLE 2 
SOLUTION USING STARS IN LONGITUDES 10°-190° 
No. of 7A K A 
Stars (Km/Sec) | (Km/Sec) | (Km/Sec*Pc) 

Se ne 177 105 + 1.88 + 0.50 +0.011 

2 361 268 5.63 + 0.92 .016 

589 164 4.81 + 1.30 .008 

868 52 12.01 — 0.01 .014 

1484 46 +21.22 —14.30 +0.014 
TABLE 3 


DETERMINATION OF A, ASSUMING K =0 (707 STARS) 


Mean oA A 
(Km/Sec) | (Km/Sec* Pc) 
(Pc) 
182 + 1.38 +0.007 
589 6.32 O11 


It thus seems certain that the value of A obtained from the present material (spectral 
types O-B9 taken together) is approximately 0.012 km/sec + pc. Two possible explana- 
tions can be given for the small value which has been derived for A. Either it is real, 
which does not seem probable, or it is due to the circumstance that the distances given 
by Stebbins, Huffer, and Whitford are too large. The latter explanation may be the valid 
one for two reasons: The adopted ratio of general to selective absorption is too small and 
the mean absolute magnitudes adopted are too bright. 

In the Washburn catalogue of 733 stars the assumed ratio of general to selective ab- 
sorption was only 2. On the other hand, the absolute magnitudes adopted were too faint, 
compared with the recent catalogue of 1332 stars. In addition to this, the stars with c 
characteristics were not differentiated from the normal stars. Stebbins, Huffer, and 
Whitford assert that the increase of this ratio will be compensated by the increase of the 
absolute luminosities. To check on this we have made a comparison of the distances given 
by the two catalogues mentioned. A total of 643 stars is common to the two catalogues; 
they have average distances of 368 and 476 pc according to the early and later catalogues, 
respectively. It thus appears that the distances of the recent catalogue are, on the aver- 
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age, 26 per cent larger than those of the earlier one. Most of this difference is probably 
caused by the separation of the c stars, which are now assigned brighter absolute magni- 
tudes. On the scale of the earlier distances the mean value of A, given in Table 1, would 
become 0.014 km/sec - pc, still too low compared to the generally adopted value. 

To investigate the question further, as well as the peculiarity shown by the K-term, 
each of the distance groups was divided into two subdivisions, one (a) containing stars of 
types O-B2 and stars showing c characteristics and the other (6) containing stars of 
types B3-B9. This is essentially a subdivision with regard to absolute magnitude. The 
results of the least-squares solutions for 7A and K and their probable errors are given in 
Table 4. All stars in distance group 5 are either types O-B2 or c stars; therefore its cor- 
responding subgroup for B3—B9 types does not exist. The fourth subgroup for types B3- 


TABLE 4 
SOLUTIONS OF SPECTRAL SUBGROUPS 
No MEAN | 
Dis- rA K A 
TANCE (Ka/SeEc) (Ka/SeEc) (K/Sec - Pc) 
TARS 
| (Pc) | 
| 
Subdivision a: Types O—B2- and c Stars 
19 179 — 0.60+ 1.6 +2.994+1.1 —0.003+0.009 
53 371 + 7.62+ 2.1 +2.10+1.9 + .020+ .005 
Ee, os 81 520 + 9.53+ 1.3 +3.72+0.8 + .018+ .003 
48 894 +12.67+ 2.8 +0.23+2.0 + .013+ .003 
53 1443 +25.74+ 2.9 —7.95+2.3 +0.018+0.002 
Subdivision 6: Types B3-B9 
109 174 + 1.67+ 0.4 +0.0 +0.3 +.0.009+0.002 
239 355 + 4.214 0.9 +0.4 +0.2 + .012+ .003 
96 580 + 2.83+ 1.9 +0.6 +1.2 + .005+ .003 
9 757 — 8.14+15.3 —1.2 +8.9 —0.011+0.020 


B9 will not be included in the discussion, for it has too few stars to yield reliable values 
for 7A and K. 

The coefficient A is systematically larger for subgroups a compared with that for sub- 
groups b. We see, then, that the apparent increase of A in the last group of Table 1 is due 
to the fact that all stars are of types O-B2 or show c characteristics and hence yield a 
larger value of A than the average over all types. The weighted mean A for the five 
groups of subdivision a gives A = 0.016 + 0.002 km/sec - pc, omitting the first group, 
A = 0.017 + 0.002. These values are quite near to those obtained hitherto. The weighted 
mean over the first three groups for the B3—B9 stars is A = 0.009 + 0,001. We therefore 
concludé that the small value of A of the general solution is due to the later-type stars. 
It seems, therefore, that the assumed distances for the stars B3—B9 are too large, which 
results in a low value for A. The overestimation of the distances cannot be caused by an 
erroneous ratio of géneral to selective absorption, for the latter would affect all stars 
equally.® It is therefore highly probable that the adopted absolute magnitudes of B3-B9 


®It seems unlikely that the normal colors of B3-B9 stars assumed by Stebbins, Huffer, and Whit- 
ford involve errors large enough to account for the overestimation of the derived distances. 
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stars are too bright. The correction to be applied to the mean absolute magnitudes for 
the average of B3-B9 stars adopted by Stebbins, Huffer, and Whitford® is of the order of 
1.2 visual magnitude. This quantity was obtained by assuming that the mean A of the 
B3-B9 stars is of the order of 0.016 km/sec + pc. (The correction will amount to 1,4 
visual magnitudes if we assume A = 0.017.) 

The K-term of Table 1 presents an interesting aspect. It has been mentioned earlier 
that in solutions made hitherto the nearest early-type stars yielded a large positive K- 
term. In the present solution the K-term of the nearest group practically vanishes. The 
K-term for stars between longitudes 10° and 190° also shows the same peculiarity as that 
for the entire group. We shall discuss the different groups in some detail later. 

It is worth while to call attention to the fact that in Table 4 the K-term for subdivision 
a is systematically larger than that for subdivision 6 by about +2.3 km/sec on the aver- 
age for the first three groups. If significance could be attached to the K-term of all the 


TABLE 5 


AVERAGE RADIAL VELOCITIES, IN KM/SEC, CORRECTED FOR SOLAR MOTION 
FOR QUADRANTS IN GALACTIC LONGITUDE 


SuppIvision a: Types O—B2 Anp c STARS Suppiviston 6: Types B3—B9 
Group Longitude Limits 
325°-55° 55°-145° 145°-235° 235°-325° 325°-55S° 55°-145° 145°-235° | 235°-325° 

— 0.1(3)* + 2.0(5) |+ 3.8(5) |+ 4.1 (6)}+1.5 (30)|—1.1 (49)|+1.8 (26)|+3.8 (4) 
+ 5.1 (6) — 2.3 (22))+ 7.8 (23) 8.5 (2)} 2.7(75)|} 1.2 (82)|+3.3 (82) (0) 
< 9.2 (20) — 1.5 (38)/+ 9.3 (22)/4+19.0(1)| 4.0 (38)} 0.4 (29)|/+0.0 (29) (0) 
eee. + 8.9 (18) +13.8 (23)|— 0.1 (7) (0)|+3.5 (2) |—8.9 (2) |—0.0 (5) (0) 
Differen- 

tial ro~ 

tation 

accord- 

ing to 

Oort. . >0 <0 >0 <0 >0 <0 >0 <0 


* Numbers in parentheses represent number of stars. 


different subgroups, this difference of +2.3 km/sec could be interpreted as a physical 
red-shift, probably a gravitational red-shift which will be larger for the more luminous, 
hence the more massive, stars of subdivision a in comparison to the less luminous, and 
hence the less massive, stars of subdivision 0. 

In order to investigate the behavior of the K-term more closely and to account for the 
disappearance of this term for the nearest group of stars, we have made a Study of the 
radial velocities in longitude quadrants. Each of the subgroups appearing in Table 4 was 
divided into four quadrants, which were chosen in such a way that at their limits the 
effect of the differential rotation on the radial velocities is zero (those limits are at 325°, 
55°, 145°, and 235°, respectively). For each of the quadrants the mean radial-velocity 
with the corresponding number of stars is given in Table 5. To compare the present re- 
sults with a solution yielding a large K-term, we have given in Table 6 the mean radial 
velocities in four quadrants of the nearest group of Solution I in an earlier investigation 
mentioned before.” The reasons for this choice are, first, that the earlier groupings were 
also based on distances (Stebbins and Huffer) and, second, that the stars could easily be 
identified. 

In the fourth quadrants of all groups in Tables 5 and 6 the radial velocities of the stars 
are either unknown or are too few to give reliable mean values. Our discussion, therefore, 


. | | | | | | | | 
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will be based only on the three first quadrants. The first group of subdivision a and the 
last group (fourth) of subdivision b have too few stars to afford mean values reliable for 
discussion. 

It is well known that the K-term represents the shift in radial velocity of the double 
wave of the differential rotation. The K-term has significance as a general effect only 
when this shift is systematic. This means that among the four quadrants considered in 
Table 5 the mean radial velocities in quadrants 1 and 3 should agree sensibly; a similar 
agreement should exist between quadrants 2 and 4. The mean radial velocity of two con- 
secutive quadrants (or of all four quadrants) should give an estimate of the K-term. 
(Owing to the high concentration of the early-type stars to the galactic plane, we have 
assumed their mean galactic latitudes to be zero in discussing the quadrants.) 

Among the subgroups of B3-B9 stars the one that comes nearest to satisfying the 
above condition is the nearest group, the group most essential for our discussion of the 
K-term. The double wave appears very clearly in the first three quadrants. It does not 
show any systematic shift. The mean value of the radial velocities in the first two quad- 
rants practically vanishes, in agreement with the very small K-term obtained from the 
least-squares solution. There is no doubt, therefore, that the K-term of this group is a 
general effect and not due to irregularities in distribution. 


TABLE 6 
AVERAGE RADIAL VELOCITIES, IN KM/SEC, CORRECTED FOR SOLAR MOTION 
FOR QUADRANTS IN GALACTIC LONGITUDE 
(47 Stars; Mean Distance, 189 Pc) 


Limits 


325°=55° | 55°-145° 145°-235° 235°-325° 


43.7 (11)* | —2,8(14) | +750) | 1.7@) 


* Numbers in parentheses represent number of stars. 


In Table 6 we see, on the other hand, that the large positive K-term which the least- 
squares solution of this group has yielded (+5.2 + 1.6 km/sec; mean distance of group, 
189 pc) is brought about by the excess of positive velocities in the third quadrant, which 
has the largest number of stars. If the number of stars in the other quadrants were com- 
parable, one could decide whether the large K-term is due to a systematic shift or is 
caused by the streaming tendency of the stars in the third quadrant. But, as it stands 
now, the K-term in this group is a local effect. Out of the 47 stars of this group, 43 are in- 
cluded in the material of the present investigation. It is interesting to mention that these 
43 stars are scattered over the three nearest groups of Table 1. Only 15 stars are of spec- 
tral types O-B2, and two others show c characteristics. 

Of the remaining groups in Table 5, distance group 2 of B3—B9 stars indicates that 
probably its K-term is real—the mean velocities in the quadrants show a systematic 
positive shift. The third group shows asymmetry. Of the O-B2 and c-star groups, the 
second and the third distance groups indicate a systematic shift toward positive veloci- 
ties. The fourth group shows an asymmetry and the fifth one a systematic shift toward 
negative velocities, pointing out that the A-term of this group is very probably a real 
effect. 

In conclusion we may emphasize again that the K-term has significance as a general 
characteristic for a group only when the double wave of the velocities shows a systematic 
shift. It will not be safe to consider the K-term as a general effect unless the material on 
which it is based is well distributed in galactic longitude and is free of local effects, such 
as group motions. 
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ON THE RADIATIVE EQUILIBRIUM OF A STELLAR ATMOSPHERE. VI 
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ABSTRACT 
In this paper the solutions to the problem of line formation in stellar atmospheres, in the third approxi- 
mation, are tabulated for various values of the parameters which are involved and at close enough inter- 


vals to allow interpolation. 
Numerical forms of the solution for the problem of the radiative equilibrium of a planetary nebula 


are also obtained. 


This paper is of the nature of a supplement to two earlier papers’ of this series and pro- 
vides some further information relating to the problems of the formation of absorption 
lines in stellar atmospheres and the radiative equilibrium of a planetary nebula. 


I. THE FORMATION OF ABSORPTION LINES; THE TABULATION OF THE 
SOLUTIONS IN THE THIRD APPROXIMATION 


As we have seen in papers II (§ 7) and IV, the solution to the problem of line formation 
under certain ‘“‘standard”’ circumstances and in the nth approximation depends on the 
roots ka, (a = 1,....,#), of the characteristic equation (II, eq. [91]) and on the solutions 
of the linear equations 


which determine the constants of integration Za, (a = 1, .,n). 


It is perhaps of some interest to note that the solution of the foregoing system of equa- 
tions can be written down explicitly. Thus, in terms of the matrix 


gall =(«:-*), (3) 


and the vectors 


the system of equations which equation (1) represents is equivalent to 


GL=M, (4) 


of which the formal solution is 


L=G'M. (5) 


That we can express the solution for Z in this form in terms of the inverse of the matrix 
G is, of course, obvious. But what is noteworthy in our present context is that for the 


1 Ap. J., 100, 76, 355, 1944. These papers will be referred to as “II” and “IV,” respectively. 
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particular matrix (2) it is possible to write the inverse quite simply. For, considering 
the determinant | G| of G, it is evident that it vanishes when any two of the y’s are set 
equal; similarly, it also vanishes when any two of the &’s are set equal. Accordingly, | G| 


cannot differ from 
1 
[] [] (6) 


lyn a¥p 


[] 


i, 2 


by more than a constant factor if we include in the two products in the numerator all 
distinct combinations (4; — mu;) and (ka — kg). (There are thus n(m — 1)/2 terms in each 
of the two products in the numerator, while there are 2” terms in the product in the de- 
nominator.) A simple inspection now shows that we have, in fact, 


if the cyclical order among the y;’s and the &,’s is maintained.” And, since the matrix of 
every minor of G is also of the same form as G itself but with fewer rows and columns, it 
is apparent that the inverse of G can explicitly be written down. Thus, for the case n = 3, 
we clearly have 


2 For x > 3 the cyclical order has to be understood in the sense that if 1,.... , and # are regarded as 
the verticies of an inscribed polygon we enumerate the various sides (7, i + 1) and the diagonals (7, i + 2), 
(i, i + 3), etc., in their cyclical orders. Thus, for » = 6 the required order is maintained in the sequence 
(1, 2), (2, 3), (3, 4), (4, 5), (5, 6), (6, 1), (1, 3), (2, 4), (3, 5), (4, 6), (5,1), (6, 2), (1, 4), (2, 5), and (3, 6). 
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For ” > 3 we have similar representations. 
With G~ as given by equation (8), the constants Z;, Lz, and L; of the third approxi- 
mation can be found according to equation (5) and are 


(1—piki) (1—p3h1) Avo 


(1—peke) (1—ypsk2) 
(ko—ks) (ki— ke) 
(1—ypyks) (1—poks) (1—y3ks) 
(Rs—k,) (ko—ks) 


Substituting the foregoing expressions for the L’s in the equation for the emergent flux 
(LV, eq. [9]), we find that 


3h kiko 
ky koks 


/ 


ky 
(wit 


F, (0) =4b,[1- 


— + + + +2 (10) 


kiko yt]. 
kikoks 


The roots hj, ke, and k; for various values of \ have already been tabulated in IV (Table 
1). With these values for the k’s, the constants Z;, Z2, and L3 were determined for various 
values of and x (IV, eq. [18]). The results of the calculations are given in Table 1. In 
Table 2 the computed residual intensities are given. 

It may finally be noted that in the second approximation the results analogous to 
equations (9) and (10) are 


(1 bok 1) [1 ke | 
1 k 1 k — 
I, = 2) [ 1 (ui tae) +> | ’ 
and 
(12) 


Avg ky + ke 

b, Ki 
II. THE FIELD OF THE ULTRAVIOLET RADIATION IN A PLANETARY NEBULA 


The equation of transfer for the ultraviolet radiation beyond the head of the Lyman 
series consistent with Zanstra’s theory of the hydrogen-line emission in planetary nebu- 
lae has been considered in IT, and the solution found in the form (II, eq. [83]) 


-s> 


4 1 +1,...., +9), 


4(1-n) (1-975) 


|| 
i 
j 
f 
| 
q 
Bos 
Te ic 
T 
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where Lia, (a = 1,....,), are 2m constants of integration, where ka, (a = 1,....,m), 
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are the positive roots of the equation, 


7 


=1 


aj 
1 —yik 


and where the rest of the symbols have the same meanings as in paper II. 


The boundary conditions which determine the constants of integrations are (cf. II, 


eqs. [84] and [85]). 


(14) 


325 


at r=7, for i=1,.....m (15) 
and 
I-;=0 at r=0 for i=1,....,”. (16) 
TABLE 2 
VALUES OF THE RESIDUAL INTENSITY GIVEN BY THE THIRD APPROXIMATION 
x 

0.1 | 0.125] 0.15 | 1/6 | 0.175 | 0.20 | 0.225 | 0.25 | 0.275} 0.3 | 0.325] 1/3 | 0.35 

0.1722| 0.1876] 0.2012| 0.2093} 0.2132) 0.2238| 0.2333! 0.2419] 0.2496] 0.2567] 0.2631) 0.2651] 0.2690 

0.2544} 0.2724] 0.2882] 0.2977] 0.3021) 0.3145) 0.3255] 0.3355] 0.3445] 0.3526] 0.3601] 0.3625] 0.3670 

0.3794) 0.3983) 0.4148) 0.4248) 0.4294) 0.4424) 0.4540) 0.4645) 0.4739} 0.4825) 0.4904) 0.4928) 0.4976 

0.4815} 0.4995} 0.5152} 0.5246) 0.5290) 0.5413) 0.5523) 0.5622) 0.5712] 0.5794} 0.5868) 0.5891) 0.5936 

0.5717| 0.5879] 0.6021] 0.6106) 0.6146] 0.6257) 0.6356] 0.6446] 0.6526] 0.6600! 0.6667| 0.6688] 0.6729 

0.6542] 0.6683] 0.6305) 0.6879| 0.6914] 0.7010} 0.7096] 0.7174] 0.7244] 0.7307] 0.7366] 0.7384] 0.7419 

0.7313} 0.7429] 0.7531) 0.7592] 0.7621/ 0.7701) 0.7772] 0.7836) 0.7894] 0.7947] 0.7996] 0.8011] 0.8040 

0.8042) 0.8133] 0.8213] 0.8260} 0.8283] 0.8345) 0.8401] 0.8451] 0.8496] 0.8538] 0.8575] 0.8587| 0.8610 

0.8739) 0.8803) 0.8860) 0.8893; 0.8909} 0.8953) 0.8993) 0.9028) 0.9060) 0.9090) 0.9116) 0.9125) 0.9141 

0.9409) 0.9446) 0.9479) 0.9498) 0.9507) 0.9533) 0.9556) 0.9576) 0.9595) 0.9611) 0.9627) 0.9632) 0.9641 

0.9736) 0.9759} 0.9779 0.9791} 0.9797; 0.9813) 0.9828) 0.9841) 0.9852) 0.9863) 0.9872} 0.9876) 0.9881 

1.0057| 1.0066} 1.0075 10032 1.0088) 1.0094} 1.0100] 1.0104] 1.0109] 1.0112! 1.0114] 1.0116 


With the solution for the 7;’s given by equation (13), these conditions require that 


— tar: ) = p G=1. n) (17) 
=] 
( 
and 
(¢e1,..... 2). (18) 
(1-975) 
j 


The field of the diffuse ultraviolet radiation in the planetary nebula is best described 
in terms of its net flux, proportional to F, and its density, proportional to J. In our pres- 
ent scheme of approximation 


(19) 


and 
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TABLE 3 
THE FIELD OF THE DIFFUSE ULTRAVIOLET RADIATION IN A 
PLANETARY NEBULA (p=0.5 AND 71=1) 
F/S J/S F/S J/s 
TABLE 4 
THE FIELD OF THE DIFFUSE ULTRAVIOLET RADIATION IN A in 
PLANETARY NEBULA (p=0.5 AND 71=2) Ir 
ti 
F/S | ys || F/S J/S to 
0.0677 | 0.0295 || 1.000............ 0.0713 | 0.0905 
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or, substituting for 7; according to equation (13), we find, after some minor reductions, 
that 


pa + L_ge that) — (20) 
p a=l1 Ra 1 aj 


and 


>> + +1 6S (21) 
i=1 


As an illustration of the foregoing solutions, we shall consider the case p = 3 and as- 
sume for 7; the values 1 and 2. For p = 3 the characteristic roots are (cf. IV, Table 1, the 
entry for \ = 3) 


k, = 3.69160 , 2=0.964672, and kz = 1.34962. (22) 


For 7; = 1, respectively 2, it was found that 


L,= —0.0071701; —0.0030984 , 

—0.0072975; = —0,0031239, 

—0.0047223; L,;= —0.0020265 , (23) 


L_-,= —0.00025276; L-,= —0.00000619 , 
L-.= +0.310465 ; +0.118659 , 
—0.014003 ; L-3;= —0.0035584 . 


With these values for the constants, F and J were computed according to equations (20) 
and (21) for various values of + (< 7;). The results of the calculations are summarized 
in Tables 3 and 4, and the solutions for the case r; = 2 are illustrated in Figures 1 and 2. 
In these figures we have further compared our present solutions on the third approxima- 
tion (full-line curves) with those given earlier by Chandrasekhar’ on the Milne-Edding- 
ton type of approximation (dashed-line curves). 


3 Zs.f. Ap., 9, 266, 1935. 
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ABSTRACT 


In this paper the problem of the radiative equilibrium of a stellar atmosphere in local thermody- 
namic equilibrium is reconsidered with a view to examining in detail the effects of a continuous absorp- 
tion coefficient, xv, dependent on frequency. It is shown that when the material departs from gray- 
ness the corrections which have to be made to the temperature distribution are of two kinds: first, 
because the integrated Planck intensity B is not, in general, equal to the average integrated intensity 
J (= §§Iidwdv/4r) and, second, because the energy density of the radiation does not have the gray 
atmospheric value. To evaluate the corrections arising from these two sources, a systematic method of 
approximation has been developed in which the first approximation is assumed to be given by the solu- 
tion for a gray atmosphere with a «x, equal to a certain appropriately chosen mean absorption coefficient x. 
It is found that the best way of choosing x is to define it as a straight mean of x», weighted according to 


the net monochromatic flux F (1) (r) of radiation of frequency » in a gray atmosphere: 


k= dy (F = integrated flux) . 


Moreover, in a second approximation, 


where J(2) is the solution appropriate for J in this approximation. The solution for J ‘2) has, in turn, been 


obtained in an nth approximation. 
To facilitate the use of the solutions in the higher approximations obtained in this paper, the mono- 


chromatic fluxes F‘) (7) and their derivatives dF\) /dr have been evaluated for a range of values of r 
and a = hv/kT, (T, denoting the effective temperature). 


1. Introduction.—The solution to the problem of radiative transfer in an atmosphere 
in local thermodynamic equilibrium is fundamental to all investigations which are in 
any way related to the continuous spectrum of the stars. More particularly, the basic 
problem is one of solving the equation of transfer 

dl, 


under conditions of a constant net flux of radiation 
foe) fo) +1 
ff Luddy, (1) 
0 0 


in all the frequencies and where B,(7) is the Planck function for the temperature 7 at «. 
The complete solution to this problem depends only on the distribution of temperature 
in the atmosphere; for, once this is known, the source function B,(7) for the radiation of 
frequency v is known at all points in the atmosphere and the determination of the in- 
tensity 7, at any point and in any given direction is immediate (see eq. [25] below). 
However, when the continuous absorption coefficient x, is allowed to be an arbitrary 
function of v and x, the determination of the temperature distribution in the atmosphere 
is by no means a simple matter. Indeed, there exists at the present time no satisfactory 
attempt to solve the basic problem with any degree of generality. And what is already 
known in this connection! can be summarized quite simply. 


1 Cf., for what follows, A. Unséld, Physik der Sternatmos pharen, pp. 113-116, Berlin: Springer, 1938. 
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’ would appear to provide a more satisfactory way of defining x, as it has the further advantages of being 
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From the equation of transfer (1) we readily obtain 


(2) 
where, as usual, we have written 
+1 
(3) 
Defining the mean absorption coefficient k by means of the relation 
1 dK, 
pdx dv= f 
we can re-write equation (2) as an equation for integrated radiation in the form 
dK , 
ar ring (5) 


where 7 denotes the optical depth in the mean absorption coefficient x, measured from 
the boundary of the atmosphere normally inward. Remembering the constancy of the 
net integrated flux, we can integrate equation (5) to give 


K =}Fr-+constant . (6) 


Using the conventional Milne-Eddington type of approximations, we readily derive 
from equation (6) the formula 


fo) +1 


But we cannot, in general, pass directly from this equation for J to an equation for the 
integrated Planck intensity B (= o7*/m). For, according to the equation of transfer 
(1), the constancy of the net integrated flux implies only that 


and we cannot infer (except when the material is “gray’’) that J = B. In other words, 
the formula for the temperature distribution which is in common use, namely, _ 


= 37!(1+4r). (9) 


where 7 is the optical depth measured in the “Rosseland mean” absorption coefficient, 


2 In practice this is defined by 
Jo mar’ 


But this is not the same ag the mean absorption coefficient defined by equation (4). It would therefore 
seem that the use of the conventional Rosseland mean absorption coefficient in the theory of stellar at- 
mospheres is of questionable value. On the other hand, an alternative form of equation (4), namely, 


10 
dp , 


a straight mean. Further remarks relating to this matter are made in § 5. 
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cannot, strictly speaking, be regarded as the solution of the fundamental equations in 
any well-defined scheme of approximation. Thus, while the use of equation (9) cannot be 
strictly defended, there does seem to be some evidence for the belief that the tempera- 
ture distribution derived on the gray-body assumption provides a “first approxima- 
tion’’ to the true distribution in some sense or another. For example, it is well known 
from the studies of Milne and Lindblad that the law of darkening in the different wave 
lengths over the solar disk and the intensity distribution in the continuous spectrum of 
the sun agree quite well with what can be predicted on the gray-body assumption. But 
attempts to improve on this agreement by allowing for a variation of x, with v ona 
“second approximation” generally indicate that the departures from grayness of the 
material is quite large. It would, accordingly, appear that the solution for the tempera- 
ture distribution on the gray-body assumption does provide a good first approxima- 
tion even for substantial variations of x, with v. But, as to in what sense it is a first ap- 
proximation can be understood only when we develop a systematic well-defined scheme 
of approximation and explicitly work out a second approximation. It is the object of this 
paper to provide such a scheme and to estimate the errors which are involved in the 
use of the solution (9). 

2. An iteration method for solving the equation of transfer for an atmosphere in local 
thermodynamic equilibrium and with a continuous absorption coefficient depending on the 
wave length.—According to our remarks in the preceding section, it would appear that the 
temperature distribution in the atmosphere derived on the assumption of grayness of 
the material provides a first approximation. Assuming that this is the case, let x be a 
certain mean absorption coefficient (undefined for the present) such that with x, = x = 
constant we obtain a first approximation. Let 7 denote the optical depth measured in 
terms of this absorption coefficient. Further, let ° 


(10) 
K 


so that 6, is a measure of the departure from grayness of the material. With these 
definitions the equation of transfer (1) takes the form 


(11) 

drt 

We suppose that this equation can be solved in two steps. First, we find the solution of 
the equation 


(1) 


appropriate to the problem on hand, and use this solution in the term which occurs as 
the factor of 6, in equation (11). Thus, the second approximation will be given by the 
solution of 


(2) (1) 


Formally, there is, of course, no difficulty in extending this method of iteration to ob- 
tain solutions to as high an approximation as may be needed. Thus, in the mth approxi- 
mation we shall have 


(n) 
dl, n). n— 
pS = 1 — Be + Br". (14) 


_ 3 Cf. Unsdld, op. cit., pp. 107-109 (see particularly Tables 24 and 26). 
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where the m constants L{”, (a = 1,...., 2 — 1), and Q are determined by the 
equations (II, eq. [21]) 
a=l —pika 
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This method of iteration will obviously converge if 6, is sufficiently small. But the point 
to which attention may be drawn at this stage is that in practice the success of this 
method is not impaired even when 6, takes moderately large values of the order of 
2 or 3 (see § 5 below). 

In this paper we shall consider only the first two approximations. 

In their integrated forms, equations (12) and (13) are 


(1) 
dr 
and 
dI() 
I B tu f by dv. (16) 


The foregoing equations have, of course, to be solved under the conditions of a constant 
net integrated flux. In the first two approximations this condition requires, respectively, 
that 


BO) = JQ) (17) 
and 


1 p> ct gy 
(2) — 
BO) =J +5f ududy. (18) 


Equation (18) can be re-written in the form 
(1) 


In other words, in the second approximation the integrated Planck intensity B differs 
from J by an amount which depends on the nonconstancy of the monochromatic fluxes 
F,. From general considerations we may expect (and this is confirmed by the calculations 
to be presented in § 5) that, unless 6, varies too widely over the relevant ranges of », 
the departures from constancy of the monochromatic fluxes F, will be of the second order 
of smallness. It is precisely on this account that the temperature distribution derived 
on the assumption of grayness of the material is as satisfactory as it has been found to be 
in practice. 

We now proceed to the solutions of equations (15) and (16). 

3. The solution in the (1, n) approximation.—Equations (15) and (17) together reduce 
to an equation of transfer of a standard type which has been treated in sufficient detail in 
earlier papers of this series. Adopting in particular the scheme of approximation de- 
veloped in paper II we can, in the nth approximation, write (II, eqs. [19] and [26]) 


(1) 
1) =3F > La + +7} G=+1,.... 2). (20) 


*Ap. J., 99, 180, 1944, and 100, 76, 1944. These papers will be referred to as “I” and “II,” respec- 
tively. 
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and the k,’s are the m — 1 distinct nonzero positive roots of the characteristic equation 
(II, eq. [10}) 


2= (22)8 
this the for J‘ is (IT, eq. [29]) 
J) =3F (r+q[r]) =3F +O 4 (23) 
The corresponding formula giving the temperature distribution is 
(T))4= 378 (24) 
a=l 


Corresponding to the distribution of temperature (24), there is a determinate dis- 
tribution of the intensities at various frequencies and at various levels. Thus, the in- 
tensity of the radiation of frequency v at an optical depth 7 and in a direction making an 
angle 3 with the positive normal is given by 


8) =f e- (t=) ee 9B, (T) sec ddt (0< < 


(25) 
0 


In terms of this solution for /{! (7, 8) we can readily derive explicit formulae for quan- 
tities such as J,, F,, etc. In fact, it can be shown quite generally that 


+1 , 
4 


+ B, (T\?) (7 dt. 


where £;,; (x) stands for the exponential integral 
@ 
Ein = fo dw. (27) 
In our further work we shall find that we also need expressions for the quantities 
= ‘ 28 
(28) 
From the equation of transfer (15) we find that 


+41 g 7) +1 2 
—1 dr Jj 


where 
€j) odd — 1 if j is odd 30 
= Q otherwise. (30) 


5 In the stfmmation on the right-hand side of this equation there is no term with i = 0. 
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Using equation (26), we can re-write equation (29) as 


+1 (1) fo) 
f dl, vidu= f B, (T) E; (t— dt 


1 dt t 
0 J 


(31) 


Equations (24), (25), (26), and (31) may be said to represent the solution to our prob- 
lem in the (1, 2) approximation referring to the fact that the equations of the first ap- 
proximation have in turn been solved in an mth approximation. 

4. The solution of the equations in the (2, n) approximation—Turing next to the 
solution of the equations (16) and (18), we have to consider the equation 


(1) 


-3f 1 dutuf w—5f dudv. (32) 


In solving this integrodifferential equation we shall follow the methods developed in the 
earlier papers of this series and replace the integrals which occur on the right-hand side 
of this equation by sums according to Gauss’s formula for numerical quadratures. 
Thus, in the mth approximation, the equivalent systems of linear equations are 


pi tui f 5, dy 
(33) 
1° 


where we have used /{?) and J‘. to denote, respectively, (7, and J{(r, and 
the rest of the symbols have the same meanings as in paper II. 

The system of equations represented by (33) is most conveniently solved by the 
method of the variation of the parameters. Thus, since the homogeneous part of the 
system of equations (33) is of the same form as that considered in II, equation (6), we 
seek a solution of our present system of the form (cf. II, eqs. [18] and [26]) 


a=l1 


(34) 
etkar (2) (2) 
where, as we have indicated, L®’, L®, (a = 1,...., — 1), and Q® are all to be con- 


sidered as functions of 7. It will be noticed that in writing the solution in this form we 
have treated as variable only (2m — 1) of the 2m constants of integration which the 
general solution of the homogeneous system associated with equation (33) involves. 
This is, however, permissible in view of the fact that equation (33) admits the flux 
integral 


I” = constant , (35) 


and we can arrange that this constant of integration has the same value as in the (1, 2) 
approximation corresponding to the circumstance of a given constant net integrated 
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Substituting for /{?) from equation (34) in equation (33), we obtain the variational 
equation 


+kar d d (2) 
1+ ike dr font 1—ypika dr drt 
(1) 


(1) 


Of the 2n equations melee in the foregoing equation, only 2 — 1 are linearly 
independent, since the equation derived from (36) by multiplying with a; and summing 
over all i’s is identically satisfied (cf. II, eq. [25]). The rank of the systems (36) is, ac- 
cordingly, (2n — yl agreement with the fact that we have only (2n — 1) functions 
L®), (a = 1,. , n — 1), and Q® to determine. 

The order of the system of equations (36) can be further reduced to (2x — 2) bya 
proper averaging of the continuous absorption coefficient x, to yield a x to be used in 
the first gray-body approximation. Thus, multiplying equation (36) by a;u; and sum- 
ming over all 7’s, we obtain 


(36) 


(2) (1) 
dQ 2 dl, i 

1 = ’ 

use having been made of the relations (II, eq. [24] and eq. [58] below) 
2 
> Gigi = and 0. (38) 


Accordingly, if we arrange that 
(1) Kk 


: v = 0 ’ 39 
y=2f 6 dv (39) 
we shall have the integral 
= constant = +AQ (say). (40) 
But equation (39) implies that the mean absorption coefficient x has to be defined accord- 
ing to 
dK,” 
cf dv= f dy, (41) 
or, alternatively (cf. eqs. [2] and [5}), 
(42) 


for with this choice of x the departures of x, from x, when similarly averaged, will be 
zero: 


(1) 
dy = if dy =0 (43) 
0 


Attention may be drawn here to a consequence of the foregoing method of averaging 
x, for the important special case when 4, is independent of 7, i.e., when x,/x is inde- 
pendent of depth. In this case we can re-write equation (19) in the form 
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and equation (43) now implies that 
Bt?) = J(?) (x,/k independent of 7). (45) 


With x defined as in equation (42), the variational equations become 


dr p> drt | 
wif 6 dv jp dy i=+1, tn). | 


Though 2n equations are represented in the foregoing equation, only 2m — 2 of these are 
linearly independent, corresponding to the (2m — 2) functions L‘?) and L®), (a = 1, 
.. — 1), which are to be determined. 

In view of that fact that the rank of the system (46) is less than the number of equa- 
tions, it appears that the most symmetrical way of treating the variational equations 
is the following: 

Multiply equation (46) by (m = 1,...., 2m), and sum over all We 
obtain 


n—1 (2) n—1 (2) 


(47) 
m 
»where we have written 
and 
(1) 
m 


and where €m,oaa has the same meaning as in equation (30). It may be noted that in 
deriving equation (47) use has been made of the relation 


2 


There is a simple recursion formula which Dp, satisfies and which enables a direct 
evaluation of this quantity. We have 


or, using equation (50), 
(= ott (m= ,4n), (52) 
ka \m 


which is the required recursion formula. For odd, respectively even, values of m, the 
formula takes the forms 


1 2 
Ds ke D2 j-2, :) (53) 
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and 
1 
a= Dej-1, «. (54) 
Combining the relations (53) and (54), we have 
1 2 1 
Doj-1, Gite D2 5-3, J (55) 
and 
1 2 
On the other hand, in our present notation the equation for the characteristic roots ka 
can be written as (cf. eq. [22]) 
From the recursion formula (52) we now conclude that 
1 1 
and 
2 2 
Ds, a= and Ds, .= (59) 


The formula (55) will now enable us to determine D,,, . successively for all odd values 
of m greater than 3; similarly, the expressions for the even values of m follow from equa- 
tion (56). In this manner we find that 


2 2 2 , 
Dej-1, a= (j=2,....,%) (60) 
and 
2 2 
D2j, a= (j=2,.... m). (61)$ 


—1) (2G —3)R 


6 In terms of these Dz;,a’s it is possible to write down an alternative form of the characteristic equa- 
tion which, in contrast to equation (22), does not require an explicit knowledge of the Gaussian weights 
and divisions. For, if p2;, (j = 0,...., ), are the coefficients of the polynomial P2,(u), so that 


Pan(u) = print , 
7=0 
it is evident that 
n a; n 2; 
7=0 i 1 + wika 7=0 
since the y;’s are, by definition, the zeros of P2,(u). Hence, with the definitions 


2 2 2 


and 
Do = 2 and D.=0, 


the equation for the characteristic roots takes the form 


p2jD2; =(Q, 


j=0 
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Returning to equation (47), we first observe that, since D;, D2, and 6; are all zero 
(cf. eqs. [39] and [49]), then, of the 27 equations represented by equation (47), those for 
m = 1 and 2 are on ‘satisfied. Accordingly, we need consider equation (47) only 
for values of m = 3,. , 2n. For odd, respectively even, values of m, equation (47) 
takes the forms 


Dej-1, a (enter — ether d = (j = 2, n) (62) 


dr T 

and 
(2) (2) 

a=1 


or, substituting for Do;-1,. and D;;,, according to equations (60) and (61), we have 


3 1 1 1 1 


8) 


and 


33 1 1 1 


where, for the sake of brevity, we have written 


pS (« —kar etkar 
dr ’ 


drt 
(66) 


= 


The linear systems represented by equations (64) and (65) can be brought to forms more 
convenient for their solutions by the following procedure. 
Considering, for example, equation (64), we find that the system of equations which 
this represents is 
3 n—l 1 


a=l 
33/1 1 - 
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The form of these equations suggests that we construct the following simpler system of 
equations: 


| = 
a=] 
4 
2 Xa= Us, 
a=l “a 
1 1 3 
= 
= U,-| 


where the U’;’s are defined as in equation (68). 
Similarly, starting from equation (65), we can construct the simpler system 


1 
where the quantities V; (7 = 1,....,— 1) are defined by the recurrence relation 
3 
V j-1 = 60; — 3 Vj-2 — 2 Vi-3— Vi (j = 3, (71) 
and 
Vi= 04. (72) 


We shall now show how it is possible to write down the solution of the equations (69) 


and (70) quite generally. 
Considering first equation (69), we find that the system of equations which this repre- 


sents can be written down as a single vector equation 
KX = 2U , (73) 


where X and U stand for the vectors 
x= (X,, and U= Un-1) (74) 
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and K for the matrix: 


1 ke 
1 
(75) 


Equation (73) expresses a relation between U and a linear transform of X. This rela- 
tion can be inverted into a relation between X and the linear transform of U obtained 
with the inverse of the matrix K.7 Thus 


X=2K-U, (76) 


where K~! denotes the inverse of K. For the particular matrix K the inverse can be 
explicitly written down. But before we do this, it is first necessary to introduce certain 
definitions. 

Consider the expansion of the product 


in descending powers of x. We have 


| | (x? k?) ( 2 ) ( k2 k2 ) 


i#r 
1) n—1 (78) 
aa: 2 p2 n—2 2. 


or, introducing the (w — 1) independent symmetric functions in the (m — 2) variables 
1), 


5, ,=1, 
ly 
1 
ly 
My (79) 
1yn—l1 
My 


’ The existence of the inverse can be inferred from the linear independence of the equations (67), or 
more directly from the nonvanishing of the determinant of K. The determinant of K can be written 
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we can re-write the expansion on the right-hand side of equation (78) in the form 


lyn—1 n—2 
I] (x? — k*) = .. (80) 
ifr 


In terms of the functions S),,,(A=0,1,....,”2—2,andr=1,....,m— 1), defined 
in this manner, the inverse of the matrix K can be written down. We have 


( sip Si, Sh, 1 = 1 ) 
1#1 
[] (2.,-4) (k?_, — (k2_, — k?) 


That the foregoing matrix does, in fact, define the inverse of the matrix K can be veri- 


fied directly. Thus, 
0 
(K-'K) ly m = k 
h=0 I] (k? k?) 
(82) 
Rm I] k?) A=0 
or, according to equation (80) 
2n—3 
m — k?) 
7 


down as it is related to the well-known Vandermonde determinant (cf. O. Perron, Algebra, Vol. 1, § 22, 
pp. 92-94, Leipzig: Gruyter, 1932). We have 


K| 1 TL 1 1 75 
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If / ~ m, the product which occurs in the numerator of the right-hand side of equation 
(83) includes the term j7 = m and consequently vanishes. On the other hand, if / = m, 
the right-hand side of equation (83) clearly reduces to 1. Hence, 


(K-'K) ;, ,,=1 if 1=m and zero otherwise , (84) 
proving the inverse relationship between the matrices (75) and (81). 


Hence, with K~ as defined in equation (81), equation (76) determines X. In terms of 
the components X, of X this equation becomes 


(2) (2) 2n—3 n—2 
Xa=F otter 8) + = 
I] (k2 — (85) 
a 


(a=1,....,#—1). 
Similarly, the solution of equation (70) is found to be 


(2) (2) 2n—2 n—2 
—a 2Ra 
II (k2 — k2) (86) 
(a=1,....,m—1). 
From equations (85) and (86) we obtain 
(2) 2n—3 n—2 
Sr, a [Un41 (7) — RaVn4i (7) | (87) 
I] (Rk? k?) h=0 
and 
Sa, (7) + kaVn41 (7) J (88) 
I] k®) h=0 


In their integrated forms the foregoing equations are 


= > etkar | (7) — RaVaai(r) dr 
0 
and 
n—2 
FL® (1) = e~*at | (7) + (7) | dr 
[] 2-2) (90) 
; (a=1,....,"—1), 
where in equation (89) AL,, (a = 1,....,— 1), are (m — 1) constants of integra- 


tion. However, it will be noticed that in integrating equation (88) in the form (90) we 
have made a particular choice of the constants of integrations so as to bé compatible 
with the requirement that none of the quantities tend to infinity exponentially as 
t— (cf. eq. [34)]). 
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Finally, the constants of integration AZ., (a = 1,....,m— 1), and AQ (cf. eq. 
[40]) which occur in the solution of the variational equations are to be determined from 
the boundary conditions 


I_-;=0 at for i=1,....,m. (91) 


Explicitly, these conditions reduce to (cf. II, eq. [21]) 


(1) 7 (2) 
or, since the L{”’s and Q™ satisfy the relation 
(1) = i= 
Mi GG=1,....,m), (93) 
the equations which determine AL, and AQ are 
n—-1 7 (2) 


This completes the formal] solution to the problem in the (2, 2) approximation. 
5. The solution in the (2,1) approximation.—The discussion of this case is of par- 
ticular interest, as in this approximation 


J = J) ‘ (95) 


and equation (19) becomes 


dF,” 
B® =J +z de. (96) 


Substituting for J“ its solution in the first approximation (II, eqs. [30] and [35]), we 
have 


(1) 
B dv). (97) 


It may be recalled that in the foregoing equations the optical depth 7 is evaluated 
in terms of the mean absorption coefficient x defined by. (cf. eq. [42]) 


(1) 
c= f (7) a. (98) 


From equations (97) and (98) it is apparent that, when 6, is independent of depth, 


BO) = BO = 3F é (6, independent of r). (99) 


8 It is possible to write down explicitly the inverse of the linear transformation which appears on the 
left-hand side of this system (cf. Ap. J., 101, 320, 1945, eq. [8]). But, as the practical use to which these 
solutions can be put is limited, we shall not develop further these formalities here. 
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In other words, we have shown that im a Milne-Eddington type of approximation the 
temperature distribution obtained on the gray-body assumption continues to be valid in a 
second approximation if the mean absorption coefficient x (in terms of which + has to be 
defined) is a straight average of x,, weighted according to the net monochromatic flux of 
radiation of frequency v in a gray atmosphere, and if, further, x,/k is independent of depth. 
This answers the question that we asked at the outset (§ 1), namely, as to the conditions 
under which an equation for the temperature distribution of the form (9) may be ex- 
pected to hold even when the material departs from perfect grayness. We may draw 
particular attention to the fact that our answer has not provided any justification for 
averaging x, according to the conventional method of Rosseland. Indeed, the Rosseland 
mean, as customarily defined by 


1 dB, 
dv= f — ay, (100) 


weights x, falsely and overemphasizes the violet region of the spectrum. On the other 
hand, our present discussion would rather suggest that in all preliminary considerations 
it would be preferable (and simpler) to average x, directly, weighting with the Planck 
function corresponding to the effective temperature of the star. One further remark con- 
cerning the requirement of the constancy of x,/x with depth may be made. With the 
negative ion of hydrogen firmly established as the principal source of absorption in 
stellar atmospheres of the solar and neighboring types, the requirement of the constancy 
of x,/k with depth is likely to be met fairly satisfactorily. 

Returning to equations (97) and (98), we see that a practical use of these equations 
will require a knowledge of the monochromatic fluxes F{") at various depths in a gray 
atmosphere. 

Now, according to equation (26), we can write 


0 


where 
Writing 
hy 
(103) 


we can express equation (101) in the form 


30a 


F 232754 1 fa(r) (104) 
where 
o E(t-—1r)dt 
and 4 
1/4 
1.23275.... (106) 


Using for g(7) the solution on the fourth approximation given in II, equation (54), 
Miss Frances Herman and the writer have numerically evaluated the function f, (7) de- 
fined as in equation (105) for various values of a and +r. The results of the numerical 
integrations are summarized in Tables 1 and 2. In Table 1 we give the values of f(r), 
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and in Table 2 these are converted according to equation (104) to give the monochro- 
matic fluxes directly. And finally, in Table 3, we give values of dF‘! /dr obtained by 
direct numerical differentiation of the fluxes tabulated in Table 2. 


TABLE 1 
THE FUNCTION fa(7) 


a 
T 
0 1 2 3 4 6 8 10 12 
ere’ eee 0.601 0.682 0.803 0.982 1.236 2.10 3.82 7.42 $5.2 
_ .529 .619 .757 .960 1.249 2:23 4.20 8.32 17.3 
.472 . 568 .938 1.257 4.58 9.24 19.4 
cule .424 .524 .680 1.259 2.46 4.94 10.19 
. 384 .647 . 893 1.257 2.56 5.29 41-35 24.1 
. 348 .450 .869 1.251 2.64 5.62 12.11 26.6 
.318 .420 .586 .847 1.244 2°72 5.95 13.08 29.1 
.292 .393 .561 .826 1.236 2.79 6.26 14.06 31.8 
.267 . 368 .535 . 804 1.224 2.85 6.56 15.03 34.6 
.346 13 .784 1.214 2.90 6.84 16.00 37.4 
.229 .326 .493 1.203 2.95 7.12 16.96 40.3 
ity .198 .292 .729 1:77 3.04 7.64 18.85 46.2 
2173 . 264 .423 .696 1.152 3.10 8.11 20.70 S22 
wane . 396 . 666 1.126 3.16 8.54 22.49 58.4 
.221 .639 1.101 3.20 8.94 24.23 64.5 
0.124 0.204 0.351 0.613 1.076 3.23 9.30 25.90 70.7 
TABLE 2 
THE MONOCHROMATIC FLUXES F‘))(7) IN UNITS OF F 
a 
1 2 3 4 6 8 10 12 
tee a ee 0.0864 0.1837 0.2074 0.1772 0.0856 0.0314 0.01013 0.00305 
.0784 .2027 .1791 .0911 .0346 .01135 00345 
@:2 .0719 . 1640 .1981 .1801 .0961 .0376 .01261 .00388 
.0663 .1556 1933 . 1804 .1005 .0406 .01390 00434 
.0615 . 1480 1886 . 1802 . 1044 .01521 .00482 
3 .0571 . 1407 1835 .1793 .1078 .0462 .01652 .00532 
.0532 . 1342 1789 .1783 .1110 .0489 .01785 .00584 
.0498 .1282 1745 .1138 .0515 .01918 00638 
.0466 1699 .1755 .1163 .02050 .00693 
.1174 . 1656 .1740 . 1186 .0563 .02182 .00749 
.0413 .1128 . 1616 .1724 . 1207 .0586 02314 .00807 
.0370 . 1042 . 1540 . 1688 .1240 .0628 .02572 00925 
.0335 .0968 .1470 .1651 .1268 .0667 .02824 .01046 
.0305 .0906 . 1406 .1614 .1289 .0702 .03069 .01169 
.0280 .0850 . 1349 1578 . 1306 .0735 .03306 .01293 
0.0259 0.0803 0.1295 0.1543 0.1319 0.0764 0.03534 0.01417 


In Figure 1 we have illustrated the distribution in frequency of the net flux of radia- 
tion at various depths. It is of particular interest to observe how the redistribution in 
the frequencies takes place as we descend into the atmosphere. 


| 
| 
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With the data given in Tables 2 and 3 we should be able to use the solution (97) to 
determine the temperature distribution under a wide range of practical conditions. 


We shall return to such applications in later papers. 


TABLE 3 
THE DERIVATIVES OF THE MONOCHROMATIC FLUXEs: dF‘) /dr IN UNITS OF F 
a 
T 
1 2 3 4 6 8 10 12 

2 ome —0.086 | —0.11 | —0.044 | +0.024 | +0.0585 | +0.0323 | +0.0120 | +0.0039 
0.1 072 .098 .046 015 0044 
.060 .088 .047 | + .007 .0467 .0128 
052 .080 .048 | + .001 0417 .0130 .00465 
0.4.. .046 074 .0485 | — .005 .0367 0281 0131 .00490 
0.5.. 041 .068 .048 | — .008 0327 .0270 0132 .00510 
036 .063 046 | — .010 0295 .0260 0133 00530 
.033 058 045 | — .013 .0264 .0250 0133 .00544 
.030 054 043 | — .015 .0240 .00558 
0.9.. .027 .050 041 | — .016 0215 .0233 0132 .00573 
024 040 | — .017 0193 .0222 0131 00583 
.020 .040 .036 | — .018 .0202 .0128 00598 
.016 034 | — .0185 0121 0186 0124 .00610 
1.6.. 014 .029 .030 | — .0185 .0095 .0170 .0120 .00617 
O11 .026 .028 | — .018 0075 0156 .0116 .0062 
a ae —0.010 | —0.022 | —0.025 | —0.017 | +0.0060 | +0.0142 | +0.0112 | +0.0062 

.20 

a), 
F 

40 

05 

‘ 

i?) 2 4 6 8 10 12 14 


Fic. 1.—The frequency distribution of the net flux of radiation at various depths in a gray atmos- 
phere. The abscissa measures the frequency (in units of kTe/h), and the ordinate measures the flux (in 
units of the constant net integrated flux F). The curves /, 2, 3, and 4 refer to the depths r=0, 0.5, 1, 


and 2, respectively. 


6. The solutions for the temperature distribution in higher approximations.—Turning 
next to solutions in higher approximations, we have (cf. eq. 34 and II, eqs. [28] and [29]) ° 


a=} 
J? = S* 12) (7) enter + S122 (7) (107) 
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or, substituting for QO, L{?) and L®) according to equations (40), (89), and (90), we have 


= a Sx, . enter [ether (7) dr (108) 
0 


where it may be recalled that the constants of integration AQ and AL,, (a = 1, : 
n — 1), are to be determined according to equation (94). Corresponding to the fore- 
going solution for J, there is the temperature distribution 


(1) 
4 = BO) = JO) 4 (7) (109) 
Tv 


Equations (108) and (109) show that, when the material departs from being gray, the 
corrections which have to be made to the temperature distribution’ are of two kinds: 
first, because B is not, in general, equal to J, and, second, because the energy density of 
the radiation does not have the gray-atmosphere value. It is only under very special 
conditions (which we have already discussed in § 5) that these corrections vanish in a 


second approximation. 
The solution (108) simplifies considerably for the case n = 2. For, in this (2, 2) ap- 


proximation there is only one characteristic root (II, eq. [40]): 


b= 1.97208; (110) 


and equation (108) becomes (cf. eqs. [68] and [72]) 
JO = JO (AQ + AL + (5 — 48, d (hur) | 
0 
(111) 
(55-451 + d (har), 
where it may be noted that (II, eqs. [30] and [42]) 
JO = (7 +0.694025 . (112) 


The equations which determine the constants of integrations AQ and AZ; in the solu- 
tion (111) are (cf. eq. [94]) 
(2) 
AL, __L1 (0) (113) 


and 
(2) 


From these equations we find that 


(1 — (1—- — p2k,) (2) 
and 
(2) 


: 
a 
“—R*) 
a 
ap 
— 
Bek) Mok) 
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Substituting these values for AZ; and AQ in equation (111), we obtain 
= Ju) 


3[ (1 (1 —pek,) — 2] 
4 (1+ yiki) 


0 


(63 ky 64) d (kyr) 
(117) 


From equation (117) we find that for r = 0 (cf. II, eq. [71]), 
V3 V3 ky 


(2) 
(0) 4 2 (1+ (1+ 


(118) 
xf (5; —15,+,5,) d (Air). 


The boundary temperature in the second approximation is therefore given by 


V3 k co (119) 
1 


It is seen that a practical use of the foregoing solutions in the (2, 2) approximation 
will require, in addition to the monochromatic fluxes F‘!, a knowledge also of the quanti- 
ties 


(1) (1) 
(1) +1 dl, 3 (1) +1 dl, 4 
as these will be needed in the evaluation of 5; and 64. However, since the first approxima- 
tion refers only to a gray atmosphere, M‘!) and N{ can be expressed in terms of F{. 
For, from the equation of transfer appropriate to this case (eq. [12]), we can readily 
show that . 


d Ky? 1p), 
dt dt 


= K" (2) (121) 
and 
(r); (122) 


and, accordingly, M{?(r) and N{!(z) can be obtained successively by quadratures 
involving at each stage only known functions. Numerical work relating to these quanti- 
ties with applications to practical problems will be given in later papers. 


In conclusion, I wish to record my indebtedness to Miss Frances Herman for assistance 
with the very laborious numerical work which was involved in the preparation of 
Tables 1, 2, and 3. 
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ON THE RADIATIVE EQUILIBRIUM OF A 
STELLAR ATMOSPHERE. VIII 


S. CHANDRASEKHAR 
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Received February 14, 1945 


ABSTRACT 


In this paper the methods which have been developed in the earlier papers of this series have been 
applied to the solution of the equation of transfer 


dl 
=I1(r, uw) — B=- Idp + 
dr 2 —1 


representing the radiative equilibrium of an atmosphere exposed to a parallel beam of radiation of flux 
aF per unit area normal to itself and incident at an angle 6 normal to the boundary of the atmosphere. 
General solutions in the nth approximation have been found; and it is shown that, on the method of 


solution adopted, it is true in every approximation that 


_ sec 

B(o) vs 
It is further shown that the angular distribution of the “reflected” radiation for the problem on hand is 
related, in a simple way, with the law of darkening in an atmosphere characterized by a constant net 
flux and with no incident radiation. Thus, elementary and simple proofs have been obtained for the 
results first established by E. Hopf, using a special theory involving integral equations. 


1. Introduction.—A problem in the theory of radiative transfer which has considerable 
practical interest is the one relating to the reflection effect in eclipsing binaries. This 
problem was first considered by A. S. Eddington! and later in greater detail by E. A. 
Milne.? As was shown particularly by Milne, the problem basic to the consideration of 
this effect is that of the radiative equilibrium of a semi-infinite atmosphere exposed to a 
parallel beam of radiation of flux +¥ per unit area normal to itself and incident at an 
angle 8 normal to the boundary of the atmosphere (see Fig. 1). In other words, we re- 
quire to solve the equation of transfer (cf. Milne, of. cit.) 

dI ist 

2J_, Idu—jFe (1) 
where J denotes the specific intensity of the integrated radiation derived from the mate- 
rial. The solution of the equation of transfer (1) was obtained by Milne in an approxima- 
tion equivalent to that of Schwarzschild’s and applied to the reflection problem. Later, 
certain exact results relating to the solution of equation (1) were proved by E. Hopf.’ 
But apparently there has been no attempt to obtain solutions explicitly in approxima- 
tions higher than those of Milne’s. In this paper we propose to apply to the solution of 
equation (1) the methods which have been developed in the earlier papers of this series.‘ 
And, as we shall see, the application of these methods to the solution of equation (1) 
leads also to simple and elementary proofs of the exact results of Hopf. 


1M.N., 86, 320, 1926. 

2 M.N., 87, 43, 1926; also Handb. d. Ap., 3, No. 1, 134-141, 1930. 

3 Mathematical Problems of Radiative Equilibrium (Cambridge Mathematical Tract No. 31), pp. 55- 
59, Cambridge, England, 1934. 

4 See particularly Ap. J., 100, 76, 1944. This paper will be referred to as “II.” Familiarity with the 
methods and results of this paper will be assumed in the present investigation. 
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2. The solution of equation (1) in the nth approximation.—As in paper II, we replace 
the integral on the right-hand side of equation (1) by a sum according to Gauss’s formula 
for numerical quadratures and obtain an equivalent system of linear equations which in 
the nth approximation is 

dl; 
Bi 


4+1...... +m). (2) 


where the various symbols have the same meaning as in II. 

It is seen that the homogeneous system associated with equation (2) is the same as 
that considered in II, §§ 2-6. Accordingly, the complementary function for the solution i 
of equation (2) is the same as the general solution (II, eq. [18]) of the homogeneous sys- 


TF 


I 


Fic. 1 


tem. To complete the solution we need only find a particular integral. This can be found 
in the following manner. 
Setting 
in equation (2) (the /,’s are certain constants unspecified for the present) we verify that 


we must have 
h;(1+ sec 8B) = +1 (G=41....., 4m). (4) 


Equation (4) implies that the constants /; must be expressible in the form 


h; (i=+1, +n). (5) 
where the constant vy has to be determined by the condition (cf. eq. [4]) 
a 
= y 
27 (6) 
In other words, 
with 
1 
Y= (8) 


{ 
— 
q 
x 
> 
Ss 
t 
: 
4 
) 
3 4 
= 
q 
q 
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represents a particular integral of the equation (2). Adding to this particular integral 
the general solution of the homogeneous system which is compatible with the bounded- 
ness of the solution for r— ~, we have 


n—1 
M.e 
where Ma, (a = 1,....,2 — 1), and X are constants of integration, to be determined 
by the boundary conditions at 7 = 0, and the k.’s (a = 1,...., m — 1) are the ( — 1) 
distinct nonvanishing positive roots characteristic of the equation (cf. IL, eq. [11]) 

(10) 


At 7 = 0 we have no radiation derived from the material. Accordingly, we should re- 
quire that 

I_-;=0 at 7=0 andfor 7=1,....,m”. (11) 
From equation (9) we now conclude that the equations which determine the M,’s and 
X are 


ft 1 — ike 1 — sec B 


In terms of the solution (9) for the 7,’s we can obtain a convenient formula for J, de- 


fined b 
y 1 
(13) 
-1 


in our present scheme of approximation. We have 


Using equations (6) and (10), we can reduce the foregoing equation to the form 


fr 
The temperature distribution in the atmosphere is given by (cf. eq. [1]) 
B(r) =J 41 (16) 
Substituting for J from equation (15), we obtain 
B(r) =1F tee xp (17) 
a=l 
As + — 0, respectively ~, we have 
a=1 


and 
=4FX. (19) 


= 
— 
a 
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The law of darkening in the atmosphere corresponding to the temperature distribu- 
tion (16) is readily found. We have ~ 


a=] 
M. 


This completes the formal solution. In the following section we shall show how the 
explicit forms of the solutions in the various approximations can be written down di- 
rectly in terms of the respective solutions of the associated homogeneous equation 


dl 1 
(21) 


given in paper IT. 

3. The relation of the solution of equation (1) with that of the associated homogeneous 
equation (21).—In the preceding section we have obtained the solutions for the tempera- 
ture distribution and for the law of darkening for an atmosphere exposed to an incident 
parallel beam of radiation. The explicit numerical forms which these solutions will take 
in any particular approximation will depend on the solution of the linear equations for 
the constants of integration M,, (a = 1,....,” — 1), and X (cf. eq. [12]). However, we 
shall now show how a straightforward solution of these equations can be avoided by ex- 
pressing the constants M, and X in terms of the solutions of the homogeneous equation 
(21) in the respective approximations. The analysis which follows is closely related to 
that of II, § 6. 

Consider the function G (u) defined by 


« 
Cw) = * 


According to the boundary conditions (12), 
G(p,;) =9 (23) 


This fact enables us to determine G () explicitly. For, by multiplying equation (22) by 
the function (cf. II, eq. [58]) 


(1 —p sec B)R(u) = (1 —p sec B) [] (1 ken), (24) 


we obtain a polynomial of degree 7 in u which vanishes for u = wi,i = 1,...., m. The 
function G (u) R (u) (1 — wsec 8) cannot, therefore, differ from the polynomial (cf. II, 
eq. [59]) 


P(u) = [] (25) 
7=1 


by more than a constant factor, which can be determined by comparing the coefficients 
of the highest powers of u in the two expressions. In the former the coefficient of yu” is 
unity, while in the latter it is 


(—1)"kyko.... secB. (26) 
Hence, 


secB P(p) (27) 


G(u) =(—1) "hike... 


1—y sec B 
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Comparing this with the function S(u) introduced in II, equations (56) and (61), we see 
that 
sec B 


G(u) sec B 


S(u). (28) 


Several consequences follow from equations (27) and (28). 
First, we observe that, according to equation (28), 


G(0) = X §(0)sec B. (29) 
But it has been shown in II that (cf. II, eqs. [62] and [70]) 


5(0) = : (30) 
Hence (cf. eq. [22]) 
_ y secB 
G(0) = (31) 
From equations (18), (19), and (31) it therefore follows that 
B(O) _ sec B 


This is an exact result, since it is valid in all approximations. Relation (32) has been 
proved earlier by Hopf (of. cit.) but by entirely different methods of more advanced 
character. 
By turning next to the constants M,, it is apparent from equation (22) that 
limit 


Using the explicit form of G(u) (eq. [27]) in the foregoing equation, we find that 
sec B 


Ma= (—1) "hike... X 1— (sec B/ha) Ra (u=k-') .n—1), (34) 
where we have written 
= (a=1......"-1). (35) 


On the other hand, the constants Z, in Il, equation (19), are given by (cf. II, § 6) 
P(u=ka’') 


Hence, 
Ra 


which relates the constants of integration for our present problem with those of the solu- 
tion of the homogeneous equation (21) to which the Z,’s refer. 


oy 
4 
. 
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Again, from equations (22) and (27) we infer that 


limit 
P (cos B) i 
cos 
= sec B B) 
or, more explicitly, 
[| (1 cos 8) 
a=l 
On the other hand, we know that (cf. eq. [8}) 
1 
y= (40) 
1 — yp? sec? B 


Equations (39) and (40) therefore determine XY. However, the evaluation of X is effected 
more conveniently by expressing the right-hand side of equation (40) differently, as fol- 
lows: 

Consider the function 


n a; 
T(x) =1- 41 
(x) i= (41) 
Multiplying this equation by 


[] (1-23), (42) 


j=1 
we obtain a polynomial of degree » in x, which clearly vanishes for 


Accordingly, 
T(x) T] (44) 
cannot differ from 
x Il (x —k?) (45) 
a=l 


by more than a constant factor. The constant of proportionality can be determined by 
comparing the coefficients of the highest powers. In this manner we find that 


a=l1 
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But (cf. eqs. [40] and [41]) 
(47) 
Hence, 
(1 — yj sec? B) 


n 
I] (sec? B —k?) 


a=1 


Comparing the two expressions for y given by equations (39) and (48), we find that 


[] (cos 


cos 
X= 49 


Il (1+. cos B) 
a=l 


But (cf. II, eqs. [58], [59], and [61]) 


[] (cos 8 


S(—cos B) =kyko.... Ra-i (50) 
a=} 
(1+. cos B) 
a=1 


Hence, combining equations (49) and (50) and remembering that (II, eq. [68]) 
(51) 


we obtain 
X =3S(-—cos B)cosp. (52) 


where it might be recalled that S(u) determines the law of darkening for the standard 
problem considered in II (cf. II, eqs. [32] and [56]). 
Substituting for X from equation (52) in equation (28), we have 


G(u) =3 S(u)S( —cos B). (53) 


Now the law of darkening given by equation (20) can be expressed simply in terms of 
G(u). We have (cf. eqs. [20] and [27]) 


I(0,u) (54) 
Hence, 
1(0, n) = 17 cos 8), (55) 


which relates the angular distribution of the “reflected” radiation for the problem under 
consideration with the law of darkening in an atmosphere characterized by a constant 
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net flux and with no incident radiation. This relation between these two problems ex- 
pressed by equation (55) is an exact one, inasmuch as it is true in all orders of approxi- 
mations. Our present derivation of equation (55) thus provides an alternative elemen- 
tary proof of a result first established by Hopf by means of the theory of integral equa- 
tions. It is remarkable that our present method of solving the equation of transfer should 
retain in every approximation all the known exact features of the problem. 
Applications of the theory presented here to problems of the reflection effect in eclips- 
ing binaries will be found in a later paper. 


s This is not the case in the method of approximation used by Milne. His solution (M.N., 87, 48, 
eq. [14], 1926) is not, for example, strictly compatible with the result (eq. [55]). 
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THE WOLF-RAYET SPECTROSCOPIC BINARIES 
HD 186943, HD 193928, AND HD 211853* 


W. A. HILTNER 
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ABSTRACT 


These three Wolf-Rayet stars have been investigated as spectroscopic binaries. 

For HD 186943, spectral type WN5, hydrogen absorption lines were measured, as well as the emis- 
sion bands He 11 4686, N v 4603, and NV v 4619. The absorption lines belong to the early B-type com- 
panion. The following elements have been derived: P = 9.550 days, e = 0.0 (assumed), Kwr (He 1 
4686) = 235 km/sec, Kp = 65 km/sec, Mwr sin? i = 5.8 ©, and Mx sin? i = 21.0 ©. There are pro- 
nounced differences in the gamma-velocities of the three emission bands, none agreeing with the gamma- 
velocity for the hydrogen absorption. The amplitude of the velocity-curve for V v 4603 and N v 4619 
is smaller than that for He 1 4686. Also, the V v velocity-curve shows a phase shift of 0.8 day with re- 
spect to He 11 4686. 

The spectral type of HD 193928 is about WN6. Only the Wolf-Rayet component of the system is 
observable. The following elements have been derived from 30 McDonald spectrograms: P = 21.64 
days, e = 0.0 (assumed), K = 130 km/sec, y (He 11 4686) = +59 km/sec, y (NV v 4603) = —119 
km/sec, and f(m) = 4.94. Absorption lines near \ 3873 and \ 4453 have been measured and identified 
as Hei displaced to the violet by approximately 1200 km/sec. These lines were abnormally strong on 
July 14, 1944, and from observed displacements and changes in dilution we conclude that the Her 
atoms recede from the star with approximately uniform velocity. A similar conclusion was suggested 
from earlier observations of HD 214419. 

The spectrum of HD 211853, spectral class WN6, has well-defined absorption lines as well as emis- 
sion bands. For the emission bands alone, the following elements are obtained: P = 6.6864 days, e = 
0.0 (assumed), K = 245 km/sec. The gamma-velocities of He 1 4686 and N v 4603 are in good agree- 
ment, both bands being displaced to the red, with respect to V rv 4058, by 200 km/sec. None of the 
emission gamma-velocities agrees with that of the absorption lines. The measurements of the absorption 
lines give a strange velocity-curve. If one interprets the velocity-curve as that of the early B companion 
with K = 95 km/sec, we have the following minimum masses for the two components: /wr sin’ 7 = 
7.6 © and Mg sin* i = 19.6 ©. The line He 11 4686 undergoes variations in structure. 

The available data on discrepant gamma-velocities are collected and briefly discussed. 


In 1941 O. C. Wilson! announced the discovery of three spectroscopic binaries in which 
one component of each system was a Wolf-Rayet star. Up to that time only one other 
known spectroscopic binary possessed a Wolf-Rayet component, HD 193576,? and one 
visual double, HD 219460.* Since this announcement by Wilson, only two other Wolf- 
Rayet binaries have been discovered.‘ In 1940 Wilson’ called attention to a number of 
Wolf-Rayet stars that have undisplaced absorption lines, which he suggested might be- 
long to a second component. Stars HD 186943 and HD 211853 were in this list. To this 
list the author wishes to add HD 193077, which is probably a spetroscopic binary. The 
star HD 168206, on Wilson’s list of Wolf-Rayet stars with absorption lines, has been 
found to be a spectroscopic binary from spectrograms taken on July 25, 1941, 602™ 
U.T., and on June 16, 1944, 4"56™ U.T. The relative displacement of carbon blend near A 
4635 on the two spectrograms is about 250 km/sec. 


* Contributions from the McDonald Observatory, University of Texas, No. 108. 
Pub. A.S.P., 53, 295, 1941. 
2 Wilson, Pub. A.S.P., 51, 55, 1939. 3 Tbid., 52, 404, 1940. 


4HD 214419, Dean B. McLaughlin and W. A. Hiltner, Pub. A.S.P., 53, 328, 1941, and HD 152270, 
O. Struve, Ap. J., 100, 199, 1944. 


5 Pub. A.S.P., 52, 405, 1940. 
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THE SpeEctrRUM OF HD 193928 


Note that on July 13.44, He 1 3889 is normal and that He 1 4471 is present as a weak sharp line. On 
July 14.42, He 1 3889 and He 1 4471 are abnormally strong and sharp. On July 15.45, He 1 3889 has not 
changed in intensity from the previous day, but He 1 4471 is appreciably weaker. This phenomenon is 
apparently not periodic, as is shown by the spectrogram taken on June 22.34 at phase 12.23. 
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No detailed spectroscopic investigation of thesé three stars had been undertaken, and 
the binary nature of only the last one had been confirmed.’ Consequently, the three stars 
were placed on the observing program of the 82-inch reflector of the McDonald Observa- 
tory in the summer of 1944. The Cassegrain spectrograph with two quartz prisms and 
500-mm camera was used. The dispersion is 42 A/mm at \ 4000 and 70 —_ mm at d 4700. 
Eastman 103a-0 emulsion was employed throughout. 


HD 1869437 


The spectrum of this star is WNS. The predominating feature in the photographic re- 
gion of the spectrum is He u 4686. Its width is approximately 56 A. The only other meas- 
urable emission features are N v 4603 and N v 4618. These bands are relatively narrow— 
only 17 A in width. There are weak diffuse emission features of V tv 4058 and He u of the 
Pickering series, but all are too indefinite to measure for radial velocity, The absorption 
features of the spectrum are extremely weak, characterized primarily by hydrogen. The 
line He 1 4471 is observed only infrequently. Wilson suspected absorption features in this 
star and attributed them to the companion star. 

The line He 11 4686 was measured for radial velocity on 32 spectrograms. Settings 
were made only on the center of the band, and no attempt was made to measure the posi- 
tion of the edges in order to obtain an improved value of the displacement of this band. 
There is some indication that the band undergoes periodic variations in shape, but the 
evidence is not conclusive. This point merits careful photometric observations. 

The period of 9.550 days was determined with the help of Wilson’s published relative 
velocities. Spectrograms taken 0.2 day apart do not suggest a period of about 1 day. 

Radial velocities as determined by He 11 4686 are recorded in Table 1 and are plotted 
in Figure 1. The wave lengths employed in computing radial velocities for this star and 
the subsequent stars studied in this paper are recorded in Table 2. Phases were computed 
in days with the formula 


Phase zero = JD 2431262.20 + 9.550E . 


The measurements are well distributed in phase, and a solution with the eccentricity 
near zero is suggested. The following elements were obtained graphically, assuming that 
the eccentricity is zero: 


P = 9.550 days K = 235 km/sec 
e = 0.0 (assumed) y = +115 km/sec 


The only other emission features measured in HD 186943 are N v 4603 and WN v 4619. 
When plotted separately, they give very similar velocity-curves except for a difference 
in the gamma-velocity. The means for the two lines are recorded in Table 1 and are 
plotted in Figure 2. For 31 spectrograms we find that 


Nv 4603 — Nv 4619 = —150 km/sec. 


The velocity-curve for the Wolf-Rayet component of HD 186943 as determined by NV v 
again suggests an eccentricity near zero but an amplitude smaller than that obtained for 
He 11 4686. The difference is large (160 and 235 km/sec) and is certainly real. Another 
striking difference is brought out when the velocity-curves of He 11 and N v are shown to- 
gether (see Fig. 2). A phase difference of about 0.8 day is strongly suggested. A phase 
difference of N 1v and He 11 of the Pickering series was suggested by the radial-velocity 
measurements of HD 214419* but was not so striking as in HD 186943. 


6 McLaughlin and Hiltner, of. cit., p. 330. 
7@ = 19h42™2; § = +28°01’ (1900); sp., Oa; mag. ptm., 9.98. 
8W. A. Hiltner, Ap. J., 99, 278, 1944. 
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Hydrogen absorption lines were measured on 28 spectrograms. At best, the lines were 
difficult to measure. A magnification of only seven times was employed during measure- 
ment. The lines Hy, Hé, H 8, H 9, H 10, and H 11 were measured when possible, but the 
number of absorption lines measured per plate varied from one to six. The mean velocity 
for each plate is recorded in Table 2 and is plotted in Figure 3. The scatter of the measure- 


TABLE 1 
RADIAL VELOCITIES OF HD 186943 
Vetocity (Km/SEc) 
PLATE Date, U.T. PHASE 
He 11 4686 Nv Habs 
1944 June 15.314 3.764 +180 +228 — 44 
18.285 6.735 —152 + 9 + 1 
XS Bee 20.256 8.706 — 18 — 69 + 2 
21.241 0.141 +152 — 26 +133 
21.297 0.197 + 76 + 43 + 13 
22: 242 1.172 +304 +168 — 22 
3349. 27.246 6.146 — 69 + 55 — 49 
So! Seer 27.303 6.203 — 78 + 37 + 37 
July 1.290 0.640 +221 + 90 + 19 
2.313 1.663 +355 + 85 — 75 
4.327 3.677 +296 +230 —144 
5.371 4.721 + 90 +217 + 64 
6.311 5.661 — 14 + 82 + 36 
7.220 6.570 —122 + 75 + 70 
11.244 1.044 +280 +111 — il 
<a 12.226 2.026 +307 +217 — 23 
12.350 2.150 +343 +217 6 
4 13.274 3.074 +356 +277 — 93 
13.329 3.129 +352 +246 — 32 
15.203 5.003 +118 +153 + 19 
15.246 5.046 + 42 +120 — 30 
16.184 5.984 — 26 + 68 + 4 
17.203 7.003 —151 — 38 + 87 
<a ae 17.357 7.157 —103 — 90 + 90 
18.179 7.979 — 98 — 99 +123 
<n 20.285 0.535 +174 + 76 + 7 
Dec: 6.667 —108 — 6 + 18 


ments is large. However, a variation in velocity in a phase opposite to the emission is 


strongly suggested. An approximate solution follows: 
e = 0.0 (assumed) 
K = 65 km/sec 
vy = 10 km/sec 


Assuming that the absorption lines result from the companion star, which we classify 
only as an early B-type star, and accepting the He 11 4686 velocity-curve as representa- 
tive of the Wolf-Rayet star, the following minimum masses are derived: 


My sin*i = 5.80 


Myx sin? = 21.00. 
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HD 193928° 


The spectrum of this star has a type between WN5 and WN6. Its emission features, 
especially He 11 4686, are stronger than in any of the five known Wolf-Rayet spectro- 
scopic binaries of the nitrogen sequence. In the photographic region, in addition to He 
4686, there are emission bands of He 11 of the Pickering series, N m1, N 1v, and N v. Hy- 


km/sec T T T T T T 
+400 
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Phase in days 
Fic. 1.—The velocity-curve of HD 186943 from He 11 4686 
TABLE 2 
STELLAR WAVE LENGTHS FOR THE DETERMINATION 
OF RADIAL VELOCITIES 
Line Wave Length Line Wave Length 


drogen emission must be weak. Only one absorption feature is normally present: a dif- 
fuse line at approximately \ 3873. Details of the absorption in this star will be given in a 
later paragraph. 

In addition to the strong band He 11 4686, N v 4603 and N tv 4058 were measured for 
radial velocity. Other emission bands were not measured, since they are fainter and ‘are 


%a = 20517™8; 5 = +36°36’ (1900); sp., Oa; mag. ptm., 9.43. 
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Fic. 2.—The velocity-curve of HD 186943 from N v 4603 and - v 4619. The curve for He 11 4686 is 
drawn in to show the differences in amplitude and phase. 
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Fic. 3.—The velocity-curve of HD 186943 from hydrogen absorption 


dif 


* 
0 
2 4 6 10 12 14 
F 
\ 
is 
r 
Br 
: 
| 
fy 


SPECTROSCOPIC BINARIES 361 


relatively difficult to measure on narrow spectrograms. Widened spectrograms were diffi- 
cult to obtain, for the magnitude is 10.5 and the period was thought to be about 1 day 
after the first few spectrograms were obtained. The line N 1v 4058 appears to vary ap- 
preciably in structure. Very frequently the band possessed a narrow band superimposed 
on the wide band and appeared near either the violet or the red side of the broad band, or 
both. The narrow band seems to appear more frequently on the violet edge, but the evi- 
dence is not conclusive. The velocity-curve for this band is distorted, and this distortion 


TABLE 3 
RADIAL VELOCITIES OF EMISSION FEATURES IN HD 193928 


VeLocity (Km/SEc) 
PLATE Date, U.T. PHASE 
He 4686) N v 4603 | N 1v 4058 
GK! Seca 1944 May 9.456 |. 10.626 + 5 — 60 — 54 
16.325 6.215 +200 + 49 — 27 
18.367 8.257 +119 + 44 44 
ce ae 21.350 11.240 + 40 —152 — 189 
22.342 12.232 + 55 — 206 —134 
27.361 17.251 — 25 —294 — 187 
July 0.627 +125 — 140 — 57 
6.369 4.619 +242 — 76 + il 
Co”) ee 7.340 5.590 +185 + 14 — 51 
12.285 10.535 + 51 —124 — 46 
13.388 11.638 + 10 — 164 —171 
13.444 11.694 + 6 —164- —220 
{re 14.370 12.620 — 9 —103 —197 
a 14.423 12.673 — 16 — 186 —184 
15.296 13.546 — 75 —181 —115 
3483 15.348 13.598 — 16 — 207 —132 
15.401 13.651 — 99 —155 —165 
3485 15.452 13.702 — 60 — 194 —145 
16.231 14.481 — 38 —194 —221 
re 16.280 14.531 — 10 — 263 —237 
3493 16.332 14.582 — 73 — 237 —221 
17.256 15.506 — 80 —255 —221 
17.411 15.661 — 48 —251 —182 
17.458 | 15.708 — 54 — 268 —201 
18.234 | 16.484 — 67 —256 —175 
CE 18.292 | 16.542 — 69 — 220 — 162 
hh re 20.339 | 18.639 — 72 — 256 —179 
20.434 | 18.684 — 45 — 204 —179 
Dec. 7.072 | 6.92 | 414 


is probably related to the variation in structure of the band. The behavior of N rv 4058, 
as well as of the other bands, certainly deserves careful analysis in the future. A spectro- 
photometric analysis can be more profitably undertaken now than when the spectro- 
grams for this investigation were taken, since the period and, consequently, the phases 
are now known. Photometric spectrograms would require an exposure of, say, 3 hours, in 
contrast to the present exposure time of 1} hours. 

The radial-velocity measurements of the emission features are recorded in Table 3. 
The lines He 11 4686 and N v 4603 gave almost identical velocity-curves except for a 
difference of gamma-velocity of 


He 4686 — N v 4603 = +178 km/sec. 


é 
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The mean velocities of these two lines are plotted in Figure 4. A period of 21.64 days 
was derived with the help of three relative radial velocities published by Wilson. A period 
of about 1 day has been eliminated by spectrograms obtained in one evening over a pe- 
riod of 0.2 day. The observations are not well distributed in phase, and consequently the 
elements are not well defined. However, the following elements are indicated: 


P = 21.64 days ¥ (\ 4686) = +59 km/sec 
e = 0.0 (assumed) ¥y ( 4603) = —119 km/sec 
K = 130 km/sec f(m) = 4.94 © 
This large mass function is similar to that derived for HD 214419, which in this case is 


f(m) = 4.38. 
The strongest absorption feature in HD 193928 is a diffuse line with an approximate 
wave length of 3873 A. This line was measured on 24 spectrograms (all that were suffi- 
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te] 5 10 15 20 25 30 35 
Phase in days 
Fic. 4.—The velocity-curve of HD 193928 from He 11 4686 and N v 4603 


ciently well exposed in this region), and the measured wave lengths are recorded in Table 
4. The measurements, when plotted (see Fig. 5), indicate that this line possesses a dis- 
placement similar to that obtained for the emission features of the spectrum. Infrequent- 
ly a second line was measured at a wave length near \ 4453. The measurements are also 
recorded in Table 4. The behavior of this line with phase is uncertain because of scatter in 
the measurements. On July 13, phase 11.7 days, this line appeared unusually sharp but 
still faint (see Pl. XXXIV), but \ 3873 appeared normal on this date. On July 14 (phase 
12.6 days) both lines were abnormally strong and sharp. On the following day, \ 4453 
had weakened appreciably, but \ 3873 remained nearly as strong as during the previous 
night. On July 16, \ 4453 was no longer visible with certainty, and \ 3873 was only 
slightly stronger than normal. The behavior of these lines does not appear to be related 
with phase, since on June 22, phase 12.2 days, the line \ 4453 was very weak and diffuse 
and \ 3873 had a normal intensity. If \ 4453 is identified as He 1 4471 and \ 3873 as Het 
3889, we find that the resulting radial velocities are almost identical for 


V (3889) — V 4471 = 6km/sec. 


Furthermore, it is observed that for the two spectograms obtained on July 14, when 
the two lines were sharp and consequently could be measured with greater precision, the 
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TABLE 4 
RADIAL VELOCITIES OF Het ABSORPTION IN HD 193928 


Velocity Velocity 
Date, U.T. Phase r Het 3889 r Het 4471 
(Km/Sec) (Km/Sec) 
1944 May 9.456........ 10.626 3873.12 — 1298 4453.52 — 1204 
JUNG 11.240 3872.58 1240 4452.93 1245 
12.232 3872.92 1214 4451.78 1322 
5.590 3874.22 1113 4454.67 1128 
10.535 3873.59 1162 4452.62 1265 
. 11.638 3873.51 1168 4451.47 1342 
11.694 3872.71 1230 4453.51 1206 
I) 12.620 3873.16 1195 4453.69 1194 
12.673 3873.16 1195 4453.69 1194 
13.651 3873.00 1207 4453.02 1238 
13.702 3873.13 1197 4453.79 1187 
16.484 3872. 64 1235 4453.02 1238 
18.684 3871.62 1314 4452.58 — 1268 
-1000-— * 
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Fic. 5.—The radial velocity observations of He 1 3888 and He 1 4471 
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radial velocities agree within 1 km/sec. If the above identifications are correct, then we 
should expect to observe He 1 4026. The behavior of this line is uncertain. On a widened 
spectrogram taken on July 14, when He 1 4471 was at maximum intensity, there is cer- 
tainly a line present at \ 4010; but a narrow spectrogram obtained the same evening 
gives no indication of a line at this wave length, although it shows \ 3873 and A 4453 
quite strongly. However, in spite of this peculiar behavior of He 14026, it appears logical 
to conclude that we are observing He 1 absorption lines resulting from light from the cen- 
tral core of the Wolf-Rayet star passing through that part of the shell which is ejected 
toward the observer. This conclusion is supported by the total widths of the measurable 
bands. The total width of He 11 4686 is 2400 km/sec. The average velocity of the Het 
absorption lines is 1200 km/sec. The variation in intensity must be interpreted in terms 
of abundance, and any variation in dilution that may have occurred while the lines were 
unusually strong can be interpreted by a variation in distance of the material from the 
central core of the Wolf-Rayet star. 

The variation in intensity and the radial velocity measures of \ 3889 and \ 4471 may 
provide some information on the velocity variation of the ejected He 1 atoms with dis- 
tance from the central core of the Wolf-Rayet star. On July 14 both He 1 3889 and Het 
4471 were strong. Their velocities were, respectively, 1195 and 1194 km/sec. The next 
night, July 15, four spectrograms were obtained, but only two were well exposed. Line 
He 1 3889 was of approximately the same intensity as on the previous night; but Het 
4471 was much weaker, indicating a change in dilution during the interval of twenty-four 
hours. The radial velocities of \ 3889 and \ 4471 are —1202 and —1178 km/sec. The 
velocity given by \ 3889 should be given greater weight. It is evident that in an interval 
of twenty-four hours, during which time a marked change in dilution occurred, the veloc- 
ity of the ejected gas remained constant within the errors of the measurements. In this 
connection we should call attention to a similar phenomenon in HD 214419.° This star 
normally shows the Wolf-Rayet component of He 1 4471, 4026, and 3889. However, a 
second component of these lines was observed on November 13, 1943, with a displace- 
ment of —800 km/sec. On the following night, November 14, only the displaced compo- 
nent of A 3889 was present, still with a displacement of —800 km/sec. Here, again, a 
large change in dilution was observed, but the velocity of the He I atoms remained con- 
stant as the material receded from the core of the Wolf-Rayet star. 


HD 211853! 


The emission features in the spectrum of this star are relatively weak. The line He 11 
4686 predominates in intensity in the photographic region. Its width is 2500 km/sec, and 
it shows variation in structure with phase. The most pronounced asymmetry occurs at 
phase 3.3 (the Wolf-Rayet star on the far side of the system), when the band is shaded 
toward the red. There is some suggestion that V tv 4058 undergoes periodic variation in 
structure, but the band is too weak to make the observations conclusive. The other prin- 
cipal emission features are He 11 of the Pickering series, V 111, VN v 4603, and possibly He 
I, especially \ 4471. Plate XX XIII shows the spectrum of this star in the photographic 
region. Its spectral type is near W NO. The lines He 114686, N 1v 4058, and N v 4603 have 
been measured for radial velocity. 

Absorption features in this star are quite evident, as is shown on Plate 000. The lines 
Hy, Hé, H 8, H 9, and H 10 are always present on well-exposed spectrograms. They are 
usually diffuse, indicating a rotational velocity of the order of 230 km/sec. 

The measurements of the radial velocities of the emission bands give very similar 
velocity-curves except for differences in the gamma-velocity and have been combined to 
give the following elements: 


10 q = 22515™0; 6 = —55°30’ (1900); sp., Ob; mag. ptm., 9.0. 
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6.6864 days 
0.0 (assumed) 
245 km/sec 


(4686) 
¥ (4603) 
(4058) 


TABLE 5 


= +50 km/sec 
= +62 km/sec 
= —125 km/sec 


RADIAL VELOCITIES OF HD 211853 


The period was improved by referring to Wilson’s published relative velocities. A period 
of about 1 day has been eliminated by obtaining as many as 5 spectrograms in the same 
evening over an interval of 0.14 day or less. The measured radial velocities are collected 
in Table 5, and the mean velocities for the emission features are plotted in Figure 6. 


Vetocity (Km/SeEc) 


PHASE 
PLATE Dare, U.T. (Days) 
Emission Hei 4471 
CQ 3175. 1944 May 7.465 1.483 +312 
3284. 17.433 2.391 +185 — 87 — 49 
18.437 3.335 + 6 —140 —158 
19.321 4.219 —176 — 89 — 36 
ee 19.441 4.339 — 166 — 56 — 81 
3317 20.380 5.278 —213 
ee 20.433 5.331 — 202 + 16 + 57 
| ere 21.397 6.295 + 13 — 42 — 61 
21.433 6.331 — 8 —107 — 61 
22.410 0.622 +187 —117 —182 
27.412 5.624 —179 + 35 + 78 
Ls 27.450 5.662 — 137 + 35 + 24 
ce 28.448 6.660 + 25 —114 —121 
3357 30.389 1.915 +275 
CY Seer July 1.358 2.884 + 92 — 32 — 44 
ce 1.435 2.961 +121 — 26 — 76 
2.433 3.959 — 28 — —117 
eS eee 3.322 4.848 —212 — 62 — 19 
ks nS 3.358 4.884 — 134 — 27 — 84 
3388. 3.392 4.918 —175 
3.427 4.953 —148 — 24 + 17 
i. ner 3.458 4.984 —198 — 20 — 100 
CT ee 4.408 5.934 —115 + 14 + 50 
| ee 4.437 5.963 — 68 + 12 + 5 
6 ih 5.461 0.301 + 85 —150 —142 
6.414 1.254 +295 — 133 —117 
8.451 3.291 + 59 ww OF 
i re 11.381 6.221 — 63 + 50 + 5 
Ce. 12.407 0.561 +212 —126 — 146 
ee 12.444 0.598 +198 —116 — 190 
19.474 0.942 +292 — 133 —207 
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Hydrogen absorption, consisting usually of Hy, H8, H9, H10, and H 11, was measured 
on all spectrograms. The mean radial velocity for each spectrogram is recorded in Table 
6 and plotted in Figure 7. The form of the velocity-curve is not well defined because of 
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Fic, 6.—The yelocity-curve of HD 211853 from He 11 4686, N v 4603, and N rv 4058 


TABLE 6 
RADIAL VELOCITIES OF Het 3889 ABSORPTION IN HD 211853 

Velocity Velocity 

Plate Phase. (Km/Sec) Plate Phase (Km/Sec) 
3.291 3873.5 — 1160 OO 3872.9 | —1210 
4.848 3872.3 1260 3871.6 1310 
ie 4.884 3872.9 1210 a 5.662 3871.7 1300 
4.984 3871.8 1300 S690 3871.3 1330 
5.278 3872.7 — 1230 6.660 3870.7 — 1380 


the large scatter. However, the hydrogen-absorption velocity appears approximately 
constant at — 65 km/sec from phase 1.6 days to 4.9 days, at which time there is a marked 
increase in velocity, reaching a maximum of +40 km/sec at phase 5.6 days. The velocity 
then decreases rapidly to — 150 km/sec at phase 0.8 day. Although He 1 4471 was meas- 
ured less frequently than the hydrogen lines, it verifies the general form of the velocity- 
curve, especially the rapid decrease in velocity. 
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The interpretation of the absorption velocity-curve is not clear, and it certainly de- 
serves further careful examination. From phase 1.6 days to phase 4.9 days, when the 
velocity of the absorption lines remains approximately constant, the B-type component 
of the system is between the observer and the Wolf-Rayet component. From phases 4.9 
to 1.6 days, when the velocity of the absorption lines is changing rapidly, the B-type 
component is on the far side of the system. This suggests that the observed variation in 
velocity may be related to a rotation effect of the B-type component. This suggestion 
does not appear adequate, primarily because it is necessary to assume that the two stars 
are nearly in contact and that the resulting mass of the B-type component would become 
unreasonably large. If one wishes to interpret this velocity-curve as resulting from the 
second component, one is faced with an eccentricity of about 0.4, in contrast to an ec- 
centricity near zero obtained for the emission Jines. This discrepancy may possibly find 
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Fic. 7.—The radial velocity variations of hydrogen and He 1 4471 absorption 


its explanation in streams of gas which are known to distort velocity-curves of certain 
eclipsing variable stars. Pending further investigations of the form of the velocity-curve, 
we shall assume that we have a measure of the variation in velocity of the companion 
star, which we will classify only as early B, with K = 95 km/sec. This permits us to com- 
pute the following minimum masses: 


760, 
My sin? i = 19.60. 


Myr i 


In HD 193928 an absorption line near \ 3873 was observed, and evidence was pre- 
sented indicating that this line should be identified with He 1 3889. In HD 211853 a-simi- 
lar line was measured whenever detected on the spectograms. The results of the meas- 
urements are listed in Table 6. It is to be observed that the line was usually measured be- 
tween phases 4.8 and 6.6 days (on 13 of 16 spectrograms), but only once between phase 
0.0 and 4.8 days and then at phase 3.3 days. This indicates that the line may be variable 
with phase and is probably related to the general variations in the spectrum of the Wolf- 
Rayet star with phase. The average velocity of this line when identified as He 1 3889 is 
—1260 km/sec. 
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DISCUSSION 


In all the four Wolf-Rayet spectroscopic binaries thus far investigated by the author, 
the so-called ‘‘center of mass,” or gamma-velocity, has shown considerable divergence, 
from one emission line to another. The data thus far available on this question are col- 
lected in Table 7. Wilson measured three emission bands (He 11 4686, N v 4600, and N y 
4619) in HD 193576 and, upon accepting Edlén’s computed wave lengths for the two N y 
lines, as we have done in the present investigation, observed that 


V (4603) = V (4603) —V (4686) = +8 km/sec 
V (4619) = V (4619) —V (4686) = —19 km/sec . 
TABLE 7 
OBSERVED y-VELOCITIES 
Star Line y (Km/Sec) Remarks 

6, He 11 4686 emission +115 No violet absorption 
N v 4603 emission + 15 No violet absorption 
N v 4619 emission +165 No violet absorption 
Habs + 10 

He 11 4686 emission + 59 Violet absorption? 
N v 4603 emission —119 Weak violet absorption 
N tv 4058 emission —100 No violet absorption 

He t1 4686 emission + 50 No violet absorption 
N v 4603 emission + 62 No violet absorption 
N rv 4058 emission —125 No violet absorption 
Habs — 55 A 

He 11 4686 emission +137 No violet absorption 
N tv 4058 emission — 75 No violet absorption 
He 11 emission Pick- +180 Violet absorption 

ering 


He also observed that the gamma-velocity for the Wolf-Rayet star was more positive by 
90 km/sec than the gamma-velocity of the absorption-line component. The conclusion 
was made that the emission bands were shifted to the red by 90 km/sec. Struve!! con- 
cluded that in HD 151932, which is located in the cluster NGC 6231, N Iv 4058 gave the 
velocity of the cluster and that He 11 4686 possesses a red shift of 85 km/sec. The spec- 
trum of this star is remarkably similar to HD 214419. The emission bands are of similar 
width, and each star shows violet absorption edges of certain bands which are variable 
in intensity. In HD 214419, He 11 4686 is displaced to the red by 212 km/sec with respect 
to N 1v 4058. In this connection attention should be drawn to HD 211853, in which He 1 
4686 is again displaced to the red with respect to N Iv 4058; but neither band gives a 
velocity comparable to the absorption lines which probably belong to the B-type compo- 
nent. 

The existence of different gamma-velocities for various emission bands is certainly 
real and cannot be explained in terms of inaccurate laboratory wave lengths. It is diffi- 
cult to draw any general conclusions concerning the behavior of the various bands, other 
than that, thus far, He 1 4686 is always displaced to the red with respect to NV 1v 4058. 


11 Op. cit., p. 192. 
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Violet absorption edges do not, in general, appear adequate to explain all the displace- 
ments of emission bands with respect to each other. For example, in HD 186943, He 1 
4686, N v 4603, and NV v 4619 do not show any visible violet absorption, although the 
observed gamma-velocities differ by as much as 150 km/sec. This phenomenon of dis- 
crepant gamma-velocities is probably associated with the observed variations in profiles 
of certain emission bands, for which a nonsymmetrical shell about the Wolf-Rayet type 
component offers, at present, the most plausible explanation. 

In Figure 2 it is clearly shown that the velocity-curve as determined by the emission 
bands of N v differs, not only in amplitude, but also in phase from that of He 11 4686 by 
approximately 0.8 day. These differences are certainly real, and an explanation must be 
sought. To explain the differences in phase, the “transit-time-effect,” as developed by 
Wilson,” immediately suggests itself; but the assumption that the shell of the Wolf- 


TABLE 8 
MASSES OF WR STARS 


Star Type WR Star Absorption Star Ratio 
1522708 WC6+ M sin? i=1.9© M sin? i=6.9© 0.28 
WN5 M sin? i=5.8 © M sin? i=21.0 © .28 
Hi 193576... WN5 M=12.40 M=31.80 .39 
Hi} WN6+ M sin? i=7.6 © M sin? i=19.6 © 0.39 


Rayet star is not appreciably influenced by the more massive B-type companion does 
not appear to be justified. 

At the present time six Wolf-Rayet stars have been investigated as spectroscopic 
binaries. In four of the six stars, an absorption-line component is present; and in the 
other two, only the Wolf-Rayet component has been observed. This gives us some infor- 
mation concerning the masses of these peculiar objects. All the data known at present is 
summarized in Table 8.'* The ratio of the mass of a Wolf-Rayet star to that of an early- 
type absorption star appears to be fairly well defined, with a value of 3. Less can be said 
with certainty concerning the mass of a Wolf-Rayet star, but the evidence suggests that 
the mass of a Wolf-Rayet star is of the order of ten times that of the sun. With the com- 
pletion by the author of a photometric investigation of the three spectroscopic binaries 
discussed in this paper, more information should become available concerning the mass 
of a Wolf-Rayet star. 


The author wishes to express his thanks to Dr. O. Struve and Dr. S. Chandrasekhar 
for many helpful discussions. 
2 4p. J., 95, 409, 1942. 


18 HD 152270, O. Struve, Ap. J., 100, 387, 1944; HD 193576, O. C. Wilson, Ap. J., 95, 406, 1942; 
and HD 214419, W. A. Hiltner, Ap. J., 99, 274, 1944. 
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THE SPECTROSCOPIC ORBIT OF RZ ERIDANI* 


Cartos U. Cesco AND JORGE SAHADE 
Yerkes and McDonald Observatories 
Received March 29, 1945 


ABSTRACT 


The orbital elements of RZ Eridani, derived from 52 spectrograms, are as follows: P = 39.28074 days 
(photometric); K = 25.8 km/sec; w = 314.3°; e = 0.36; T = phase 28.0 days after minimum of light; 
y = +32.0km/sec;asini = 13 X 10°km; f(m) = 0.06 ©. The system is found to consist of a “‘metallic- 
line” main-sequence star, whose spectrum, in the region \\ 4000-4300, matches that of « Peg (F5 V), 
which at primary minimum is eclipsed by a subgiant G8 star surrounded by an envelope of calcium gas, 


The star BD—10°993 (8.0 mag.; a = 4°36™9, 6 = —10°57’5 [1855.0]) = HD 30050 
(Sp. A5) = AG Cbr M. 1174 (8.1 mag.), an eclipsing binary of the Algol type, was dis- 
covered as such in 1916 by Miss J. C. Mackie! on Harvard plates. Its provisional desig- 
nation was 45,1916. 


TABLE 1 


STAR LINES USED FOR THE DETERMINATION OF THE 
RADIAL VELOCITIES OF RZ ERIDANI 


During the months of August, September, October, and December, 1944, the star was 
included in the program of the 82-inch reflecting telescope of the McDonald Observatory 
for the purpose of determining its spectroscopic orbit. Fifty-two spectrograms were se- 
cured in that period, all of them on Eastman 103 a-O emulsion and with a dispersion of 40 
A/mm at \ 3933. Except at minimum light, the exposure time was of about 40 minutes 
in fair seeing. 

The spectrograms taken in the phases corresponding to maximum light show only one 
component; but on three plates taken during eclipse, namely, at phases 0.274, 0.381, and 
0.493 day, the spectrum is purely that of the secondary star. The determination of the 
spectral types of both components was very kindly made by Dr. W. W. Morgan, who ex- 
amined the spectrograms CQ 3847 (phase 11.212 days) and CQ 4080 (phase 0.493 day) 
(see Pl. XXXV). Dr. Morgan’s statement is as follows: 


CQ 3847.—In the range AA 4000-4300 the type is about F5; in this region the spectrum 
matches that of the main-sequence star « Peg (F5 V) fairly closely. The K line, however, is 
weak for this class. The ratio H/K corresponds approximately to spectral type A5?. The spec- 
trum probably belongs to the peculiar “‘metallic-line” group. The ratio \ 4045/\ 4077 indicates 
a higher luminosity than « Peg (Myis = +3.4). The classification given is uncertain because 
of unfamiliarity with plates of the dispersion used and the unavailability of suitable com- 
parison stars. 


* Contributions from the McDonald Observatory, University of Texas, No. 109. 


| Harvard Circ., No. 196, 1916; A.N., 207, 215, 1918. 


2 Probably the classification of the HDC for this star was made by taking into account the Cam K 
line. 
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TABLE 2 


RADIAL VELOCITIES OF RZ ERIDANI 


Radial 
Pate 
(in Km/Sec) 
Sept. 23 10:51 0.274 +44.9* 
Dec. 11 2:54 0.381 10.0* 
11} 5:36 0.493 23 .9* 
12 9:44 1.666 215 
er 12 10:17 1.688 15.6 
13 5:07 2.473 29.1 
ee 14 4:48 3.460 25:2 
/ Sept. 27 10:25 4.255 
28 9:33 5.219 32.6 
28 10:18 5.250 21.6 
Dec. 17 3:55 6.423 27.4 
cS Sept. 3 8:40 7.182 31.4 
30 9:18 7.209 26.1 
eee Oct. 1 8:40 8.168 25.1 
1 9:20 8.210 26.0 
y 10:40 9.266 30.0 
Dec. 20 2:52 9.379 16.0 
4 9:23 11,212 15.3 
< Aug. 26 11:34 11.584 
Ost 8:20 12.168 2.4 
5 8:58 12.195 27.3 
Aug. 27 12.583 29.8 
<r Oct. 6 8:35 13.179 13.1 
Ke 7 8:42 14.184 10.6 
7 10:49 14.272 
oo eee Aug. 30 11:16 15.571 1.4 
30 11:38 15.587 24.6 
31 11:37 16.586 14.0 
Oct. 10 8:22 17.170 19.0 
Sept. 1 10:44 17.549 15.3 
31 8:54 18.192 12.7 
| ee Sept. 2 10:43 18.548 17.6 
Oct. 12 10:44 19.268 9.3 
ee 13 9:28 20.216 14.4 
ee 14 10:10 21.245 24.0 
15 8:43 22.184 6.3 
16 9:09 23.202 
| 17 9:08 24.202 25.7 
18 9:42 25.225 37.9 
Sept. 10 11:02 26.562 45.3 
13 12:01 29.603 65.9 
co” 14 11:59 30.601 62.4 
| 15 11:30 31.581 57.0 
17 11:56 33.599 56.6 
19 11:58 35.601 56.0 
21 10:20 37.533 41.2 
1} | Dec. 9 5:06 37.753 34.3 
9 6:05 37.794 42.1 
ci eee Sept. 22 9:10 38.484 31.1 
22 9:52 38.513 444.7 


* The velocity is that of the fainter component of RZ Eridani, 
t The emission lines of Ca u give +30.2 km/sec. 
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372 CARLOS U. CESCO AND JORGE SAHADE 


CQ 4080.—The spectral type is probably a subgiant of class G8. The luminosity appears to 
be definitely fainter than that of a normal giant. The visual absolute magnitude is probably 
around +3. The same remarks apply to the accuracy of the classification. 


The spectrum of the secondary component shows Ca 1, H and K, emission lines. 

In Table 1 are listed the lines measured for radial velocities, and in Table 2 the results 
of our measurements, which are plotted in Figure 1. In this figure the crosses represent 
the values derived for the secondary star. Of the three plates obtained during the eclipse, 
the best is CQ 4080; the other two plates are rather poor and do not go beyond ) 4000. 
The radial velocity derived from the emission lines on spectrogram CQ 4080 is almost the 
same as the one derived from the absorption lines. Therefore, we conclude that the G8 
component is surrounded by an envelope of calcium gas. 
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In plotting the measurements we have adopted as the period the value given by 
Gadomski,* which was a mean between his results and those of Hertzsprung, namely, 
39.28074 days. The phases were computed by considering as the time of minimum the 
date JD 2423854.057, given by Schneller in his Katalog und Ephemeriden Verdnderlicher 
Sterne fiir 1941. 

The orbital elements were derived by applying the method of Lehmann-Filhés, and 
the results are shown in Table 3. No attempt was made to perform a least-squares solu- 


TABLE 3 
ORBITAL ELEMENTS OF RZ ERIDANI 


P = 39.28074 days (assumed) vy = + 32.0 km/sec 
K = 25.8 km /sec asint= 13 X 10°km 
o= 314.9 my 0.066 

e= 0.36 aia 

T = phase 28.0 days 


tion because of the scattering of the measurements. Before this computation was made, 
we plotted, tentatively, the measurements derived only from the H lines and those de- 
rived from the metallic lines; but, since the trend of the curves was quite similar, we have 
considered the straight means of the velocities derived from all the lines, in the determi- 
aation of the elements. 


3 Circ. Obs. Cracovie, No. 24, 1927. 
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The time of mid-eclipse, computed from the condition v + w = 90°, is phase 39.0 
days, a value which is in good agreement with the photometric elements. 
The only light-curve of RZ Eridani which has ever been published, so far as we know, 
is the one by Gadomski,* who found: 


Gadomski’s observations seem to indicate’ the probable existence of a secondary mini- 
mum of amplitude 0.2 mag., which is advanced by 0.4 day. 

The eccentricity of 0.36 which we have obtained seems to be real, as the spectrum of 
the primary star does not show any change in the different phases; hence, the possibility 
of a distortion of the velocity-curve was rejected. Furthermore, in a system with such a 
long period as that of RZ Eridani we would expect an eccentricity of the order which we 
have found. It is to be hoped that new photometric observations of this star will give more 
exact information about the secondary minimum and will provide a check of our spectro- 
scopic observations. 


We are indebted to Dr. C. Payne-Gaposchkin for suggesting the study of this star, to 
Dr. O. Struve for taking some of the spectrograms and for discussion of the material, and 
to Dr. W. W. Morgan for helping us in classifying the spectra. 

4 Bull. Acad. polonaise d. Sci. et d. Lettres, Ser. A, p. 303, 1926. 

5 Tbid., p. 310. 
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NOTES 


SERIES LINES OF MAGNESIUM IN THE SOLAR SPECTRUM 


In an earlier note! our knowledge of two series of Mg lines (3'D — n'F° and 3'D — 
n°F°) previously recognized in the laboratory by Paschen? was extended and improved 
with the aid of solar data. 

The purpose of the present note is to correct two errors which occurred in connection 
with = 5, to which we formerly assigned solar lines \ 9257.01 and A 10813.14 in the 
singlet and triplet series, respectively. Better solar spectrograms now fail to show 
X 10813.14, which had been observed only on a plate on which other weak lines were sub- 
sequently found to be spurious. In place of these early assignments we now introduce for 
n = 5 two other solar lines, \ 9255.79 and \ 10811.14, of more appropriate intensity. 
We have also revised the solar wave lengths, intensities, term values, n'*F°, Rydberg 
denominators, ”*, and series limits, as shown in Table 1. 

Our erroneous assignments resulted from an attempt to fit the solar series to Paschen’s 
laboratory wave lengths, \ 9257.9 and A 10812.9; but experience now shows that his data 
for Mg require larger and more diverse corrections than we at first supposed. For this 
reason and since less than half of the lines in Table 1 have been found in the laboratory, 
we base our identifications on improved series relations. 

The identity of the 'F° and *F° terms is now established beyond question, and two 
new members are added to the singlet series. 

The limits of the series, 3'D and 3*D, shown at the top of Table 1, have been derived 
by adding 1.5 cm and 1.37 cm“, respectively, to Paschen’s values. As we state the lim- 
its, their difference is equal to 1553.97, the weighted mean value of Av. In obtaining the 
mean Av we have excluded the first and have given weight 1 to the next four, and weight 
} to the last four, values. 

Wave lengths of the two new singlet lines are predicted by assuming that n — n* 
converges to the value +0.0478, as indicated by the data for m < 12. With this assump- 
tion n* takes the values 12.9522 and 13.9522; and the corresponding term values are 
computed from the relation n'*F° = Ra+ n**, where Ro is the Rydberg constant for a 
nucleus of infinite mass. We take Ro = 109737.303 cm, as given by Birge,* and find 
the terms to be 654.13 and 563.73, respectively; their combinations with the series limits 
give predicted wave lengths \ 6841.17 and A 6799.10. Since the agreement of the ob- 
served lines shown in Table 1 is satisfactory, their terms and Rydberg denominators are 
derived and tabulated. 

The best lines of these series are difficult to measure; the weakest are at the limit of 
visibility. It is, therefore, notable that between \ 6900 and A 10,800 the accordance 
among the values of Av corresponds to errors in the measurements no greater than 0.02 
or 0.03 A. Since the measured intervals, Av, are independent and overlap in wave length, 
it is probable that throughout this region our scale of wave lengths is free of systematic 
errors large enough to affect the identification of solar lines. 

Estimates of intensity for wide diffuse lines are generally difficult—the more so when, 
as for these Mg lines, a great range of wave length and of intensity is involved. Beyond 


! Russell, Babcock, and Moore, Phys. Rev., 46, 826, 1934. 

2 Sitz.d. Preuss. Akad. d. Wiss., Phys.-math. Kl., 32,709, 1931. 

3 Reports on Progress in Physics (London Physical Society), 8, 129, 1941, 
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d 9000.our estimates in Table 1 may be systematically too low; at shorter wave lengths 
an error of two units may readily occur. In spot spectra the Zeeman effect further in- 
creases the uncertainty. The data show, nevertheless, a tendency of these lines to weaken 


TABLE 1 
INFRARED SERIES OF Mgt IN THE SUN 


3'!D2=15267.50 
(EP 5.73 VotrTs) 


33D3,2,1 = 13713.53 
(EP 5.92 Votts) 


RESULTS 


Int. Int. 
—n'F° 33D —n3F° n* 
Disk Spot Disk Spot 

8273.28 6719.9 (1553.4) | 6994.22 | 3.9610 
pee 9255.79 4N | 3 10811.14t 5N| 4N 

10801 .09 9247.18 1553.91 | 4466.38 | 4.9568 
heed 8213.041 4N | 3N | 9414.95 5N| 5N 

12172.42 10618 .49 1553.93 | 3095.06 | 5.9545 
7691. 569 9N | 9N | 8736.040 10N| 5N 

12997 . 68 11443. 69 1553.99-| 2269.83 | 6.9531 
7387 .700 5N | 5N | 8346.131 6N| 2 

13532.29 11978. 31 1553.98 | 1735.22 | 7.9524 
7193. 183 1N | | 8098.746 5N| 5N 

13898 .23 12344. 20 1554.03 | 1369.30 | 8.9522 
7060.446§||| 4N | 4 7930.819 2N| 2N 

14159.52 12605 .59 1553.93 | 1107.96 | 9.9521 
ees 6965.408|| | —1N | ob 7811.16 1N| ob? 

14352.71 12798 .68 1554.03 | 914.82 | 10.9524 
6894.89|| | —3N |...... 7722.64 —1N| —2N 

14499. 50 12945 .38 1554.12 | 768.07 | 11.9530 


t Uncorrected laboratory value by Paschen. 
t Not fully resolved from strong solar line \ 10810.89. 
§ Blend with atmospheric line. 
|| Spot effects for five weakest singlets are confused by molecular absorption. 
§ ob = obliterated. 
** Blend with weak solar line. 


in spots. The apparent anomalies shown in the progression of intensity along these series 


are possibly analogous to the abnormal weakening noted by Fowler‘ in \ 6965 of K and 


\ 2722 of Ca. 


4 Report on Series in Line Spectra, pp. 102, 124, London, 1922. 
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376 NOTES 


The Photometric Atlas of the Solar Spectrum by Minnaert, Mulders, and Houtgast® 
shows that in the disk equivalent widths progress along the series, though none too 
smoothly. A calibration of estimated intensities for solar lines beyond \ 6600 in terms 
of equivalent width is in progress. 

D. Bascock 
CHARLOTTE E. Moore 
Mount WILSON OBSERVATORY 


PRINCETON UNIVERSITY OBSERVATORY 
February 1945 


NOTE ON THE FUTURE ORBIT OF COMET DELAVAN (1914V) 


Following E. Strémgren’s work! on the original orbit of near-parabolic comets, several 
investigators have carried perturbations backward in similar and more recent cases, al- 
ways finding that, when such objects are followed far enough back from the time of 
perihelion, the hyperbolic excess of the eccentricity invariably becomes negative. The 
obvious deduction that comets are permanent members of the solar system has, how- 
ever, been weakened by the case of an orbit which, though originally elliptic, was ‘found 
to extend to nearly stellar distances at aphelion. A typical example is that of Comet 
Delavan 1914V, where the orbit was determined within narrow limits.” 

Chandrasekhar? has called attention to the interest of such a case in connection with 
the question of the stability of binary systems and to the desirability of integrating the 
cometary motion forward so as to find whether the orbit reverts in the future to an eccen- 
tricity as nearly parabolic as when the comet approached the outskirts of the solar sys- 
tem. Starting from the slightly hyperbolic elements? deduced from the discussion of all 
the material available on this comet, I have computed the co-ordinates over a period of 
nearly twenty years, taking into account the perturbations due to the four major planets. 
From the values obtained for 1933, June 17 and November 14, separated by an interval 
of 160 days, I found, after reduction to the center of gravity of the solar system, the 
following dynamical elements: 


e = 0.999861 and 10.000126. 


The comet was at that time at a distance of 39 A.U. from the sun, and this circumstance 
reduces to a maximum of 8 per cent the further possible changes in the major axis. The 
elliptical character of the orbit is here much more pronounced than at the corresponding 
interval before perihelion. The aphelion distance now becomes about 16,000 A.U., as 
against 170,000 A.U. found for the same interval of time before perihelion. 


G. VAN BIESBROECK 


YERKES OBSERVATORY 
February 28, 1945 
5 Amsterdam, 1940. 
1 Pub. Kobenhavn, No. 19, 1914. 
2 Pub. Yerkes Obs., V, Part II, 31, 1927. Ap. J., 99, 58, 1944, 
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NOTE ON THE ECLIPSING BINARY SX HYDRAE 


Velocities of SX Hydrae have been determined from 14 spectrograms taken by Dr. O. 
Struve at the Cassegrain focus of the 82-inch reflector of the McDonald Observatory. 
The star is unusual for the great depth of primary minimum, the light dropping from 
My, = 8.6 at maximum to mpg = 12.6 during total eclipse. The photometric orbit has 
been computed by Mrs. M. B. Shapley,' who finds a period of 2.895697 days, semidura- 
tion of primary eclipse 0.171 day, semiduration of totality 0.021 day (uniform solution), 
and no detectable secondary minimum. Mrs. Shapley gives radii of 0.104 and 0.298 and 
an inclination of 79.1°. The relative intensities in photographic light are 0.965 and 
0.035. 

From 12 spectrograms taken at the Mount Wilson Observatory during 1934-1936, 
R. F. Sanford? has derived the approximate spectroscopic elements. Table 1 gives the 


TABLE 1 
RADIAL VELOCITIES OF SX HYDRAE 


R.V. (Ka/SEc) 
PLATE DATE PHASE 
€ O-C 
1944 

May 2 5h4gm 65° — 100 —74 —26 
79 81 —79 —18 
8 6 11 93 81 —86 + 5 
3231 14 5 50 117 80 —90 +10 
5 ee Apr 30 3 39 165 62 —72 +10 
i: ee 30 4 41 170 54 —68 +14 
(7 30 7 07 183 69 —54 —15 
1S) ae May 3 6 18 192 68 —44 —24 
‘i ee 5 44 202 26 —32 + 6 
L.A 9 6 15 218 26 —17 — 9 
15 6 04 242 16 0 —16 
1 4 07 292 6 + 5 —11 
| ae 1 5/25 299 13 + 2 —15 
7 6 22 329 — 29 —13 —16 


velocities from the McDonald spectrograms under the heading ‘“‘O” and the velocities 
from Sanford’s velocity-curve, at corresponding phases, under the heading ‘C.” The 
phases are counted from primary minimum. In consideration of the wide and diffuse 
character of the lines in the spectrum of the primary, which is of HD type A3, the agree- 
ment is perhaps as close as might be expected. The McDonald spectrograms indicate a 
slightly smaller semiamplitude A than Sanford has given, but the evidence is inconclu- 
sive. The dispersion on the McDonald plates is 55 A/mm at Hy, as compared with 75 
A/mm on Sanford’s plates; however, the scatter of the velocities ““O” in the table is no 
less than the scatter of Sanford’s velocities, so that the computation of a new orbit has 
been deemed unwarranted. 

Dr. Struve has also secured several spectrograms during eclipse, using a Schmidt f/2 
camera with glass prisths. This combination gives a dispersion of about 75 A/mm at 
Hy. Several of the films taken during the partial phase have been measured for radial 


1 Harvard Bull., No. 797, 1924. 
2 Ap. J., 86, 153, 1937; Mt. W. Contr., No. 574. 


= 
st 
* 
] 
a 
of | | | | 
V- 
id 
th 
1e 
n- 
ill 
f 
| 
e 
| | | | ae 
> 
— 


378 NOTES 


velocity, but the scatter is too great for conclusions to be drawn regarding rotation. The 
character of the lines appears not to change with the phase. On one spectrogram, at 
phase 2°, the spectrum of the secondary component appears weakly. From comparison 
of the region of the D lines with the same region on a comparable spectrogram of 8 UMi, 
it is concluded that the type of the secondary is near K5. This result is in approximate 
agreement with the type implied by Mrs. Shapley’s ratio of surface brightnesses (220 
for the uniform solution). If the mass of the A3 star is 3©, Sanford’s mass function, to- 
gether with Mrs. Shapley’s inclination and radii, gives the unusual ratio of mean densi- 
ties of 0.0018. 


I wish to thank Dr. Struve for putting at my disposal the spectrograms of SX Hyd 
and Dr. Hiltner for securing several spectrograms of 8 UMi. 
ARMIN J. DEUTSCH 


YERKES OBSERVATORY 
April 10, 1945 
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REVIEWS 


“The Arc Spectrum of Iron (Fer): Part I, Analysis of the Spectrum,” by HENRY Norris 
RUSSELL and CHARLOTTE E. Moore, and “Part Il, The Zeeman Effect,” by Dorotuy W. 
WEEKS. Transactions of the American Philosophical Society, 34, Part I1, 111-207, 1944. $2.25. 


All astrophysicists interested in the identification of stellar lines or in the study of curves 
of growth will welcome this outstanding contribution concerning the atom that reveals 
the largest number of absorption lines in the sun and in most stars of spectral classes from F to 
M. Moreover, the paper will also be beneficial to astrophysics from a very different point of view. 
There are still many atomic spectra of astrophysical interest that await detailed analysis—as, 
for example, most of the doubly ionized metals and rare earths. There are good reasons to be- 
lieve that elements such as Cr 111, Ni 11, and Ti 1 will eventually be found to play as important 
a role in stellar spectra as Fe 11. But in recent years very few physicists have been willing to 
tackle this work on analysis, and newcomers in the field will have to be found either among 
the astronomers themselves or, preferably, among young students of spectroscopy. Analyzing a 
spectrum requires a good knowledge not only of theoretical spectroscopy but also of a number of 
semiempirical rules (dare I say tricks? Spectrum analysis has often been likened to an exalted 
crossword puzzle!) that do not appear clearly in the textbooks. In earlier publications Dr. 
Russell and his pupils have often mentioned helpful devices for estimating the behavior of 
certain terms, and other spectroscopists have also done their best. But, so far as I know, no 
example has ever been able to provide such a complete and clear. exposition of the technique 
of analysis of a spectrum as does Fe 1 in the hands of Russell, Moore, and Weeks. Most certainly, 
many newcomers in the field will, in the years to come, turn to this paper for inspiration and 
help. 

It is quite remarkable that a spectrum which has been observed for more than seventy years 
and whose first twenty multiplets were found as early as 1923 could still reveal 60 new terms 
to the authors in 1944. But the spectrum of Fe 1, although quite orderly, is indeed a very com- 
plicated one. Moreover, the investigation has been delayed by political conditions, since a large 
number of terms (89 unpublished) found by Catalan—who in peacetime would have been a 
joint author, as Russell and Moore emphasize—could be communicated only after long delays. 
Additional progress will now require new observations. Laboratory sources may possibly be 
found that reveal more Fe 1 lines than the sources used thus far; this is apparent from the fact 
that the sun exhibits many more Fe I lines than any iron arc or spark. Moreover, the analysis 
is based on an approximate wave-number range of only 45,000 cm~ (from \ 12000 to \ 1850). 
The spectral region covered may still be extended by about 10,000 cm in the ultraviolet and 
by some 8000 cm~ in the infrared, and such an extension may reveal a few new terms; the ex- 
tension will be especially important in the infrared. But, on the whole, the Fe 1 spectrum is just 
about as completely analyzed as it can be or as its companions, Fe 1 and Fe Im, are. 

The authors list 4860 classified lines arising from combinations among 464 energy-levels. 
The wave lengths are the best available; yet the material is not homogeneous. While certain 
wave lengths have been determined very accurately for use as International secondary stand- 
ards, a few are known only to 0.1 A. The situation is even worse for the intensities which are still 
in an extremely unsatisfactory state. It is to be hoped that the monograph by Russell, Moore, 
and Weeks will induce laboratory spectroscopists to work toward a more critical list of wave 
lengths, and especially toward a set of homogeneous intensities, extending the determination of 
f-values by R. B. and A. S. King. 

Several of the latest additions to the analysis of Fe 1 were obtained by comparison with the 
previously published analyses of Fe m1 and Fe 11 or with the help of the Zeeman spectrograms. 
The principal low terms of Fe 1 arise from the configurations 3d*4s? and 3d’4s and are thoroughly 
intermingled. The singlet terms belonging to 3d°4s? were found by comparison with Fe 11, in 
which the terms corresponding to 3d° lie much lower than all others and have been identified 
with certainty. The comparison between d? of Fe m and d’s of Fe1 was similarly instructive. 
The assignment of certain low terms to a definite configuration was based on a comparison with 
other atoms in which the terms based on d®, d®~'s, and d®~*s? are known. All along for the medium 
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and high terms, the analysis proceeds by frequent comparisons with other atoms having similar 
electron configurations. Obviously, the Fe1 analysis contains, now, a wealth of information 
that will help in analyzing other spectra later on! Very few strong lines of Fe 1 observed in 
the laboratory remain unclassified, but a good many lines of low intensity or of poor wave 
lengths are still unclassified. The new value of the ionization potential is 7.858 volts, which is 
0.02 volt higher than the previously accepted value and is probably reliable to 0.01 volt. At 
the end of page 113 the authors make the erroneous statement that Fe 1 is, up to the present, 
the only spectrum in which terms of four different multiplicities are known. Both Cr1 and 
Fe 111 now share this privileged status with Fe 1, which first came to possess this position when 
Catalin found the singlet terms. 

A long section of the paper deals with the problem of the predicted iron lines in the solar 
spectrum. In 1928 Russell and Moore had discovered that many Fe 1 lines, predicted from the 
term values but not yet observed in laboratory sources, were actually present in the sun. On 
the basis of this example, it has now become common practice in stellar identification work to 
consider certain predicted permitted lines of abundant metals as well as the laboratory lines. The 
very complete data on Fe I, now available, have enabled Russell and Moore to make a thorough 
statistical discussion of the real and accidental coincidences of predicted and solar lines. Their 
discussion sets an example that should be followed in all problematical identification cases, 
especially when molecular bands are involved. The evidence obtained is quite decisive: there 
should be some 1200 lines of iron which, though not yet produced in the laboratory, should be 
observable in the solar spectrum between A 2975 and A 6600. A table is given of those predicted 
lines for which the evidence of presence in the region from A 2980 to A 11000 was adjudged to be 
good or fair—1254 in all. Doubtless, some accidental coincidences are present, and some real 
iron lines have been omitted in the rather drastic process of exclusion. Moreover, many lines 
(some 300 between 2975 and 6600) have perforce been omitted, owing to masking in the 
sun. But, with these limitations in mind, any astrophysicist will often be greatly helped by the 
table of predicted lines present in the sun, since the great majority of the tabular lines must be 
real. 

The Zeeman spectrograms, obtained with the Bitter electromagnet at the Massachusetts 
Institute of Technology and studied by Miss Weeks, have been of great help in the analysis. 
The field strengths ranged from 83,000 to 87,000 oersteds; and the plates were measured upon 
Harrison’s automatic measuring machine, which recorded on a film the wave lengths and rela- 
tive photographic intensities of the components. Miss Weeks gives a table of the observed Zee- 
man patterns for 1038 lines and the g-values deduced from each, and another table containing 
the g-values for all levels for which they have been determined. A detailed comparison of ob- 
served and theoretical g-values is included. Deviations from theory, because of perturbations, 
are discussed. As is well known, perturbations may be anticipated between neighboring levels 
of the same parity and with the same J, irrespective of multiplicity, Z-value, and electron con- 
figuration. In the absence of Zeeman data, perturbations are indicated by the mutual re- 
pulsion of the levels involved (interrupting the usual sequence of components of these terms) 
and by strong intersystem combinations. Similarly, the g-values are also conclusively affected by 
perturbations. 

The authors must be congratulated on having combined in this memoir an abundant source 
of data precious to astrophysicists and a general sketch of the analysis procedure which will be 
an inspiration to newcomers in this field. The American Philosophical Society has spared no 
effort in making the monograph a beautiful piece of printing; indeed, the tables are very easy to 
read, compared with most spectroscopic papers. Only four minor misprints have been found in 
the course of work based on the tables: the word ‘‘components” is misprinted in the table head- 
ings on pages 189 and 201, and the two wave lengths A 3820.572 (p. 155) and \ 4638.687 (p. 175) 
should be replaced by \ 3827.572 and \ 4628.687, repectively. 

P. Swincs 


Pasadena, California 


Stellar Variability, Vol. 11: Eclipsing Variables. By D. J. MArtiNnorrF. (In Russian.) Moscow and 
Leningrad: ONTI, 1939. Pp. 150. R. 4.25. 
Although this book was originally written in 1935 and 1936 and was revised in 1939, it has 
only recently arrived in this country. The author intended it to form the second volume in a 
series of monographs on variable stars by B. W. Kukarkin, P. P. Parenago, W. P. Zessewitsch, and 
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D. J. Martinoff. The first volume, by Kukarkin and Parenago, was devoted to the physical vari- 
ables and was published in 1937. The third volume will treat a number of practical questions en- 
countered in the study of variable stars. Martinoff’s book is an introduction to the study of 
eclipsing variables. In the first three chapters it presents the theory of orbits of eclipsing stars 
and a number of new diagrams to illustrate the influence of the principal elements upon the 
light-curve. The treatment follows the classical theory by Russell. It does not reproduce the 
various tables required in actual computations (these will presumably be given in the third vol- 
ume of the series) but attempts to clarify the essential points upon which the theory depends. Al- 
though much recent work in this field has not been included and the discussion is not complete 
in all respects, the book constitutes an unusually lucid account of the more elementary aspects 
of the theory and is, therefore, especially valuable as a textbook for students. Chapter iv deals 
with periodic and secular changes in the periods of the eclipsing variables. Chapter v, on the 
physical characteristics of eclipsing variables, contains some new material, as well as material 
published by Russian astronomers in recent years, not all of which has been available to readers 
who are unfamiliar with the Russian language. The last chapter, on the origin and evolution of 
close binaries, is quite short and limited in scope. 

Although the Russian language constitutes a barrier to the effective use of this book outside 
of Russia, it must be recognized that there is nothing even approaching it in the English lan- 
guage. The original articles written by Russell and by Russell and Shapley, published in 1912, in 
Volumes 35 and 36 of this Journal, are now out of print. There is a condensed account of the 
simplest case of Russell’s theory in Aitken’s book, The Binary Stars (2d ed.; New York: Mc- 
Graw-Hill Book Co., 1935). But there are no textbooks which treat this theory in such a manner 
as to enable the student to start practical work on eclipsing variables without an enormous 
waste in time and many sad experiences. It is important that we should realize our lack of ad- 
vanced textbooks in astrophysical subjects. The astrophysical monographs sponsored by the 
Astrophysical Journal are intended to relieve this need; but there is, as yet, not enough support 
for this venture on the part of readers and potential authors. The principal difficulty consists in 
the fact that our universities and research institutions feel no responsibility in the matter of 
textbooks, while commercial publishers know that they cannot make a profit from advanced . 
books and therefore prefer to publish elementary texts. Since 1936, when the late B. P. Gerasi- 
movié edited two volumes written by a number of leading Russian astronomers,! B. A. Voront- 
zov-Velyaminav has published a Textbook of Practical Astrophysics,? 648 pages in length, and a 
Collection of Problems in Astronomy,’ comprising nearly 200 pages—a total of 1,246 questions 
and their answers, together with tables, reproductions of stellar spectra, etc. Not available to me 
are a Textbook of Stellar Astronomy, by P. P. Parenago, and a book on Variable Stars and Methods 
for Observing Them, by Parenago and Kukarkin. To the extent that future progress in astronomy 
depends upon the training of advanced students, the Russians will have a marked advantage 
over other nations because of their policy to produce the necessary books. Perhaps their success 
in this respect will be tempered by the careless external appearance oi the books: the paper is 
very poor, the illustrations of photographs are almost useless, and even the writing is sometimes 
defective (though not in the volume under review). 

O. STRUVE 


Yerkes Observatory 


The Mathematics of Physics and Chemistry. By HENRY MARGENAU and GEORGE MOSELEY 
Murpuy. New York: D. Van Nostrand Co., 1943. Pp. xii+581. $6.50. 


The subject matter of this book covers an unusually wide range of topics. In fact, the au- 
thors’ intentions have been to present a fairly complete account of the basic mathematical tech- 
niques in modern physics and chemistry. The exposition follows a moderate level of mathemati- 
cal rigor with an occasional omission of proofs. A number of problems are included. 

An account is given of ordinary differential equations, with frequent applications to interest- 
ing physical problems. Special emphasis is given to second-order equations, and particularly to 
the common special equations in physics, together with the elementary properties of their solu- 


1 Reviewed in this Journal, 86, 106, 1937. 
2 Moscow and Leningrad: ONTI, 1940. 
3 Moscow and Leningrad: ONTI, 1939. 
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tions. The classical partial differential equations likewise find a place in the text with applica- 
tions to eigenvalue problems. Good elementary treatments of vector analysis, matrices, integral 
equations, and group theory are given. One chapter is devoted to numerical calculations involv- 
ing interpolation, numerical integration, the theory of errors, and other questions. 

In separate chapters are included topics with more specific reference to definite physical sub- 
jects. They are: (1) “Thermodynamics,” (9) ““Mechanics of Molecules,” (11) “Quantum Me- 
chanics,” and (12) ‘Statistical Mechanics.”’ In these chapters some physical background is 
given, in addition to the purely mathematical subject matter. The most extensive of these is the 
one devoted to quantum theory. Here, following a postulatory introduction, a large variety of 
the fundamental problems find a treatment. 

The book is remarkable in the high level of clarity and general scientific quality attained in 
the treatment of so extensive a range of topics. One should note that this range includes much 
that is not found in most texts of comparable purpose and level. For these reasons it should prove 
a useful reference volume, in addition to its primary purpose as an introductory textbook. How- 
ever, one may note that, precisely because of the book’s conciseness, the student may have diffi- 
culty in appreciating the physical aspects of the subjects without a good deal of additional ma- 
terial in the form of supplementary reading or lectures, or without prior knowledge. 


L. G. HENYEY 
Yerkes Observatory 


Tables for Converting Rectangular to Polar Co-ordinates. By J. C. P. MILLER. New York: Dover 
Publications, 1939. Pp. 16. $0.75. 
This useful table was published in 1939 and is recommended for such problems as the trans- 
formation of harmonic constants a and 6, obtained by harmonic analysis, to amplitude ¢ and 
phase angle €, in accordance with the relation 


c sin (nt + €) = asin bcosat, 
or the evaluation of the magnitude and direction of a vector from rectangular components and 


the conversion of complex numbers from the form x + iy to the form rei#. The tables are intend- 
ed for use with a calculating machine or a slide rule. 
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